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SUMMARY 

This  report  documents  the  findings  of  a  3  year  study  to  asses  the  impact  of  Large  Woody 
Debrio  (LWD)  upon  channel  evolution  and  morphology  in  unstable  sand  bed  rivers  in  northern 
Mississippi,  The  aim  of  this  research  is  to  gain  an  improved  understanding  of  the  basin-wide 
impact  of  LWD  dynamicft  in  unstable  and  stable  channel  environments  and  to  develop  a  set  of 
coherent  debris  management  strategies  for  erosion  control,  habitat  enhancement,  and 
maintenance/design  considerations  for  run-of-river  structures,  based  upon  sound  geomorphic 
and  engineering  analysis. 

Data  from  the  US  Army  Corps  of  Engineers  Demonstration  Erosion  Control  (DEC)  survey 
program,  conducted  in  May  1995,  was  used  to  locate  significant  debris  jams  with  respect  to 
planform  and  long  profile  data  on  23  river  reaches  in  northern  Mississippi.  The  reaches 
surveyed  are  between  4000  and  12000  feet  long  and  range  in  upstream  basin  area  from  3.5  to 
150  square  miles.  A  comprehensive  understanding  of  debris  dynamics  can  be  established  by 
surveying  these  channels  because  the  reaches  fell  into  several  categories  including, 
stable/unstable  reaches,  straight/meandering  reaches  and  reaches  which  have  either  a 
predominantly  agricultural  or  wooded  riparian  zone.  The  debris  jams  in  each  reach  were 
surveyed  in  detail  to  determine  the  mechanisms  and  locations  of  debris  input,  jam  impact  upon 
channel  morphology  and  sediment  routing  and  jam  stability  over  time,  The  last  of  these 
objectives  was  assessed  by  comparing  the  survey  results  of  the  current  study  with  those 
obtained  in  a  1994/1995  research  project  (see  Wallerstein  &  Thorne,  1995), 

An  up  to  date  review  of  the  literature  concerning  the  geomorphic  impacts  of  in  channel 
LWD,  and  current  LWD  management  strategies  is  presented. 

A  conceptual  debris  jam  classification  mode!  is  presented,  b.'  sed  upon  ‘key’  debris  dimensions 
relative  to  channel  width  and  found  to  correspond  well  with  the  field  data,  Thif  .odel  also 
explains  the  spatial  distri’mtion  of  sedimentation  and  scour  processes  associated  with  debris 
jams.  Results  show  that  the  net  balance  between  debris  induced  sediment  retention  and  debris 
induced  channel  scour  is  in  favour  of  sedimentation,  indicating  that  debris  jams  help  to 
accelerate  sedimentation  processes  in  these  channel  systems.  Comparison  of  survey  results 
between  summer  1995  and  summer  1996  shows  there  to  be  little  change  in  the  number  and 
position  of  debris  jams  present  in  each  reach  It  appears,  therefore,  that  debris  jams  arc  stable 
features  in  the  short-term,  although  a  better  understanding  of  jam  longevity  can  only  be 
achieved  through  a  long-term  monitoring  program. 
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The  LWD  Management  Program  (DMP)  has  been  tested  and  improved  based  upon  the  new 
findings.  The  Debris  at  Bridge  Pier  Prediction  Program  (DBP3)  has  also  been  improved  and 
now  calculates  the  pier  scour  adjustme.it  factors,  whereas  these  had  to  be  entered  by  hand  in 
the  previous  version. 

Both  the  DMP  and  DBP3  programs,  and  user  manuals  (written  in  MS  Word  ver.  7)  are 
enclosed,  inside  the  rear  cover  of  this  report, 

The  GIS  user  interface  for  the  DMP  is  also  now  available  for  use  and  can  be  obtained  upon 
request  from  nick@geography, nottingham.ac.uk 

The  program  runs  on  the  UNIX  version  of  ARC  INFO  and  requires  approximately  98  mb  of 
memory  to  install. 

Final  conclu.sions  regarding  the  management  of  LWD  in  unstable  channel  systems  are  made 
based  upon  the  results  of  this  three  year  research  program. 


UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  report  can  be  converted  as  follows: 


To  convert 

To 

Multiply  by 

inches  (h) 

millimetres  (mm) 

25.4 

feet  (ft) 

meters  (in) 

0.305 

yards  (yd) 

meters  (m) 

0.914 

miles  (Mi) 

kilometres  (km) 

1.61 

square  miles  (sq.  miles) 

square  kilometres  (km  ') 

2.59 

cubic  feet  per  second  (cfs) 

cubic  meters  per  second  (cms) 

0.0283 
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1  INTRODUCTION 

There  has  been  increasing  interest  in  the  role  of  vegetation  in  fluvial  geomorphology  in  recent 
years  because  it  has  been  recognised  that  river  dynamics  cannot  be  fiilly  understood  without 
taking  into  account  the  impact  that  vegetation  has  upon  bank  stability,  flow  velocity,  and 
riverine  habitat. 

As  a  consequence,  the  study  of  in-channel  Large  Woody  Debris  (LWD)  or  Coarse  Woody 
Debris  (CWD)  as  it  is  sometimes  referred  to  (that  is  trees,  branches  and  other  larger  organic 
matter,  operationally  defined  as  material  with  a  length  greater  than  1  metre)  and  its 
accumulation  as  jams  or  dams  and  impact  upon  the  channel  environment  have  become  topics 
receiving  increasing  research  interest  over  the  past  5  to  1 0  years. 

in  a  review  of  relevant  literature  undertaken  in  an  earlier  study  (Wallerstein  1994)  it  was 
established  that  a  large  proportion  of  the  research  performed  to  date  has  been  carried  out  in 
upl.iid  areas,  and  in  stable,  gravel  bed  rivers  such  as  in  the  Pacific  North  West  (Hogan  ct.  al, 
1995,  Fetherston  ct.  al.;  1995)  to  determine  the  impact  of  LWD  on  salmon  habitat  and 
migration,  and  in  relation  to  logging  operations  and  forest  management.  Very  little  is  known 
about  the  impact  of  LWD  in  unstable  or  sand  bed  rivers.  Much  of  the  work  is  fairly  qualitative 
and  observational  in  nature  and  there  has  been  little  emphasis  on  determining  the  key  variables 
at  play  in  LWD  dynamics,  and  the  modes  of  their  interaction.  Most  studies  have  also  been 
undertaken  in  isolated  reaches,  rather  than  covering  basin-wide  debris  processes,  although 
there  are  one  or  two  notable  exceptions  (see  Gregory  et  al.,  1993). 

At  present,  LWD  management  is,  therefore,  conducted  from  an  incomplete  understanding  of 
debris  impacts  and  dynamics  and  operational  maintenance  of  debris  is  carried  out  on  ati  aJ  hoc 
basis. 

The  aim  of  this  research  etYort  has  been  to  asses  the  catchment  wide  impact  of  LWD  over  a 
wide  range  of  channel  sizes  in  unstable,  rapidly  evolving  rivers  with  sand,  clay  and  loess  bed 
and  banks.  The  research  in  this  project  ha.s  been  centred  on  streams  in  the  DEC 
(Demonstration  Erosion  Control)  watersheds  draining  the  Bluff  Line  hills  of  Northern 
Mississippi,  which  are  known  to  be  evolving  rapidly  in  response  to  complex  response  in  the 
fluvial  system  following  catchment  land  use  changes  and  past  engineering  interventions. 

The  specific  aims  of  this  research  are; 

1)  To  collect  a  large,  representative  data  set  concerning  the  reach  scale  and  basin-wide 
influerces  of  I  ,WD  on  channel  morphology  in  a  different  type  of  channel  environment  to  that 


which  has  been  studied  to  date;  specifically,  in  unstable,  rapidly  evolving  sand-bed 

rivers. 

2)  To  assess  whether  there  are  preferential  sites  of  debris  input  and  accumulation  within  the 
channel  environment  and  investigate  the  stability  of  debris  jams  in  term  of  their  longevity  in  a 
given  reach. 

3)  To  investigate  how  effectively  debris  jams  inhibit  or  promote  bed  scour,  sediment  transport 
and  storage  to  determine  whether  they  are  net  stabilising  or  destabilising  elements  in  the  fluvial 
system. 

4) To  assess  the  impact  of  debris  at  run-of-river  hydraulic  structures  such  as  grade  controls, 
bridges,  bendway  weirs,  locks  and  dam  sluices. 

5)  To  develop  a  set  of  guidelines  for  in-channel  LWD  management  that  can  be  used  by 
engineers  and  river  managers  as  an  aid  to  assessment,  design  and  maintenance  of  stable 
channels. 

6)  To  develop  guidelines  for  LWD  management  technologies  at  run-of-river  structures. 

This  report  presents  state-of-the-art  review  of  literature  concerning  the  gee  norphological 
impact  of  LWD  and  LWD  management  strategies. 

Data  collection  was  initiated  through  a  three  week  survey  program  in  May  1991  with  the 
assistance  of  the  Colorado  State  University  DEC  monitoring  project  survey  crew.  A 
comprehensive  reconnaissance  survey  was  made  of  all  debris  jams  in  the  23  study  reaches  and 
each  site  was  also  surveyed  to  gather  the  long  profile  and  cross  section  data  needed  to  enable 
comparison  with  the  data  obtained  in  an  earlier  survey  of  May  1994.  Analysis  of  these  data  is 
include  in  this  report 

The  geomorphological  characteristics  of  jams  in  each  reach  have  been  analysed  and  plotted 
against  independent  catchment  variables  including  drainage  basin  area,  stream  power  and 
average  channel  top  width,  to  determine  whether  the  geomorphological  effects  of  LWD  have  a 
coherent  and  predictable  spatial  relationship.  Debris  jam  sediment  budgets  have  also  been 
calculated  and  related  to  spatial  parameters  to  determine  whether  the  net  impact  of  debris  jams 
is  through  sediment  retention  or  sediment  scour  and  mobilisation.  An  understanding  of  this 
impact  is  important  as  it  indicates  whether  LWD  is  a  net  stabilising  or  destabilising  agent  in 
unstable,  sand-bed  rivers. 

The  long-term  aim  of  this  research  is  to  improve  understanding  of  the  basin-wide  impacts  of 
I  ,WD  in  unstable  and  stable  channel  environments,  and  to  develop  coherent,  basin-wide  debris 
management  strategies  for  erosion  control  and  habitat  enhancement  and  to  propose  new 
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m'dintfinance/desi^n  procedures  for  DEC  and  run-of-river  structures,  based  on  sound 
geomorphic  and  engineering  analyses.  The  report  presents  final  conclusions  regarding  the 
management  of  LWD,  based  upon  the  results  of  this  three  year  research  program. 

Both  the  DMP  and  DB1P3  programs,  and  user  manuals  (written  in  MS  Word  ver.  7)  are 
enclosed,  inside  the  rear  cover  of  this  report 

The  CIS  user  interface  for  the  DMP  is  also  now  available  for  use  and  can  be  obtained 
upon  request  from  nick(%geography.nottingham.ac.uk 

The  program  runs  on  the  UNIX  version  of  ARC  INFO  and  requires  approximately  98 
mb  of  memory  to  install. 


2  LITERATURE  REVIEW 

2.1  INTRODUCTION 

Organic  or  woody  debris  is  an  important  channel  independent  variable  in  many  fluvial  systems 
(Hogan,  1987),  For  example,  Bevan  (1948;  quoted  in  Keller  and  Macdonald,  1995)  concluded 
that  in  the  Middle  Fork  Willamette  River,  Oregon,  woody  debris  was  responsible  for  more 
channel  changes  than  any  other  factor. 

In  a  literature  review  of  published  material  then  available,  Hickin  (1984)  suggested  that 
vegetation  may  influence  channel  processes  through  five  mechanisms: 

a)  Flow  resistance 

b)  Bank  strength 

c)  Bar  sedimentation 

d)  Formation  of  log  jams 

e)  Concave-bank  bench  deposits 

Hogan  also  identified  that  the  literature  concerning  this  subject  was  of  two  main  types;  that 
dealing  with  the  indirect  influence  relations  between  vegetation,  water,  sediment  yields  and 
river  morphology,  and  that  dealing  with  the  direct  impacts  of  channel  vegetation  on  channel 
morphology. 

Since  the  1980s  the  number  of  papers  dealing  with  vegetation  in  rivers  has  increased  markedly, 
however,  including  a  number  of  studies  concerning  Coarse  Woody  Debris  (CWD),  (Nakamura 
&  Swanson,  1993),  Large  Organic  Debris  (LOD)  (Hogan,  1987)  or  Large  Woody  Debris 
(LWD),  (Smith  &  Shields,  1992)  and  its  accumulation  as  jams  or  dams  in  river  channels. 

Studies  can  be  grouped  by  topic  into  those  dealing  primarily  with: 

a)  Input  processes,  distribution  and  residence  time  of  LWD 

b)  Geomorphic  significance  of  LWD 
c'<  Ecological  impact  of  LWD 

The  physical  processes  in  olved  in  each  topic  vary  depending  upon  the  size  of  the  stream 
relative  to  that  of  the  CWD  (Nakamura  et  al,  1993). 

Most  studies  have  been  carried  out  in  essentially  stable  channel  environments  in  the  US  and 
Canadian  Pacific  Northwest,  UK,  and  New  Zealand.  Instability,  in  the  form  of  iandsliding,  is 
cited  by  Poarce  &  Watson  (1981)  as  a  means  for  debris  to  enter  channels,  but,  more  generally, 
the  study  of  debris  impacts  in  inherently  unstable  channels  has  not  been  addressed 
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2.2  INPUT  PROCESSES,  FORMATION  AND  RESIDENCE  TIME  OF  LWD 

2.2.1  Input  Processes 

Large  Organic  Debris  enters  river  systems  by  two  main  processes,  either  from  outside  the 
channel  due  to  bank  erosion,  mass  wasting,  windthrow,  collapse  of  trees  due  to  ice  loading  or 
biological  factors  such  as  death  and  litter  fall  (Keller,  1979),  or  from  inside  the  channel, 
through  erosion  and  flotation  of  emergent  and  riparian  trees  (Hogan,  1987)  (Figure  2.1). 
Fetherston  et  al.  (1995)  suggest  that  debris  inputs  are  either  “chronic  or  episodic”.  Chronic 
inputs  are  frequent,  but  small  in  magnitude  and  occur  due  to  tree  mortality  and  bank  failure. 
Episodic  inputs  are  infrequent,  but  provide  a  large  amount  of  material.  Episodic  input 
processes  include  windthrow,  ice  storm,  fire  and  flood  events.  The  dominance  of  different 
input  processes  varies  widely  For  example,  45  percent  of  input  is  due  to  windthrow  in  the 
Lymington  Basin,  UE»  (Gregory  et  al,  1993),  while  massive  inputs  from  landsliding  of  debris  in 
a  mountain  catchment  are  reported  by  Pearce  &  Watson  (1983)  and  by  lansliding  as  a 
consequence  of  logging  operations  in  the  Queen  Charlotte  Island,  British  Columbia  by  Hogan 
et  al.  (1995).  Keller  et  al.  (1979)  suggest  that  in  low  gradient,  meandering  streams  inputs  are 
predominantly  the  result  of  bank  erosion  and  mass  bank  wasting,  windthrow  and  ice  loading, 
while  in  mountain  streams  the  main  process  is  debris  avalanche  Diehl  &  Bryan  (1994)  found 
the  dominant  input  process  to  be  bank  erosion  in  unstable  rivers  in  Tennessee  and  noted  that 
channel  instability  could  be  a  good  indicator  of  in-channel  debris  abundance  LWD  that  has 
been  input  by  bank  erosion  can  be  identified  and  distinguished  from  that  which  has  entered  by 
other  processes  because  the  trees  will  usually  have  an  asymmetrical  root  mass  due  to 
progressive  slipping  of  the  tree  from  the  bank  into  the  channel  (Diehl  ^  Bryan,  1994)  Smith  et 
al.  (1993)  found  debris  input  to  be  spatially  random.  However,  the  locations  of  zones  from 
which  LWD  is  supplied  will  vary  as  a  function  of  the  distribution  of  riparian  vegetation, 
streamside  topography,  channel  characteristics  and  the  prevailing  wind  strength  and  direction, 
(Fetherston  et  al.,  1995)  It  may,  therefore,  be  pos,>ible  to  determine  which  are  the  dominant 
input  factors  based  on  observations  of  these  factors  and,  thereby,  predict  the  distribution  of 
major  source  areas  within  the  catchment. 

2.2.2  Formation  of  Jams 

Once  in  a  channel,  debris  may  form  into  jams  oi  dams  Jams  usually  form  around  "key  coarse 
woody  debris"  (Nakamura,  1993),  which  are  usually  large,  whole  trees  that  have  entered  the 
channel  by  one  of  the  mechanisms  mentioned  above  and  which  may  be  anchored  to  the  bed  or 
banks  at  one  or  both  ends.  Smaller  debris  floating  down  the  channel  then  accumulates  against 
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the  key  elements,  v/hich  act  as  a  sieve  to  debris  and,  later  to  sediment.  If  there  is  no  fine  debris 
in  the  stream  a  mature  jam  may  never  form,  so  that  the  impact  of  key-debris  is  minimal.  The 
location  of  debris  jams  within  the  channel,  their  size  and  their  coherence  vary  as  functions  of 
position  in  the  catchment.  In  small  stieams  much  debris  will  accumulate  where  it  falls  because 
the  flow  is  not  competent  to  move  course  material  and  it  is  in  larger  streams  that  distinct  jams 
may  form  In  yet  larger  rivers  debris  may  never  accumulate  because  it  is  carried  away 
downstream.  Piegay  (1993)  observed  debris  distribution  in  a  sixth-order  river  in  France  and 
found  that  most  material  was  deposited  on  the  channel  margins,  forming  a  narrow  debris  line 
rather  than  in-flow  jams.  Wallance  &  Benke  (1984)  noted  a  similar  distribution  in  meandering 
rivers  in  the  southeast  USA  where  dense,  partial  jams  formed  at  a  an«le  to  the  main  flow.  As 
mean  channel  dimensions  and  flow  competence  increase  downstream  more  and  more  debris 
will  be  moved  from  its  position  of  input,  until  all  but  the  largest  trees  are  transported.  This 
process  relationship  may  result  in  a  trend  of  reducing  LWD  frequency  downstream  but,  at  the 
same  time,  produce  an  increase  in  the  volumetric  size  of  each  jam  (Swanson  et.  al.,  1982) 

Z.2.3  Residence  time  of  debris  jams 

The  residence  time,  or  persistence,  of  debris  jams  is  an  important  factor,  which  determines  the 
timespan  over  which  channel  morphology  at  a  jam  site  will  be  affected.  The  influence  exerted 
by  jams  on  channel  morphology  also  varies  with  time  as  the  debris  in  the  jam  structure 
deteriorates  (Hogan  et.  al ,  1995)  Assessing  residence  time  is  difficult  and  estimates  range 
between  12  months  for  a  36%  change  or  removal  (Gregory  &  Gurnell  (1985),  to  40-90  years 
(Hogan,  1987),  to  200  years  for  streams  in  British  Columbia  (Keller  &  Tally,  1979).  Residence 
times  may  vary  as  a  function  of  drainage  basin  area,  and  are  largely  dependent  upon  the  return 
period  of  a  flood  with  a  magnitude  capable  of  entraining  a  significant  proportion  of  the  trapped 
debris  or  moving  larger  key  components  of  the  jam.  Other  important  factors  affecting  jam 
persistence  are  average  tree  dimensions  and  wood  deterioration  rate.  Swanson  ci  al.  (1982) 
discovered  that  the  density  and  volume  of  in-channel  debris  are  greater  in  rivers  which  flow 
through  coniferous  forests,  than  they  are  in  rivers  that  flow  though  deciduous  forests.  This  is 
because  conifers  are,  on  average,  taller  and  have  slower  decay  rates  than  deciduous  trees. 
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Figure  2.1  Dyaaitiics  of  woody  debris  (adapted  from  Keller  &  Sw'anson,  1979) 
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Figure  2.2  Energy  transformations  in  streams  with  log  steps  (adapted  from  Marston, 
1982) 
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2.3  GKOMORPHIC  SIGNIFICANCE  OF  LWD 

2.3.1  Eflects  of  channel  scale 

It  is  important  to  recognise  that  processes  are  scale-dependent  and  that  the  influence  of  LWD 
on  channel  and  valley  morphology  may  change  systematically  downstream  through  the 
network  (Abbe  &  Montgomery,  1993)  Zimmerman  et  al.  (1967)  found  that  debris 
accumulations  in  a  very  small  stream  completely  obscured  the  usual  hydraulic  geometry 
relations,  while  Robinson  &  Beschta  (1990),  and  Keller  &  Tally  (1979)  suggest  that  bris 
loadings  increase  with  stream  size.  Gregory  et  al.  (1985),  have  characterised  jams  into  three 
types  : 

1 )  Active  (form  a  complete  barrier  to  water  and  sediment  movement,  and  create  a 

distinct  step  or  fall  in  the  channel  profile) 

2)  Complete  (a  complete  barrier  to  water/sediment  movement,  but  no  step  formed) 

3)  Partial  (only  a  partial  barrier  to  flow) 

They  suggest  that  these  types  become  more  prevalent  sequentially  as  channel  size  increases. 
Once  trees  fall  into  a  stream,  their  influence  on  channel  form  and  process  may  be  quite  different 
to  that  when  they  v.'ere  on  the  banks,  changing  from  a  stabilising  to  a  destabilising  influence 
due  to  their  effect  in  causing  causing  local  bed  scour  and  basal  erosion  of  the  banks.  Thus,  jams 
represent  a  type  of  auto-diversion,  that  is  a  change  in  channel  morphology  triggered  by  the 
fluvial  process  itself  (Keller  &  Swanson,  1979). 

The  type  and  degree  of  impact  on  channel  morphology  depends  primarily  on  the  channel 
width/tree  height  ratio  and  on  debris  orientation  relative  to  the  flow  Mean  discharge  and  the 
dominant  discharge  recurrence  interval  are  also  important  because  the  higher  the  flow  is 
relative  to  jam  size,  the  smaller  will  be  the  jam’s  impact  through  flow  diversion  and  channel 
roughening  The  principal  effects  of  debris  upon  channel  morphology  are  described  below. 

2.3.2  Impact  of  debris  jams  upon  channel  morphology 

LWD  influences  the  geomorphology  of  rivers  at  three  levels  (Gray,  1974):  the  overall  channel 
form;  detailed  features  of  the  channel  topography,  and,  channel  roughness. 

Heede  (1985),  Smith  (1993),  Andrus  et  al.  (1988)  and  Mosley  (1981)  have  all  observed  that 
the  spatial  distribution  and  number  of  pools,  riffles  and  gravel  bars  is  positively  related  to  the 
distribution  and  volume  of  LWD  in  the  channel.  This  relationship  has  been  explained  through 
laboraton,  experiments  by  Smith  &  Beschta  (1994),  who  found  that  the  pool-ritfle  sequence  in 
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Figure  2.3  Schematic  diagrams  of  pool  typos  (Modified  from  Robinson  and  Beschta, 
!990) 
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gravel-bed  rivers  is  maintained  by  a  combination  of  mean  boundary  shear  stress  and 
intermittent  lift  and  drag  forces  due  to  velocity  fluctuations  around  debris.  Randonfi  debris 

'4  input  will  also  distort  the  pool-riffle  sequence,  making  it  less  systematic,  so  that  the  long- 

profile  h;!s  very  little  spatial  memory'  or  periodicity  (Robinson  &  Beschta,  1990).  Robinson  and 
Beschta  (1990)  devised  a  pool  classification  system  containing  six  pool  categories  (lateral 
scour,  fluvial,  plunge,  underflow,  deflector  and  dam)  based  on  flow  and  debris  interaction  (see 
figure  2.3).  Other  studies  have  shown  that  a  considerable  proportion  of  the  vertical  fall  of 
channels  can  occur  at  the  sites  of  debris  jams,  accounting  for  a  4%  of  the  vertical  drop  along  a 

412m  reach  of  channel  in  Vermont  (Thompson,  1995)  and  60%  of  the  total  drop  in  Little  Lost 
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Man  Creek  in  Northern  California  (Keller  &  Tally,  1979).  Debris  jams,  therefore,  may  act  as 
local  base  levels  and  sediment  storage  zones  which  provide  a  buffer  in  the  sediment  routing 
system  (Heede,  1985,Bilby,  1981). 

Thompson  (1995)  found  that  LWD  causes  an  important  negative  feedback  mechanism  where, 
in  the  case  of  channel  degradation,  there  is  an  increase  in  debris  input  due  to  mass  bank  failure, 
which  ii.  turn  causes  greater  sediment  storage,  Increased  sediment  storage  causes  channel  bed 
elevation  lo  rise,  tending  to  reduce  bank  heights  and  the  rate  of  bank  failure  and  debris  input  is, 
therefore  reduced.  On  this  basis,  Klein  et  al.  (1987)  argue  that  jam  removal  can  reduce  the  base 
level  for  the  channel  upstream  and  may  trigger  bank  erosion.  However,  in  an  experimental 
study  by  Smith  et  al.  (1993a  and  b)  it  was  found  that,  while  the  removal  of  debris  from  a  small 
gravel  bed  stream  initially  caused  a  four-fold  increase  in  bed  load  transport  at  bankfull  flow,  the 
associated  loss  of  scour  turbulence  and  greater  flow  resistance  imparted  by  alternate  bars 
actually  resulted  in  a  reduction  in  stream  power  which  was  compensated  for  by  sediment 
deposition  and  net  channel  aggredation. 

Potential  energy  is  dissipated  at  jams  and  jam  energy  loss  may  be  as  much  as  6%  of  total 
potential  energy  (MacDonald  et  al.,  1982).  Shields  &  Smith  (1992)  found  that  the  Darcy- 
Weisbach  friction  factor  was  400  %  higher  at  base  flow  in  an  uncleared  river  reach  compared 
to  a  cleared  condition,  but  that  this  value  declined  to  35%  at  high  flows.  The  velocity 
distribution  is  also  fru  more  heterogeneous  in  debris-filled  reaches,  especially  at  low  flow. 
Changes  of  stream  power  distribution  due  to  flow  resistance  effects  in  turn  give  jams  the  ability 
to  influence  the  location  of  erosional  and  depusitional  processes.  Also,  the  backwater  effect 
cieated  by  jam  back-pools  may  induce  local  silting  (Keller  et  al.  1976)  Thus,  in  small,  stable 
channels  log  steps  generally  increase  bank  stability  and  reduce  sediment  transport  rates  by 
creating  falls,  runs  and  hydraulic  jumps. 

Figure  2.2  shows  how  potential  energy  is  lost  through  a  log-step  sequence,  as  outlined  by 
Marston  (1982).  The  localised  dissipation  of  energy  can,  however,  result  in  associated  local 
scour  and  bank  erosion  which  causes  channel  widening.  Bank  failure  may  also  occur  through 
flow  diversion  around  a  debris  obstruction  (Murgatroyd  &  Ternan,  1983).  Davis  &  Gregory 
(1994)  have  also  suggested  a  mechanism  whereby  bank  failure  is  induced  through  the  erosion 
of  a  porous,  gravel,  bank  subsurface  due  to  the  greater  hydrostatic  pressure  caused  by  debris- 
dammed  flow.  C’onversely,  Keller  &  Tally  (1979)  have  observed  that  flow  convergence  under 
logs  may  cause  channel  narro'ving,  with  sediment  storage  upstream  and  a  scour-pool 
downstream  of  the  log  step. 
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As  drainage  area  increases  and  the  channel  width/tree  size  ratio  exceeds  unity,  flow  is  diverted 
laterally,  to  induce  bank  erosion  through  local  basal  scour.  Hogan  (1987)  found  that,  in 
undisturbed  channels  in  British  Columbia,  organic  debris  orientated  diagonally  across  the 
channel  resulted  in  high  width  and  depth  variability  However,  in  catchments  v'here  there  had 
been  logging  operations  the  majority  of  in-channel  discarded  timber  was  orientated  parallel  to 
the  flow  and  it  subsequently  became  incorporated  into  the  stream  banks,  protecting  them  from 
erosion.  Nakamura  &  Swanson  (1993)  and  Keller  &  Swanson  (1979)  have  suggested  that 
there  is  a  wide  range  of  types  of  interaction  between  debris  jam  and  channel  processes, 
progressing  from  local  base  level  control  and  possible  local  widening  in  low-order  streams,  to 
lateral  channel  shifts  and  even  meander  cut-off  in  middle-order  channels,  where  debris  is 
moved  into  larger  more  coherent  jams  which  may  either  increase  or  decrease  the  channel 
stability  depending  upon  the  erodibility  of  bed  and  banks.  In  larger  channels  still,  bars  may 
form  and  flow  bifurcate  around  debris  obstructions  This  last  process  has  been  documented  by 
Nanson  (1981)  in  British  Columbia,  who  found  that  organic  debris  deposited  at  low  flow 
provided  the  nuclei  for  development  of  scroll  bars,  through  the  local  reduction  of  stream 
power.  Hickin  (1984)  also  observed  crib-like  bar-head  features,  but  was  undecided  regarding 
whether  the  debris  caused  bar  formation,  or  whether  the  bars  pre-dated  and  trapped  the  debris. 
In  either  case,  organic  debris  would,  at  the  very  least,  enhance  sediment  deposition  and  bar 
formation. 

2.4  ECOLOGICAL  IMPACT  OF  LWD 

LWD  dams  are  very  important  in  small  stream  ecosystems  because  they  provide  a  source  of 
organic  matter  and  retain  floating  leaves  and  twigs  in  the  dam  structure  and  backwater  pools. 
This  coarse  particulate  organic  matter  (CPOM)  is  broken  down  by  shredder  invertebrates  in 
the  low  energy  pool  environment,  creating  fine  particulate  organic  matter  (FROM)  and 
dissolved  organic  matter  (DOM),  which  are  the  required  energy  sources  of  a  succession  of 
invertebrate  species  who  are,  in  turn,  the  energy  source  of  high  fauna  species  Bilby  &  Likens 
(1980)  found  that  the  percentage  of  the  standing  stock  of  organic  matter  retained  by  jams 
changed  from  75%  in  first-order,  to  58%  in  second-order,  and  20%  in  third-ordc!  streams 
because  the  prevalence  of  dam  type  jams  declined  downstream  The  volume  of  CPOM 
therefore  declines  downstream,  while  the  volume  of  FROM  and  DOM  increases  This  gives  rise 
to  a  spatially  varied  invertebrate  community,  changing  from  shredders  in  small  channels  to 
gathers  of  FROM  downstr-eam.  .Smock  et  al.  (1989)  and  Wallance  &  Benke  (1984)  found 
similar  correlations  between  debris  volume  and  invertebrate  abund,ance  in  sand-bed  streams. 
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where  debris  provides  the  only  stable  substrate  for  organic  matter  retention  and  invertebrate 
habitat.  Higher  species,  such  as  fish,  use  debris  and  associated  pools  for  shade,  protection  from 
predators,  feeding  and  spawning  grounds.  The  pools  and  tails  created  by  log  steps  also  help  to 
oxygenate  the  flow,  and  provide  a  variety  of  different  energy  environments  which  are  can  be 
colonised  by  niche  species. 

In  addition  to  providing  essential  fauna  habitat,  LWD  is  also  a  vital  factor  in  the  development 
of  the  riparian  forest  mosaic  (Fetherston  et  al ,  1995)  Debris  deposition  in  the  channel  and  on 
the  floodplain  creates  sites  of  low  boundary  shear-stress  where  vegetation  colonisation  can 
take  place.  This  leads  to  the  development  of  vegetation-stabilised  islands  and  bars  (affecting 
the  geomorphological  development  of  the  channel)  which  may  subsequently  coalesce  and/or 
become  attached  to  the  bankline  to  form  new  areas  of  forested  floodplain  that  provide  shade, 
bank  stability  and  supply  and  stot  age  of  organic  matter,  sediment,  water  and  new  LWD,  Figure 
2,4  shows  a  modified  version  of  the  LWD-driven  model  of  Fetherston  et  al ,  (1995)  for 
riparian  forest  development,  based  upon  research  findings  from  the  Pacific  Northwest 


Figure  2.4  Conceptu.vl  model  of  montane  riparian  forest  development.  (Modified  from 
Fetticrston  et  a!.,  1995) 
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2.5  MANAGEMENT  STRATEGIES 

Until  basic  research  concerning  in-channel  LWD  began  to  suggest  otherwise,  it  was  commonly 
believed  that  LWD  was  detrimental  to  the  fluvial  system,  hydraulically,  ecologically  and 
geomorphically.  On  this  basis,  reasons  for  debris  removal  included: 

a)  To  improve  navigation, 

b)  1 0  increase  channel  conveyance  by  reducing  roughness; 

c)  To  eliminate  bank  erosion; 

d)  To  facilitate  the  migration  of  fish,  especially  salmon  (MacDonald,  1982). 

It  is  now  recognised  that  there  are  advantages  to  be  gained  by  maintaining  or  even  increasing 
in-channel  debris  accumulations  (Gregory  &  Davis,  1992;  Keller  &  McDonald,  1995). 
Management  strategies  that  are  currently  advocated  vary  widely,  however  This  perhaps 
reflects  our,  as  yet,  incomplete  understanding  of  LWD  dynamics  in  different  channel 
environments,  and  because  goals  vary  between  different  management  strategies.  In  this  respect 
effective  debris  management  depends  on  the  underlying  aims  of  the  proposed  management 
action. 

Successful  management  also  depends  upon  a  comprehensive  understanding  of  the  following 
hydrogeomorphological  factors  (Gregory  &  Davis,  1992): 

a)  The  relationship  between  river  channel  processes  and  river  channel  morphology, 

b)  Awareness  of  the  timescales  over  which  river  channels  may  adjust; 

c)  Consideration  of  channel  management  in  the  wider  context  of  river  basin 
management. 

More  specifically,  debris  management  must  consider: 

a)  Channel  stream  power  characteristics, 

b)  Sediment  movement  and  storage  relationships  (high/low;  tine/coarse  sediment, 

suspended/bedload), 

c)  Channel  stability; 

d)  Size  and  character  of  river  channel  in  relation  to  debris  size, 

e)  Spacing  and  frequency  of  jams; 

f)  Size  and  character  of  jams,  and  orientations  of  component  material, 

7)  Age  and  stability  of  component  materials. 

In  an  evaluation  of  soft  engineering  for  in-stream  structures,  including  some  using  woody 
debris,  to  mitigate  the  effects  of  highway  construction  in  British  C'olumbia,  Miles  ( 1995)  found 
that  nearly  50%  of  the  structures  had  been  severely  damaged  after  8  to  14  years  Miles 


attributed  this  problem  to  insufficient  understanding  and  consideration  of  the  stability  of  the 
structures  in  a  high  energy  river  environment.  He  concluded  that  soft  restoration  techniques 
may  not  be  appropriate  in  highly  energetic  mountain  rivers  and  that,  if  restoration  is  to  be 
performed,  funding  must  be  made  available  for  long  term  monitoring  and  maintenance. 

There  appears,  in  general,  to  be  a  consensus  of  opinion  amongst  researchers  interested  in  LWD 
regarding  appropriate  management  approaches  for  channel  restoration.  Bren  (19y3)  and 
Nunnally  (1978)  argued  that  the  riparian  zone  should  be  left  undisturbed,  in  a  natural  state 
(although  defining  natural  is  difficult  in  most  channels)  and  that,  because  debris  is  so  important 
for  the  river  ecosystem,  debris  jams  should  be  left  in  place.  Keller  and  McDonald  (1995) 
studied  catchments  which  had  been  disturbed  by  logging  operations.  They  recommended  that  a 
riparian  buffer  strip  should  be  left  to  maintain  the  natural  LWD  supply  and  warned  that 
landsliding  events  caused  by  badly  controlled  logging  operations,  can  cause  excessive  l.WD 
input  which  is  detrimental  to  stream  habitat  and  flow  and  sediment  conveyance. 

There  may  be  a  case  in  streams  lacking  a  wooded  riparian  strip  for  the  introduction  of  debris 
jams  (Keller  &  McDonald,  1995).  If  a  debris  recharge  policy  is  to  be  implemented,  however,  it 
is  important  that  debris  jam  volume  and  orientation  emulates  the  values  which  would  be  found 
under  natural  conditions  (Robinson  &  Beschta,  1990).  Wallace  &  Benkc  (1984)  concluded 
that,  in  most  instances,  the  best  management  is  probably  no  management,  except  where 
adjacent  floodplains  have  to  be  protected  from  flooding 

Comprehensive  studies  of  coarse  woody  debris  in  relation  to  river  channel  management  have 
been  carried  out  by  Gregory  and  Davis  (1992)  and  Gurnell  and  Gregory'  (1995a  and  b)  The 
collation  of  analyses  from  twenty  two  research  papers  with  primary  field  studies  carried  out  in 
the  New  Forest,  UK  (Gregory  &  Davis,  1992)  demonstrated  the  significance  of  LWD  to 
channel  morphology,  processes  and  ecology  (Figure  2.5)  and  produced  a  preliminary  table  of 
debris  management  criteria  based  upon  their  findings  (Figure  2.6).  They  conclude  that,  "...  a 
conservative  approach  to  debris  removal  should  be  adopted  for  most  areas,  but  that  different 
strategies  are  needed  according  to  the  characteristics  of  particular  localities."  (Gregory  and 
Davis,  1992,  pg.  133). 

It  should  be  noted,  however,  that  this  study,  in  common  with  most  others  cited,  was  carried 
out  in  an  essentially  stable,  equilibrium  channel  environment  where  changes  to  channel 
morphology  are  negligible  and  significant  impacts  relate  mostly  to  ecological  habitat  diversity. 
Also,  little  attention  is  paid  to  "different  strategies"  that  may  be  required  in  contrasting  channel 
environments  and  there  is  no  discussion  of  conflicts  between  practices  advocated  by  various 
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organisations  in  the  USA,  For  example,  Gregory  &  Davis  (1992)  suggest  that,  based  on  their 
literature  survey,  no  debris  should  be  removed  from  channels  exhibiting  low  stability  (Figure 
2,7).  However,  this  contradicts  the  practice  described  by  Brookes  (1985,  pg,  64),  "In  North 
America  the  concept  of  channel  restoration  was  developed  in  North  Carolina  under  the  funding 
of  the  Water  Resources  Research  Institute  of  the  State  University  ...Restoration  is  achieved  by 
removing  debris  jams  and  providing  uniform  channel  cross-sections  and  gradients  whilst 
preserving  meanders,  leaving  as  many  trees  as  possible  along  the  stream  banks,  and  stabilising 
banks  with  vegetation  and  rip-rap  where  necessary  .  . 

Similar  approaches  have  been  documented  and  carried  out  by  numerous  researchers  and 
organisations  in  the  USA  including;  McConnel  et  al,  (1980),  based  upon  work  on  the  Wolf 
River,  Tennessee;  the  American  Fisheries  Society  (1983),  in  a  publication  entitled  “Stream 
Obstruction  Removal  Guidelines”,  (see  Figure  2.8);  Shields  and  Nunnally  (1984);  and,  Palmiter 
(Institute  of  Environmental  Sciences,  1982). 

The  recommendations  of  Palmiter  (1982)  include  the  following: 

a)  Removal  of  log-jam  material  by  cutting  it  to  a  manageable  size; 

b)  Protection  of  eroding  banks  using  brush  piles  and  log-jam  material,  with  rope  and 

wire; 

c)  Removal  of  sand  and  gravel  using  brush-pile  d  jtlcctors; 

d)  Revegetation  to  stabilise  banks  and  shad  >out  aquatic  plants; 

e)  Removal  of  potential  obstructions  such  as  trees  and  branches, 

Willeke  &  Baldwin  (1984)  assessed  the  Palmiter  techniques  and  found  them  suitable  for  areas 
experiencing  chronic,  low  intensity  flooding  and  bank  erosion,  but  not  advisable  for  rivers  with 
extreme  flood  problems  They  are  also  found  to  be  largely  ineffective  for  erosion  control  where 
the  mechanism  of  bank  failure  is  that  of  mass  wasting  rather  than  tri  five  force  erosion 
(Hasselwander,  1989), 

It  is  evident  from  the  preceding  discussion  of  LWD  management  strategies  that 
recommendations  vary  considerably  ranging  from  limited  or  no  interference  to  total  clearance 
of  debris  from  the  channel.  These  apparently  contradictory  recommendations  must  be  viewed 
in  the  light  of  the  overall  management  program  that  they  were  designed  for,  as  requirements 
for  habitat  enhancement  differ  from  those  for  flood  defence. 

Finally,  and  of  great  importance,,  is  the  fact  that  the  recommendations  of  type  made  by 
Palmiter  and  others,  address  debris  management  predominantly  in,  low  gradient,  sand-bed,  and 
perhaps  iin.stable,  flood  prone  rivers  (South  tdSt  USA),  while  those  prescribed  by  Gregory  and 


Davis  (1992)  and  others  aie  based  upon  findings  from  upland,  even  montain,  gravel-bed  rivers 
and  streams  (Pacific  Northwest  USA).  Process  relationships  between  the  debris  and  the 
channel  are  likely  to  differ  between  these  two  fluvial  environments,  although,  as  yet,  these 
differences  have  not  been  recognised  or  investigated.  Indeed,  white  there  is  a  wealth  of 
research  concerning  the  geomorphological  impacts  of  LWD  in  upland  giavel-bed  rivers,  there 
has  been  little  comparable  research  in  lowland,  sand-bed,  and/or  unstable  river  environments. 
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Figure  2.5  The  significance  of  coarse  woody  debris  dams  for  channel  morphology, 
channel  processes  and  ecology  (modified  from  Gregory  and  Davis,  1992) 
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Figure  2.6  Determinants  for  a  management  strategy  for  rivers  in  woodland  areas 
(modified  from  Gregory  and  Davis,  1992) 
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Figure  2.7  Derinition  of  Stream  Obstruction  Conditions  (Modified  from  American 
Fisheries  Society,  1983) _ 


Condition  One 

llicst:  stream  segments  have  acceptable  How  and  no  work 
would  l)e  required.  They  may  contain  various  amoimts  of 
instreum  debris  and  fine  sediment,  such  as  silt,  sand, 
gravel,  rubble,  boulders,  logs  and  brush.  In  certain 
situations  How  may  be  impeded,  but  due  to  stream  and  land 
classification  or  adjacent  land  mse,  this  is  not  a  problem 


Condition  Two 

These  stream  segments  currcnitly  have  no  major  How 
impediments,  but  existing  conditions  arc  such  tliut 
obstructions  ore  likely  to  form  in  the  near  future,  causing 
unacce^itable  problems.  This  condition  is  generally 
cbaructerised  by  small  accumulations  of  logs  and/or  otlier 
debris  which  occasionally  span  the  entire  stream  width. 
Accumulations  are  isolated,  not  massive  and  do  not 
presently  cause  upstream  ponding  damage. 


Manaiicinent  Criteria 

No  work  to  be  conducted. 


Condition  Three 

fhe^c  stream  segments  have  unacceptable  How  problems. 
Obstructions  are  generally  characterised  by  large 
accumulations  of  lodged  trees,  root  wads,  anJ/or  other 
debns  Chat  frequently  span  the  entire  stream  widUt. 
Although  impeded,  some  How  moves  through  the 
obstruction,  Large  amounts  of  sediment  have  not  covered  or 
hxiged  in  the  obstr  uction 


Condition  l<'our 

llicsc  slrcain  segments  arc  characterised  by  major 
blockages  cuising  luiacceptable  How  problems. 
Obstructions  consist  of  compacted  debris  uml/or  ilebris  (hut 
severely  restricts  How 


Condition  Five 

Idie-sc  stream  segments  po.^slJ.^  imiqiie,  .sensitive,  or 
especially  valuable  biotic  resources  and  should  be  dealt 
with  on  a  case-by-case  basis,  lixamples  include,  but  are  not 
limited  to:  Arcus  harbouring  rare  or  endungcTcd  .species, 
.shell!i.sh  beds,  fish  spawning  and  rearing  ureas,  and 
rookeries. 


Manti}>cnii'nt  Criteria 

liquipment  that  will  cause  the  lea.st  diunuge  to  tile 
environment  shall  be  selected  for  perform  iig  tlte  work. 
First  consideration  will  be  given  to  the  use  of  hand  oi>ei'ateJ 
equipment  such  us  axes,  chain  saws,  .vtid  winches  to  ren.cve 
accumulations.  Bouts  witli  motors  may  he  used  where 
needed.  When  the  use  of  hand  operated  looks  is  not  feasible, 
heavier  equipment  may  be  used.  e  g.  snuill  tiactors, 
backhoes,  bulldo/ers,  log  skidders  and  low  PS'  equipment. 
I-quipment  shall  be  operated  in  u  maimer  Uiat  results  in 
least  duinugc  to  vegetation  and  soils  of  the  project  men.  hi 
some  cases  cxplosive.s  may  be  u.sed  lesuiling  in  less 
damage.  Debris  designated  for  removal  from  the  stream  or 
11ui.,lw'ay  should  be  removed  or  secured  in  such  a  inamicr  as 
to  restrict  its  rcHintry  into  the  chuimcl  tlencrally,  it  should 

be  positioned  so  as  to  reduce  Hood  How  iuipcdimeiit _ 

MHnuKeinent  Criteria 

liquipment  limitations  will  be  the  same  as  fur  euiiditiun  two 
segments.  Work  shall  be  accomplished  within  the  clumnel 
or  from  one  side  of  the  chaiuiel  whete  possible.  Selective 
tree  clearing  shall  be  limited  to  the  minimmn  necessary  for 
equipment  access  and  elficient  oi>eration  of  equipment  on 
the  worked  side  ol'the  chaiuicl.  Disposal  of  equipment  may 
be  aeeoinplished  by  removing  it  from  the  Hixxlplain  or  by 
huiiiing,  burying  or  piling,  as  appropriate,  with  the 
minimum  amount  of  disturbance  to  vegetation.  Piled  debns 
shall  be  gapped  at  Irequcnl  intervals  and  at  all  tributaries 
and  dislributarie.s. 


Miuiiigeinent  Criteria 

Blockage  removal  muv  ennploy  any  equipment  nete.ssary  to 
accomplish  the  work  in  tlie  least  damaging  miuiuer.  Work 
should  be  accomplished  I'rom  one  side  of  the  elianiiel  where 
praetieal.  Material  shall  be  disposed  in  aeeordaiiee  with 
guidelines  presented  above  for  condition  three  .segments. 
Spoil  piles  should  be  constructed  as  higli  us  sediment 
pro|x:rties  allow,  liie  placement  of  spoil  around  the  base  of 
mature  trees  should  be  avoided. 


IVIiinjivement  Criteria 

S|ieciul  provision  for  protecting  unique,  sensitive,  or 
prrxluctive  biotic  re,sourccs  .shall  be  developed  by 
appropriate  professionals  on  case  by  case  basis. 
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3  REGIONAL  CHARACTERISTICS 

3.1  Introduction 

Creeks  in  northern  Mississippi  have  received  considerable  attention  from  geoinorphologists 
and  engineers  recently  in  an  effort  to  try  and  stabilise  channels  which  have  been  degrading 
through  knickpoint  migration  due  to  lowering  of  the  base  level  (Schumm,  1984)  Degradation 
problems  have  been  combated  by  the  construction  of  grade-control  structures  which  reduce 
stream  power  by  means  of  a  hydraulic  jump.  Degradation  has  been  caused  in  the  past  forty 
years  by  several  factors  including  trunk  stream  dredging,  changes  in  farming  practices  and 
channel  straightening  which  was  originally  carried  out  to  reduce  flooding  and  improve  drainage 
of  valley  bottom  lands.  The  degradation  problem  is  manifest  mainly  in  the  bluftline  hills,  away 
from  the  Mississippi  delta  and  floodplain,  where  creeks  are  often  ephemeral  and  "flashy"  in 
nature,  but  flow  through  highly  erodible  alluvium  and  loess. 

The  study  sites  chosen  all  lie  on  this  type  of  material,  in  the  Bluff-Line  Hills.  Plate  3  2  shows  a 
geological  map  of  Mississippi  and  the  location  of  the  DEC  study  area. 

3.2  GEOLOGY 

3.2.1  Pleistocene-Holocene  Stratigraphy 

Four  major  stratigraphic  units  make  up  the  Pleistocene-Holocene  valley  fill  in  northern 
Mississippi  From  oldest  to  youngest,  these  are; 

1)  Meander  Beltl(MBl): 

A  fine-grained  deposit  which  contains  five  types  of  sediment;  point  bar,  channel,  natural  levee, 
abandoned  channel,  and  backswamp  If  natural  levee  and  abandoned  channel-till  sediments 
crop  out  in  the  toe  slope  they  provide  a  measure  of  stability,  whereas  point  bar,  channel,  and 
backswamp  sediments  are  highly  erodible.  The  level  of  channel  degradation  has  now 
progressed  to  the  point  where  toe-slope  materials  are  comprised  of  these  units. 

2)  Lacustrine  (Old  Palaeosol)  (L) ; 

Deposited  in  a  low  energy  fluvial  environment,  these  sediments  are  homogeneous  and  comprise 
loess  and  alluvially  re-worked  loess.  Characteristic  of  this  unit  is  a  well-developed  polygonal 
structure,  with  seams  up  to  2  cm  thick,  filled  with  unconsolidated  material.  These  seams 
present  distinct  planes  of  weakness  in  the  material.  The  unit  also  contains  iron  concretions  in 
the  lower  level  which  offer  some  resistance  to  knick-point  migration  (Little  et  al.,  1982)  In 
this  study  Old  Palaeosol  outcrops  were  recognised  in  both  Lick  Creek  and  Nolehoe  Creek. 
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3)  Meanderbelt  2  (Young  Palaeosol)  (MB2) : 

These  materials  are  fine-grained,  weathered,  less  cohesive,  and,  therefore,  more  erodible  than 
Meanderbelt  1.  These  sediments  form  the  upper  portion  of  many  banks  and  represent  a 
significant  source  of  suspended  and  bedload  sediments.  This  layer  is  highly  permeable 
compared  to  the  Old  Palaeosol,  and  piping  often  occurs  at  the  L/MB2  interface  Piping  is  an 
important  contributor  to  bank  failure  (Grissinger,  1982). 

4)  Post  Settlement  Alluvium  (PSA) : 

PSA  was  deposited  during  the  last  150  years  due  to  the  rapid  vertical  accretion  of  floodplain 
sediments  following  European  settlement  of  the  region.  Up  to  16  ft  of  deposition  has  occurred 
due  to  farming  malpractice  that  commented  in  the  1830s  and  continued  until  the  1930s  The 
bulk  of  sediment  has  been  deposited  in  the  upper  and  middle  reaches  of  valleys.  These 
sediments  are  unweathered  and  cohesionless  and  provide  a  principle  source  of  suspended 
sediments.  This  Pleistocene-Holocene  stratigraphy  is  shown  in  as  a  column  in  Figure  3.1a  and 
in  cross-section  in  Figure  3.1b. 

3.3  GEOMORPHOLOGY  OF  CHANNELS  IN  NORTHERN  MISSISSIPPI 

3.3.1  Bank  Characteristics 

In  all  four  creeks  studied,  the  Pleistocene-Holocene  MBl,  L  and  MB2  units  are  located  within 
the  channel  banks  to  a  depth  of  4-8  feet.  This  loess  material  is  fairly  uniform  in  color,  and 
composed  of  a  brown-gray-sandy  silt  to  silty  clay  (Schumm  ct  al.,l‘)84).  The  chief 
characteristics  of  this  loess  are: 

a)  uniformity  of  texture, 

b)  Irregularity  of  shape  and  extreme  fineness  of  particles; 

c)  Generally  massive  structure; 

d)  Lack  of  coherence, 

e)  Capacity  to  stand  as  vertical-faced  walls; 

f)  Capacity  to  absorb  water. 

3.3.2  Bed  and  Bank  Failure  Mechanisms  in  Degrading  Channels 

Knick  zones  can  be  observed  in  nearly  all  drainage  basins  in  northern  Mississippi.  I’hese 
features  reflect  an  abrupt  change  in  channel  gradient,  the  typical  knick-point  unit  being 
characterised  by  a  2-4  ft  vertical  fall  within  a  cohesive  unit,  which  reflects  an  overall  incision  of 


the  channel  in  response  to  regional  lowering  of  the  base  level  (Schumm  et  al ,  1984),  Perched 
tributaries,  drainage  ditches  which  enter  the  main  stream  at  a  higher  elevation  than  the  bed  of 
the  main  channel  and  mature  tree  growth  significantly  higher  than  the  channel  bed  are  also 
characteristic  of  overall  incision.  The  bed  of  the  middle  to  upper  portion  of  Byhalia  Creek,  for 
example,  is  dominated  by  knick  zones,  perched  tributaries  and  undercut  structures. 

The  principal  units  being  degraoed  in  northern  Mississippi  are  Meanderbelt  1  and  Lacustrine 
However,  natural  levee  and  abandoned  channel  sediments  of  MBl  are  more  resistant  to 
erosion  than  channel,  point  bar  and  backswamp  sediments  and  offer  some  resistance  to 
knickpoint  migration.  Upstream  knickpoint  migration  is  a  function  of  flood  frequency  and  bed 
material  erodibility.  Where  a  knickpoint  has  been  at  the  same  location  over  a  period  of  time,  a 
characteristic  widening  of  the  channel,  downstream,  takes  place  This  is  known  locally  as  a 
"blow-out"  (Begin  et  al ,  1981). 

A  sequence  of  channel  reach  types  associated  with  knickpoint  migration  in  northern  Mississippi 
has  been  conceptualised  by  Schumm  et  al.  (1984),  in  what  has  been  termed  the  Channel 
Evolution  Mode  (CEM)  This  idealised  sequence  is  shown  schematically  in  figure  3.1c.  Along 
a  given  channel  morphological  classification  Types  1  through  5  occur  in  sequence  downstream 
and,  at  a  given  location,  will  alsooccur,  in  series,  over  time.  The  data  used  to  determine  the 
channel  evolution  model  were  obtained  from  Oaklimiter  Creek,  and  are  presented  in  table  3.1 
The  morphological  categories,  as  described  by  Watson  et  al.  (1993)  arc  outlined  below 

•  lype  i  These  reaches  are  in  equilibrium  and  are  characterised  by  a  sediment  transport 
capacity  that  exceeds  sediment  supply,  a  bank  height  that  is  less  than  the  critical  height,  a 
saucer  shaped  cross-section,  possibly  small  precursor  knickpoints  and  a  channel  bed  with  little 
or  no  sediment  deposited  Width  to  depth  ratios  are  highly  variable. 

•  Type  2  :  These  reaches  are  located  immediately  downstream  of  primary  knickpoints  and  are 
characterised  by  sediment  capacity  that  exceeds  sediment  supply,  a  bank  height  that  is  equi  .1  to 
the  critical  bank  height,  rectangular  cross-section,  a  channel  bed  with  little  or  no  sedi’ rent 
deposited,  a  lower  bed  slope  than  Type  1  reaches  because  degradation  reduces  bed  slope  and  a 
lower  width-depth  ration  value  than  Type  1  reaches  because  the  depth  has  increased  but  the 
banks  are  not  failing. 

•  Type  3  :  These  are  locat  d  downstream  of  Type  2  reaches  and  are  characterised  by  a 
sediment  transport  capacity  that  is  highly  variable  with  respect  to  sediment  supply,  a  bank 
height  that  is  greater  than  the  critical  height  and  erosion  due  primarily  to  slab  type  failure,  an 


accumulation  of  less  than  two  feet  of  sediment  in  the  bed  channel  which  locally  can  be  much 
higher  if  there  have  been  recent  bank  failures  that  have  not  been  removed  by  subsequent  flows, 
a  channel  depth  that  is  somewhat  less  than  that  in  the  Type  2  reaches  due  to  limited 
aggradation  and  channel  widening  due  to  bank  erosion. 

•  Type  4  :  These  are  located  downstream  of  Type  3  reaches  and  are  characterised  by  a 
sediment  supply  that  exceeds  transport  capacity  cause  aggradation  of  the  channel  bed,  the 
formation  of  longitudinal  berms  along  the  channel  margins  which  eventually  define  the  edge  of 
the  effective  discharge  channel  and  the  elevation  of  a  new  floodplain  within  the  incised  channel, 
a  bank  height  approaching  critical  with  failure  rate  lower  than  in  Type  3  reaches  which  is 
enhanced  by  berms  loading  the  bank  toe,  a  near  trapezoidal  cross-sectional  form,  and  width- 
depth  ratio  values  that  are  higher  than  those  for  Type  3  reaches  because  aggradation  has 
reduced  the  channel  depth  and  bank  failure  has  increased  the  width  of  the  channel. 

•  Type  5  :  These  are  located  downstream  of  Type  4  reaches  and  are  characterised  by  a  near 
equilibrium  between  sediment  transport  capacity  and  sediment  supply  for  the  channel  cross- 
section  that  is  defined  by  the  top  of  the  berms,  a  bank  height  less  than  critical  and  reduced  bank 
angles  due  to  incorporation  of  failed  bank  material  into  berms,  colonisation  of  berms  by 
riparian  vegetation  which  increases  hydraulic  roughness  and  thus  promotes  further  accretion, 
sediment  depth  in  the  bed  of  the  channel  that  generally  exceeds  four  feet,  a  compound  channel 
with  the  low  stage  defined  by  the  tops  of  the  berms  and  high  stage  defined  by  pre-incision 
floodplain,  and  width  to  depth  ratios  that  are  higher  than  those  of  Type  4  reaches.  Type  5 
reaches  define  the  new  dynamic  equilibrium  conditions  in  the  channel. 

This  model  assumes  that  the  channel  is  responding  to  a  single  base  level  lowering  event  and 
that  the  channel  perimeter  is  readily  adjustable  both  vertically  and  laterally.  The  model  is  also 
only  applicable  at  a  system  wide  level  as  local  processes  such  as  meander  bend  erosion  will 
distort  the  sequence.  Nevei-the-less,  the  sequence  provides  a  framework  whereby  the  complex 
adjustments  in  a  channel  can  be  understood  and  the  morphology  of  individual  reaches  related 
to  processes  operating  both  upstream  and  downstream  in  the  network. 

Bank  erosion  in  Type  3  and  4  channels  can  cause  considerable  problems  in  terms  of  loss  of 
agricultural  land,  and  the  undermining  of  bridges  and  other  structures.  Sediment  production 
through  bed  erosion  and  bank  failure  from  the  adjusting  channel  segments  also  causes 
problems  further  downstream,  where  reservoirs  can  become  choked,  navigation  channels 
blocked  and  the  risk  of  overbank  flooding  exacerbated. 


Figure  3.1a  Generalised  Stratigraphic  Section  of  Pieistocene-Holocene  Valley  Fill, 
Northern  Mississippi  (modined  from  Schumm,  1984) 
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Figure  3.1b  Schematic  L.ongitudinal  Profile  of  the  Northern  Mississippi  Valleys  Showing 
the  distribution  of  Pieistocene-Holocene  Stratigraphic  units  (modified  from  Schumm, 
1984) 
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Figure  3.1c  Schematic  Longitudinal  Profile  of  an  Active  Channel  Showing  Identifiable 
Features  (modified  from  Schumm,  1984) 
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Figure  3.2  Geological  map  of  Mississippi  showing  the  DEC  survey  area  (modiFied  from 
Watson  et  al.,  1993) 


Table  3.!  Summary  of  morphometric  data  used  to  determine  channel  evolution  model  - 


Oaklimiter  Creek  (after  Schumm  et  aL,  1984) 


Stage 

L(X'.ation 

fop  Widtli 
(tl) 

l>cpth(n) 

Width;  Depth 
Ratio  (rt) 

■fhalwcg 
Slope  (ll/lt) 

Depth  of 
Sediment  (It) 

Dominant  PriKcs.s 

up.strcam  of 
hcuUcut 

82 

Xl.S 

4.7 

0,0020 

0 

traiisiwrt  of 
.sediment 

II 

IimncUialciy 

dowiistrcain 

ol'heaUcut 

82 

21.0 

.1.8 

0.0018 

variable  0-2 

degradation 

111 

DoMusUcam 
of  11 

100 

20.1 

4.9 

0.0018 

1.5 

rapid  widening 

IV 

Oowistrcatn 
of  III 

115 

19.2 

6.0 

0.0016 

2.5 

aggradation  and 
devt.  of  meandering 
thalweg 

V 

Downstream 
of  IV 

119 

1.5.1 

7.8 

0.0010 

6.1 

aggradation  and 
stabilisation  of 
alternate  bars 

'A' 


Riverbank  retreat  in  type  3  and  4  reaches  takes  place  by  fluvial  erosion  and  bank  mass 
instability  (Thorne,  1981),  The  relative  amount  of  vertical  and  lateral  erosion  (which  cause  an 
increase  in  bank  height  and  bank  angle,  respectively)  is,  therefore,  a  function  of  bank  material 
micro-scale  electrochemical  properties  (in  cohesive  material)  bank  geometry,  bed  materials  and 
excess  flow  boundary  shear-stress  above  a  critical  value  for  entrainment.  Mass  instability  on 
the  other  hand,  depends  on  macro-scale  soil  properties  and  bank  geometry.  Poor  drainage  and 
consequei  t  saturation  promotes  bank  failure  by  causing  positive  pore- water  pressure  which,  in 
turn,  reduces  effective  cohesion  and,  therefore,  the  factor  of  safety  (Fs  =  restraining/disturbing 
force).  The  most  favorable  conditions  for  poor  drainage  are  heavy  precipitation  and/or  rapid 
channel  flow  draw-down,  which  changes  the  bank  from  submerged  to  saturated  conditions 
effectively  doubling  the  bulk  weight  of  bank  material  (Thorne,  1973).  Slab-type  failure  is  most 
common  in  steep  banks  (^60  degrees,  Taylor,  1948),  where  a  block  shears  along  a  planar 
surface  and  topples  forward  into  the  channel.  This  mechanism  is  promoted  by  the  development 
of  tension  cracks  within  the  bank.  The  removal  of  woody  perennial  vegetation  from  the  top  of 
banks  enhances  the  formation  of  these  cracks  althouah  larae  veaetatinn.  .such  as  trees  mav  act 
to  hasten  failure  by  surcharging  (adding  excess  load  to)  the  bank  if  the  root  network  is 
undercut. 

Downstream  of  a  knickpoint  the  channel  widens  by  mass  failure,  followed  by  basal  clean-out  of 
slump  debris,  until  stream  power  has  fallen  to  a  level  where  failed  bank  material  can  no  longer 
be  carried  away  A  sinuous  channel  then  develops,  within  the  confines  of  the  incised  and 
blown-out  terrace  walls  of  the  incised  canyon.  Where  a  meander  encounters  the  terrace  wall, 
bank  failuic  may  continue  to  take  place,  probably  by  a  rotational  slip  of  the  lower  angle,  higher 
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bank.  The  channel  eventually  becomes  stable  once  more,  except  that  it  lies  within  an 
entrenchment  in  the  former  floodplain. 

While  bank  failure  is  related  to  bulk  mechanical  properties,  the  rate  of  retreat  is  dependent 
upon  the  fluvial  entrainment  of  failed  bank  material.  This  phenomenon  has  been  termed  Basal 
Endpoim  Control  (Thome,  1978),  and  it  explains  how  the  rate  of  channel  width  adjustment  is 
controlled  by  the  balance  of  sediment  moving  into  and  out  of  the  basal  area, 
lititially,  slumped  debris  protects  the  bank  from  further  erosion  and  failure  through  buttressing, 
bank  toe  loading,  and  by  buffering  fluvial  activity  away  from  the  bank  base.  Eventually  debris  is 
removed  however,  (it  is  much  less  resistant  to  erosion  than  the  intact  bank)  and  further  basal 
erosion  of  the  bank  takes  place.  Figure  3.3  shows  a  schematic  diagram  of  sediment  fluxes  in 
the  near  bank  zone. 

Figure  3.3  Basal  Endpoint  Control :  Sediment  flux  in  the  near  bank  zone  (modified  from 
Darby  &  Thorne,  1993) 


There  are  three  states  of  basal  endpoint  control,  and  therefore  lateral  channel  stability  : 

1)  INPUT  >  OUTPUT  (Impeded  removal):  Basal  sediment  wedge  accretion,  berm 

formation  and  increasing  bank  stability 

2)  INPUT  -  OUTPUT  (Unimpeded  removal):  Equilibrium  between  processes  delive--;  g 

and  removing  sediment,  bank  retreat 
limited. 

3)  INPUT  <  OUTPUT  (Excess  basal  capacity):  Removal  of  sediment  wedge  by  basal 

scour  causing  basal  lowering,  a  reduction 


in  bank  stability,  and  an  increased  bank 


retreat  rate 

Channel  evolution  models  such  as  that  developed  by  Schumm  et  al.  (1984)  use  the  theory  of 
Basal  Endpoint  Control  to  explain  the  sequence  of  degradation,  widening  and  stabilisation  that 
is  characteristic  of  unstable  Bluff  Line  streams. 

Grissinger  (1982)  recognised  that  bed  stability  is  also  controlled  by  velocity-dependent 
hydraulic  properties  and  describe  how  degradational  erosion  of  the  Lacustrine-Old  Paleosol 
material  is  influenced  by  its  polygonal  structure  and  is  initiated  through  the  development  of  a 
narrow  chute,  which  widens  and  deepens,  with  discrete  blocks  being  eroded  and  entrained 
It  is  channel  bank  failure  downstream  of  a  knickpoint  in  wooded  areas  that  leads  to  the  input  to 
the  channel  of  large  quantities  of  trees  and  other  organic  debris  The  size,  orientation,  and 
persistence  of  organic  debris  then  affects  stream  power,  flow  routing,  the  Basal  Endpoint 
Control  state  and  bank  erosion  rates  It  is,  therefore,  probable  that  the  processes  and 
morphology  in  the  idealised  sequence  of  the  Channel  Evolution  Model  will  be  distorted  by  the 
input  of  substantial  quantities  of  LWD  in  Type  2  and  3  reaches  if  the  riparian  zone  is  heaviliy 
wooded.  Debris  jams  may  either  act  to  reduce  local  channel  instability,  through  energy 
dissipation  and  sediment  retention,  or  may  assist  in  driving  bed  and  bank  erosion  by  causing 
local  bed  scour  and  flow  diversion  into  channel  banks  It  may,  therefore,  be  necessary  to  refine 
the  channel  evolution  model  to  describe  adequately  the  form-process  relationships  occurring  in 
unstable  channels  with  w'ooded  riparian  zones  The  affect  of  LWD  on  the  Channel  Evolution 
Model  is  an  interesting  topic  and  one  which  is  addressed  in  the  results  analysis  of  this  report 


4  RESEARCH  METHODS 

Field  ata  collection  v  as  undertaken  in  May  1995  and  August  1996.  Debris  jams  were  surveyed 
with  long-profiles  and  cross-sections  taken  at  23  reaches  which  are  monitored  by  the  U  S. 
Army  Corps  of  Engineers  as  part  of  the  DEC  monitonng  research  program.  A  full  description 
of  the  DEC  monitoring  site  characteristics  is  given  in  Watson  et  al  (1993).  The  reaches 
surveyed  were  between  1220  and  3660  metres  long  and  ranged  in  upstream  watershed  area 
from  9  to  388  Figure  4.2  shows  the  rivers  and  major  catchments  of  the  project  area  in  detail. 
Reacheswere  surveyed  on  the  following  creeks: 


Nolehoe  Creek 
Burney  Branch 
Sykes  Creek 
Worsham  Creek  (East) 
Abiaca  Creek 
Otoucalofa  Creek 
Perry  Creek 


Lick  Creek 
Long  Creek 
Fannegusha  Creek 
Worsham  Creek  (West) 
Red  Banks  Creek 
Lee  Creek 
Marcum  Creek 


Sarter  Cree>' 

James  Wolf  Creek 
Hotopha  Creek 
Worsham  Creek  (Middle) 

Harland  Creek 
Coila  Creek 
Hickahala  Creek 

These  surveys  provide  a  comprehensive  data-set,  which  not  only  covers  a  range  of  drainage 
basiii  areas,  from  9.5  to  256  Km'^,  but  also  allows  comparison  of  debris  loadings  between 
reaches  with  wooded  and  agricultural  riparian  zones,  between  straight  and  meandering  reaches 
and  between  highly  unstable  and  stabilising  or  equilibrium  reaches. 

Debris  jam  sites  have  been  surveyed  into  the  thalweg  and  cross-section  data  for  each  creek  so 
tliat  there  position  and  associated  changes  in  local  channel  geomorphology  can  be  monitored 
over  time.  Data  from  the  current  survey  has  been  processed,  overlaid  and  compared  with  that 
collected  in  May  1994  so  that  an  assessment  can  be  made  of  the  rate  of  debris  input,  the 
longevity  of  jams  and,  therefore,  their  effectiveness  as  geomorphological  channel  controls,  and 
the  changing  patt  erns  of  associated  sedimentation  and  erosion 

Debris  jam  sediment  budgets  have  also  been  calculated  to  determine  whether  the  net  impact  of 
debris  jams  is  through  sediment  retention  or  sediment  scour  and  mobilisation.  Sedimentation 
processes  are  assessed  using  hydraulic  geometry  relationships.  V  olume  of  debris  and  number  of 
jams,  per  unit  •  ach  length  are  also  related  to  processes  occurring  in  the  sequence  of  stages  in 
the  Schumm  channel  evolution  model  (CEM). 

Geomorphological  ri  onnaissance  was  also  performeo  at  each  jam  .site  to  document  the 
volume  of  debris  in  each  jam,  to  identify  its  mode  of  input  into  the  channel,  to  determine  the 
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jam  type  in  terms  of  impact  upon  flow  pattern  and  erosion,  and  to  measure  the  volume  of 
sediment  retained  in  backwater  areas  or  bars.  The  following  variables  were  assessed  at  each 
jam  site: 

1)  Debris  jam  volume:  Estimated  vol  ime  of  woody  material  (m’)  in  each  jam.  values  are  then 
summed  for  each  survey  reach 

2)  Morphological  classification:  A  debris  classification  system,  modified  from  that  developed 
by  Robinson  and  Beschta  (1989),  wbdeh  describes  the  geoniorphological  impact  of  debris  jams 
throughout  the  drainage  network.  Figure  4  2  shows  the  original  classification  scheme  and 
figure  5,2a  shows  the  modified  version  Jam  classification  types  are  as  follows: 

Underflow  jams.  In  small  catchments  where  fallen  trees  span  the  channel  at  bank-full 
level  Local  bed  scour  may  occur  under  debris  at  high  flows,  otherwise  the  in¬ 
channel  geomorphic  impact  of  the  LWD  is  minimal. 

Dam  jams  In  channels  where  the  average  tree  height  to  channel  width  ratio  is  roughly 
equal  to  one,  so  that  debris  completely  spans  the  channel  cross-section.  This 
type  of  jam  causes  significant  local  bank  erosion  and  bed  scour  due  to  flow 
constriction  and  backwater  effects  will  cause  sediment  deposition  in  the  lower 
energy  environment  upstream  Bars  may  also  form  immediately  downstream  of 
the  jam. 

Deflector  jams.  Found  where  input  debris  does  not  quite  span  the  channel  so  that  flow 
is  deflected  against  one  or  both  banks  causing  localised  bed  scour  and  bank 
erosion.  Subsequent  bank  failures  result  in  the  input  of  new  LWD  material  to 
the  reach,  so  that  the  jam  builds  up  further  Backwater  sediment  wedges  and 
downstream  bars  may  form  at  this  type  of  jam  provided  that  stream  power  is 
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suspended  sediment  transport. 

Plow  Parallel  jams.  Found  where  channel  width  is  significantly  greater  than  the  key-debris 
length  and  flows  are  competent  to  rotate  debris  so  that  it  lies  parallel  to  the 
flow.  Debris  is  also  transported  downstream  in  high  flows  and  deposited  against 
the  bank-base  on  the  outside  of  meander  bends  and  at  channel  obstructions  such 
as  hydraulic  structures.  Related  bank  erosion  and  bed  scour  will  be  minimal,  and 
bank  toes  may  even  be  stabilised  by  debris  build-up.  Flow  parallel  debris  may 
also  initiate  or  accelerate  the  formation  of  mid-channel  and  lateral  bars 
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3)  Blockage  classification:  Jams  are  classified,  using  the  scheme  developed  by  Gregory  et  al, 
(1985),  according  to  their  potential  to  block  the  downstream  movement  of  water  and  sediment 
The  classification  types  are: 

Active:  Jam  forms  a  complete  barrier  to  water  and  sediment  movement  and  also  creates 
a  distinct  step  or  fall  in  the  channel  profile 
Complete,  complete  barrier  to  water/sediment  movement,  but  no  significant  step 
Partial:  Jam  is  only  a  partial  barrier  to  flow. 

4)  Alpha/Beta  Indices:  The  alpha  angle  describes  the  predominant  alignment  of  the  debris  in 
the  channel  with  respect  to  the  flow  direction.  These  indices  were  first  used  by  Cherry  and 
Beschta  (1986)  in  connection  with  LWD  flume  experiments.  The  Beta  angle  is  a  measure  of 
the  predominant  orientation  of  the  debris  jam  material  in  the  vertical  plane,  orthogonal  to  the 
flow  direction.  See  figure  4  1  below 

Figure  4.1  Alpha  and  Beta  Indices  (modified  from  Cherry  and  Beschta,  1986) 
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Alpha  Index  (degrees) 


Beta  Index  (degrees) 


5)  Sinuosity:  A  vi,sual  estimate  of  channel  sinuosity  in  the  jam  reach  (straight,  slightly  sinuous, 
sinuous,  meandering) 

6)  Knickzones:  Presence  of  knickpoints  or  knickzones  in  the  channel  reach  (a  measure  of 
channel  instability). 

7)  Sediment:  Bedload  D;o  classification  (clay,  silt,  sand,  gravel) 
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8)  Deposition/Scour:  Estimated  volume  of  bar  deposition,  backwater  sedimentation  and 
bed^ank  scour  induced  by  each  jam,  Total  deposition  and  scour  values  are  then  calculated  for 
each  channel  reach. 

The  debris  volume  and  debris  frequency  measures,  morphological  classification,  jam 
sedimentation  and  erosion,  and  net  jam  sediment  budgets  in  each  reach  have  been  and  plotted 
against  three  independent  catchment  variables.  These  are:  upstream  drainage  basin  area,  reach- 
averaged  channel  width,  and  reach-average  stream  power  per  unit  bed  area.  Stream  power  per 
unit  bed  areais  calculated  using  the  following  equation: 


w 


(3  1) 


where,  co  =  stream  power  per  unit  channel  channel  length  (N/m/s  '),  p=  density  of  water 
(kg/m^),  g  =  gravitational  constant  (9.81  m/s'^),  Q  =  predicted  bankfull  discharge 
(cumecs),  s  =  bed  slope  (m/m),  w  -  reach  average  channel  width  (m). 

These  independent  catchment  variables  are  used  to  determine  whether  the  geomorphological 
effects  of  LWD  have  a  coherent  and  predictable  spatial  relationship. 
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Figure  4.2  DEC  project  area  site  location  map  (modified  from  Cooper  et  al,  1996) 


5  RESULTS  AND  DISCUSSION 


All  raw  field  data  and  calculations  are  listed  in  Appendix  D 

5.1  Input  Mechanisms 

Debris  jams  were  found  in  17  of  the  23  reaches  surveyed.  All  reaches  containing  jams  were 
found  to  have  wooded  riparian  zones,  while  the  6  reaches  with  no  significant  debris 
accumulations  were  found  to  have  non-wooded,  open  riparian  zones. 

‘Key  debris’  (that  is  large  trees  which  initiate  jam  formation  (Nakamura  and  Swanson,  1993)) 
input  mechanisms  were  found  to  fall  into  the  following  categories: 

37%  due  to  outer  bank  erosion  in  channel  bends; 

36%  due  to  bank  mass-wasting  in  degrading  reaches; 

1 2%  due  to  wind-throw, 

7%  resulting  from  ‘paleodebris’  (material  reintroduced  into  the  channel  from  erosion  of 
old  alluvial  deposits  containing  preserved  debris); 

5%  initiated  by  large  logs  floated  from  upstream;  and 
4%  were  found  to  be  formed  by  beaver  dams. 

Plate  1  shows  debris  input  due  to  mass  failure  of  the  bank  on  Harland  Creek  and  Plate  2  shows 
a  beaver-LWD  “dam”  located  on  the  sill  of  a  low  drop  grade  control  structure  on  Worsham 
Creek. 

If  channels  have  wooded  riparian  zones,  the  largc.st  volume  of  debris  input  can  be  predicted  to 
occur  where  the  channel  is  in  phase  two  and  three  of  the  CEM  or  is  actively  meandering.  Jams 
tend  to  form  where  the  key  debris  elements  fall  into  the  river  and.  hence,  arc  commonly  located 
at  bend  apices  or  in  unstable  i  caches  downstream  of  knickpoints. 

5.2  Distribution  of  LWD 

Correlations  were  made  between  average  volume  of  debris  (per  vol.  Per  100m)  measured  in 
each  reach  sup/eyed  against  watershed  area  and  stream  power.  Positive  relationships  were 
expected  in  these  correlations,  as  found  by  Gregory  et  al.  (1993),  because  debris  is  more  likely 
to  be  mobilised  and  transported  downstream  into  larger  channels  as  flow  competence 
increases.  In  fact,  no  statistically  significant  relationship  was  found  between  watershed  area 
and  debris  volume.  However,  the  relationship  between  stream  power  and  debris  volume  was 


found  to  be  statistically  significant  using  Pearson  Product  Moment  correlation  analysis,  at  the 
0,05  confidence  level.  This  relationship  is  shown  in  Figure  5.  la. 

Correlations  were  made  between  average  number  of  jams  (jam  frequency)  per  lOOni  reach 
length  for  each  reach  surveyed,  watershed  area  and  stream  power  Negative  relationships  were 
expected  in  these  plots  because  as  channel  size  becomes  greater  flow  competence  increases 
and  jams  are,  therefore,  more  likely  to  become  dislodged  (Swanson  et  al ,  1984)  However, 
neither  correlation  was  found  to  have  a  statistically  significant  relationship. 

The  use  of  hydraulic  geometry  and  energy  relationships  to  predict  the  distribution  of  debris  is, 
therefore,  simple  an  approach  to  explain  debris  dynamics  in  these  unstable  channel 
environments.  Reach  stability,  and  channel  sinuosity  are  probably  better  predictors  of  debris 
volume  and  the  frequency  of  jams,  because  these  factors  to  a  large  extent  determine  the  rate  of 
debris  input, 

Examination  of  the  values  for  debris  loadings  for  each  Stage  (1  to  5)  of  reach  stability  outlined 
in  the  Schumm’s  (1984)  Channel  Evolution  Model  were  examined  to  detennine  whether 
loadings  can  be  attributed  to  the  processes  occurring  in  the  stages  and,  therefore,  whether 
stage  of  evolution  is  a  good  measure  of  debris  input  rate.  The  CEM  model  is  shown  in  Figure 
3.1c. 

It  is  expected  that  there  would  be  very  low  debris  input  in  stable.  Stage  I  reaches,  where  n''\ss 
failure  of  the  bank  is  absent  or  limited.  A  rapid  rise  in  debris  input  rate  might  be  expected  in 
Stage  2,  reaches,  where  knickpoints  are  active,  because  although  the  channel  has  not  yet  begun 
to  adjust  laterally,  large  trees  on  the  incised  channel  banks  produce  surcharging  which  could 
result  in  local  bank  failures.  It  is  expected  that  there  would  be  the  highest  debris  input  rates  in 
Stage  3  reaches  due  to  widespre  id  mass  failure  of  the  banks  following  degradation  Finally,  in 
Stage  4  reaches,  where  bank  erosion  rates  decline  and  in  Stage  5  reaches,  where  bank  erosion 
ceases,  debris  input  rate  and,  therefore,  in-channel  debris  volumes  would  also  be  expected  to 
decline.  Debris  will  also  become  buried  due  to  the  increasing  rate  of  alluvial  deposition  in 
Stage  4  and  5  reaches.  For  the  same  reasons  a  similar  distribution  would  also  be  expected  for 
the  frequency  of  jams  in  each  reach 

Averaging  the  volume  and  frequency  data  for  the  survey  reaches  found  in  each  stag/;  produces 
the  distribution  of  values  shown  in  Table  5.1  (values  arc  weighted  according  to  the  length  of 
channel  in  each  stage).  It  is  evident  that  the  trends  through  the  values  closely  match  those 
described  above.  The  CEM  categories,  therefore,  appear  to  be  a  better  indicator  of  the  relative 
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magnitudes  of  debris  volume  and  jam  frequency  found  in  various  channel  reaches.  The 
percentage  differences  shown  in  T>ible  5. 1  are  relative  to  ‘equilibrium’,  Stage  1  reaches. 

Outer  bank  erosion  in  active  meanders  was  found  to  be  the  dominant  mechanism  for  debris 
input.  It  might,  therefore  he  assumed  that  channel  sinuosity  would  be  a  good  predictor  of 
debris  volume  and  the  rmmber  of  debris  jams  per  unit  reach  !eng;h.  In  fact  regression  analysis 
of  channel  sinuosity  against  average  volume  of  debris  and  average  number  of  jams  for  all 
reaches  surveyed  produced  r^  values  of  only  0.05  and  0.02  respectively,  when  power  functions 
were  fitted,  I'tiese  coefficients  of  determination  are  not  statistically  significant. 

While  the  trend  of  the  functions  is  as  expected,  with  more  debris  and  a  greater  number  of  jams 
per  unit  reach  length  as  sinuosity  increases,  the  weakness  of  the  relationships  means  that  they 
have  no  predictive  value.  The  explanation  probably  lies  in  the  fact  that  a  high  sinuosity  does 
not  necessarily  imply  active  meandering  (ie  where  the  channel  has  a  sinuous  planform  which 
actively  adjusts  through  outer  bank  erosion  and  point  bar  deposition).  Active  channel  migration 
is  difficult  to  quantify  and  requires  knowledge  of  other  parameters  such  as  available  stream 
power  and  bank  material  properties,  or  observation  of  channel  planform  change  over  time. 
There  are  insufficient  data  available  at  present  to  incorporate  a  reliable  active  meandering 
factor  into  the  analysis. 


Table  S.l  :  Channel  Evolution  Model  Stage  as  a  measure  of  lii-channel  Debris  Volume 
and  Frequency  of  Debris  Jams _ _ _ _ 


CEM 

Stage 

no.  of  reaches 
surs'cycd 

reach  length 
surveyed  (m) 

\  ol.  of  debris:  av.  per 
100m  sub-reach  (ni’) 

%  din'ercnce 
from  Stage  1 

no.  of  Jams:  av.  per 
lOOni  sub-reaeh 

%  diiTercnce 
from  Stage  1 

1 

3 

16000 

.S.4 

- 

0.22 

- 

2 

5 

17000 

6.1 

t  13 

0,48 

+  1 18 

4 

22(K)0 

V.l 

Ht') 

0,41 

+86 

4 

3 

13000 

6,7 

+24 

0.38 

+72 

2 

12000 

o.« 

-85 

0,12 

-45 

S.3  Impact  of  LWD  >n  Channel  Processes 

The  impact  of  debris  jams  upon  local  scour  and  sedimentation  processes  has  been  analysed  in 
stable,  gravel-bed  rivers  (Abbe  and  Montgomery,  19%)  but  has  not  been  studied  in  unstable 
sand  bed  rivers  Two  important  questions  that  need  to  be  answered  in  the  case  of  unstable, 
sand-bed  rivers,  are  whether  debris  jams  cause  net  scour  or  net  sediment  retention  and  whether 
there  are  any  spatial  trends  in  these  processes.  Answering  these  questions  will  enable  river 
resti  nation  programs  to  determine  whether  or  not  debris  is  of  benefit  to  incising  sand-bed 
ri  I  S,  in  terms  of  aiding  stabilisation  and  sediment  retention  and  enhancing  aquatic  habitat. 


Geomorphic  field  reconnaissance  in  the  current  research  and  in  previous  studies  (Wallerstein 
and  Thorne,  1994)  indicates  that  the  ratio  of  key  debris  length  to  average  channel  width  and, 
therefore  to  some  extent  watershed  area  (using  hydraulic  geometry  principles),  is  a  good 
indicator  of  the  impact  that  a  jam  will  have  upon  flow  field  and  hence  channel  morphology  and 
sedimentation  processes.  A  debris  classification  scheme,  modified  from  a  pool  formation  model 
(Robinson  and  Beschta,  1990)  has  been  used  to  describe  the  observed  impact  of  debris  jams 
when  moving  downstream  through  the  watershed  network.  This  classification  scheme  is 
outlined  in  Figure  5.2a.  Plate  3  shows  an  example  of  an  underflow  type  jam  on  Worsham 
Creek  (watershed  area  10.3  km^);  Plate  4  shows  an  example  of  a  0am  type  jam  on  Lick  Creek 
(watershed  area  22  km^),  Plate  5  shows  an  example  of  a  Deflector  type  jam  on  Famregusha 
Creek  (watershed  area  46  km^);  and  plate  6  shows  an  example  of  a  Flow  Parallel  type  jam  on 
Harland  Creek  (watershed  area  69  km^)  The  observed  frequency  distribution  of  jam  type  for 
each  reach  is  shown  in  Table  5.2  Watershed  area  for  each  reach  is  also  shown  to  demonstrate 
the  spatial  trend  in  jam  type. 

It  is  evident  that  the  jam  type  encountered  with  increasing  watershed  area  matches  to  a 
reasonable  degree,  the  distribution  outlined  in  the  simple  schematic  model  in  Figure  5.2a 
Anomalies  in  the  distribution  in  Table  5.2  can,  for  the  most  part,  be  explained  by  variations  in 
local  channel  stability  which  affect  channel  dimensions  and  distort  simple  hydraulic  geometry 
relationships.  For  the  simple  spatial  model  to  be  a  good  predictor  of  jam  type  it  must, 
therefore,  incorporate  other  predictive  variables.  A  channel  stability  factor  is  necessary  to 
account  for  unstable  channel  geometries  that  cannot  be  predicted  by  stable  hydraulic  geometry 
relationships.  A  further  factor  is  also  necessaiy  to  account  for  the  degree  of  active  channel 
migration.  These  factors  are  to  be  studied  in  future  research. 

Given  that  there  is  a  weak  spatial  trend  in  jam  type  with  distance  dov/nstreain  in  the  watershed, 
it  is  also  reasonable  to  hypothesise  that  local  scour  and  sediment  retention  at  debris  jams  is  also 
spatially  predictable  Sediment  retention  and  local  scour  caused  by  debris  jams,  averaged  over 
a  100m  channel  length,  are  plotted  against  watershed  area  in  Figure  5.1b. 

The  trends  in  Figure  5  lb  are  best  described  by  quadratic  regression  relationships.  These  trends 
are  statistically  significant  using  one-tailed  Pearson  Product  Moment  tests,  at  the  0.05 
confidence  level.  No  statistically  significant  trends  were  found  when  stream  power  was 
correlated  with  sediment  retention  and  scour 

The  trends  in  Figure  5.1b  can  be  explained  because  they  mirror  the  impact  of  each  debris  jam 
classification  type.  Underflow  jams,  in  small  watersheds,  interfere  very  little  with  the  flow  and, 


therefore,  do  not  have  a  high  scour  and  sediment  retention  potential.  Dam  type  jams,  found 
further  dov/nstream,  cause  large  volumes  of  sediment  to  be  stored  in  backwaters,  but  also 
cause  plunge  pool  scour.  Downstream  from  Dam  type  jams,  at  Deflector  jams,  sediment  is 
stored  in  the  lee  of  jams  as  bar  deposits,  but  the  jams  also  cause  flow  to  impinge  on  one  or 
both  banks,  resulting  in  scour  and  bank  collapse.  Further  downstream  still,  Flow-Parallel  debris 
block  the  flow  much  less,  so  that  energy  dissipation,  sediment  retention  and  scour,  are  less 
effective. 

The  net  impact  of  debris  jams  on  the  sediment  budget  in  each  reach,  calculated  as  volume 
stored  minus  volume  scoured,  is  shown  in  Table  5.2.  These  data  were  plotted  against 
watershed  area  to  deteimine  whether  there  were  any  spatial  trends  between  these  variables. 
No,  statistically  signiucant  relationship  were  found.  It  is  important  to  note,  however,  that  the 
data  indicate  that  the  balance  between  sediment  scour  and  sediment  retention  caused  by  debiis 
jams  favors  of  net  sedimentation.  Thirteen  of  the  channel  reaches  have  a  positive  budget,  one 
has  a  zero  budget  and  only  two  have  negative  sediment  storage  budgets. 

5.4  Jam  Stability  Through  1'ime 

Comparison  of  survey  results  between  summer  1995  and  summer  1996  shows  there  to  be  little 
change  in  the  number  and  position  of  debris  jams  present  in  each  reach  A  total  of  99  jams 
were  surveyed,  in  17  of  the  23  reaches  established  in  1995.  The  six  other  reaches,  all  of  which 
had  agricultural  riparian  zones,  were  found  to  contain  no  significant  debris  accumulations  By 
summer  1996,  4  jams  had  been  destroyed,  of  which  one  had  been  buried  by  sediment  in  a  stage 
4  reach,  one  had  been  manually  removed  for  unknown  reasons,  and  two  had  been  broken  up 
and  transported  downstream  by  the  flow.  Six  new  jams  were  identified  in  the  survey  of  1 996, 
two  of  which  were  new  beaver  dams,  while  the  other  four  were  found  to  have  formed  by  input 
of 'key  '  debris  elements  due  to  bank  mass  failure 
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Figure  5.1a:  Volume  of  Debris  as  a  function  of  Stream  Power 
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Figure  4.2  Debris  Jam  Classification  Model 
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5.5  DEBRIS  JAM  RESIDENCE  TIMES:  A  THEORETICAL  ANALYSIS 


If  is  apparent  from  the  preceding  analysis  of  debris  jam  stability  between  1994  and  1996  that, 
despite  the  fact  that  many  of  the  channels  are  highly  unstable  and  that  the  discharge  regime  is 
also  very  flashy,  the  majority  of  jams,  have  survived  intack  and  remain  at  the  location  where 
they  were  first  surveyed.  It  should  be  noted  that  several  storm  runoff  events  during  the  period 
1993-94  in  each  of  the  three  creeks  exceeded  the  2-year  discharge  values  calculated  by  Watson 
et.  al,  (1993).  If  these  events  are  taken  as  representative  of  an  average  year  in  the  region  it 
follows  that  debris  jams  are  subjected  tc  greater  than  channel-ibrming  flows  several  time  per 
year  and  yet  remain  in  place  This  finding  can  be  verified  in  a  more  analytical  fashion  by 
examining  the  buoyancy  and  drag  forces  which  act  upon  elements  comprising  a  debris  jam. 

For  a  debris  element  to  be  lifted  fiom  the  bed  and  floated  downstream  the  buoyant  force  acting 
on  the  element  must  exceed  the  weight  of  the  element. 

The  flotation  force  due  to  the  pressure  on  the  under  surface  of  a  submerged  or  partially 
submerged  body  is  given  by: 


Ff-  pwgLA  (5.1) 

where  ,  Ff  -  flotation  force  (N),  pw  =  density  of  water  (1000  kg/m'),  g  =  gravitational  constant 
(9,81  m/s'^),  LA  =  volume  of  the  body  submerged  in  water. 

The  weight  force  resisting  flotation  is  given  by; 

W,  =  p,gla  (5,2) 

where,  W,  ^  weight  (N),  p,  =  density  of  the  body  (kg/m'),  la  =  total  volume  of  the  body  (m') 
To  bring  a  object  to  the  point  of  floating,  these  forces  must  be  balanced: 

Ff  -  W,  01 

p,vgLA/pJa-l  (5.3) 

The  density  of  wood  varies  greatly  and  values  between  400  and  800  kg/m'  are  suggested  by 
Ashby  &  Jones  (1980).  The  density  of  debris  in  a  river  vriH  vary  as  a  function  of  time  and  How 
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history  md  as  for  example  the  wood  becomes  water-logged  and  then  dries  out.  Hence,  it  is 
difficult  to  ascribe  a  precise  density  to  debris  jams  per  sg.  For  the  purposes  of  this  analysis  a 
value  of  800  kg/m’  is  used  on  the  assumption  that  the  debris  element  has  recently  fallen  into 
the  river  and  is,  therefore,  not  waterlogged. 

The  volume  of  the  debris  element  in  a  jam  varies  greatly  with  time,  but  of  greatest  importance 
are  the  “key”  debris  pieces,  which  initiate  the  formation  of  jams  and  are  largely  responsible  for 
their  stability.  An  “average”  large  tree  trunk  in  the  region  is  of  the  order  of  25m  long  with  a 
maximum  diameter  of  Im  yielding  a  volume  of  25m\  this  of  course  ignores  the  volume  (and 
added  density)  of  the  root  wad  and  of  the  smaller  branches  and  leaves.  For  these  para-neters 
the  flow  depth  required  to  float  the  element  is; 
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1000x9.81xLA  7  800x9.81x25  -  1  =9810xLA=  196,200 
hence,  LA  =  20m’,  and  20/plan  area  (25m^)  =  0.8m  water  depth. 

This  depth  of  flow  is  shallower  than  that  occurring  at  all  the  creek  sites  at  flow  and  so,  in 
theory,  even  quite  large  t'ees  could  be  floated  several  times  per  year.  However,  the  filed 
evidence  evidently  indicates  that  this  does  not  occur  actually  that  often.  Factors  reducing  the 
potential  for  tree  floatation  include,  as  already  mentioned,  increased  wood  density  due  to 
water-logging,  the  great  weight  of  the  root  wad,  grounding  of  trees  on  shallows,  and  snagging 
of  the  root  wad  and  branches  on  the  bank  and  other  vegetation. 

The  second  motivatinng  force  to  consider  is  drag  force  on  debris  Abbe  and  Montgomer>' 
(1996)  use  an  analytical  formula  for  describing  flow  past  a  bridge  pier  as  an  analogue  to 
produce  first  order  approximations  for  the  force  acting  on  debris  elements.  This  formula  is 
given  as  (see  Chadwick  &  Morfett,  1986): 


^  F„  =  0  5p,vCnlfA  (5.4) 

where,  Fd  =  drag  force,  pw  =  density  of  water  (1000  kg/m’),  Cp  =  drag  coefficient,  U  =  mean 
incident  flow  velocity,  A  =  submerged  area  of  the  obstruction,  normal  to  the  incident 
^  flow. 

Selection  of  a  drag  coefficient  value  is  problematical.  Petryk  &  Bosmajia,  (1975)  use  (  d  =  1, 
for  flow  through  flexible  vegetation.  Boundary  conditions  have  the  effect  of  increasing  Cd,  as 
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does  the  blockage  ratio,  Br,  defined  as  the  ratio  of  the  width  of  the  obstruction  to  the  channel 
width,  and  Gippe!  et  al.  (1992)  suggest  that  the  effect  on  Cd  becomes  significant  when  Br  » 
0.05.  Conversely,  surface  roughness  on  the  debris  increases  turbulence  and  thus  reduces  drag. 
Abbe  &  Montgomeiy  use  an  average  value  of  Cd  =  1.55,  as  suggested  by  Rouse  (1946),  and 
this  value  will  be  used  in  the  following  analysis. 

If  a  debris  element  is  to  be  moved,  the  drag  force  must  exceed  frictional  resistance  provided  by 
the  immersed  normal  weight  of  the  tree  and  contact  with  the  channel  bed  Interlocking  between 
debris  elements  also  greatly  increases  the  resisting  forces. 

For  the  sake  of  simplicity  the  forces  on  a  simple  “key”  debris  element  will  be  assessed,  having  a 
length  of  25m,  diameter  of  Im,  and  average  density  of  800kg/m\  Interlocking  forces  are  not 
examined  due  to  the  complexity  of  calculation,  and  the  following  analysis  therefore  provides  a 
“worst  case”  analysis  of  force  balance. 

The  drag  force  is  determined  assuming  that  the  element  lies  in  the  channel  with  long  axis  at 
orthogonal  to  the  flow,  as  has  been  found  at  most  debris  jam  sites  in  the  smaller  creeks.  Hence, 
it  presents  an  area  of  25m"  normal  to  the  incident  flow  It  is  worth  noting  that  even  if  we 
consider  trees  in  a  larger  creek  which  have  been  rotated  parallel  to  the  flow  and  have  a  root 
wad  normal  to  the  incident  flow,  the  area  of  this  root  wad  will  rarely  exceed  25m^  and,  in  any 
case,  it  is  unlikely  that  the  entire  root  wad  would  be  submerged,  even  during  over-bank  flows. 
The  assumptions  made  in  this  analysis  should  be  tenable  for  large  debris  elements,  whatever 
their  orientation  to  the  flow  may  be.  The  drag  force  on  such  an  element  has  been  calculated 
using  mean  flow  velocity  for  the  predicted  2-ycar  discharge  (02),  (or  bankfull  discharge  if  Q2 
is  an  out  of  bank  flow),  for  each  on  the  creek  reaches,  where  data  are  available.  The  results  are 
presented  in  Table  5  3. 

Frictional  resisting  forces  can  be  approximated  to: 


4  F,  =  Wstani|>  (5.5) 

where,  F,  =  frictional  resisting  force  (N)  W,  =  submerged  weight  of  the  log  (N),  <|)  -  angle  of 

.nter.nal  fri‘  ion  for  the  bed  material  (degrees) 

i 

The  angle  of  internal  friction  for  the  bed  material  can  be  estimated  from  analysis  presented  in 
civil  engineering  texts  such  as  Ashby  and  .lones  (1996).  The  submerged  weight  of  the  log, 
based  on  the  previous  analysis  of  buo>ancy  force  is  however  less  than  zero,  making  the 
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estimation  of  a  resisting  force  by  this  method  impossible.  It  is  beyond  the  limits  of  this  study  to 
njake  a  more  in-depth  analysis  of  force  balance  on  LWD  elements,  however,  from  the 
observation  of  debris  jams  over  the  thr  ee  year  study  period  it  can  be  concluded  that  jams  which 
are  initiated  by  large  key  debris  elements,  fonn  semi  permenant  structures  in  the  channel  and 
are  therefore  significant  agents  of  local  channel  geomorphological  genesis  and  change. 

Table  5.3  Drag  Force  on  Simple  Debris  Elements 


Creek 

Mean  flow  velocity 
at  bankfull  flow  (m/s 

Drag  Force  (kN) 

Hurland 

1.19 

27.55 

Paniicgusha 

1.85 

66.64 

Abiaca  3 

1.42 

39.25 

Abiaca  4 

1.22 

28.95 

Coila 

1.22 

28.70 

Abiaca  6 

1.04 

20.82 

Nolchoc 

1.98 

75.92 

Red  Banks 

1.77 

00.46 

Lee 

0.40 

3.17 

Hickalrala 

1.66 

53.62 

Lower  Hotoplia 

1.25 

30.45 

Upper  Holopha 

0.95 

17.63 

Marcum 

1.82 

64.46 

CXoucaloi'a 

1.84 

65.42 

Sartor 

1.58 

48.59 

Pcrt> 

1.2.3 

29.53 

Sykes 

0.87 

14.72 

East  Worsham 

169 

.55.27 

Middle  Worsham 

1.50 

40  96 

West  Worsham 

1.85 

66  10 

James  Wolf 

1.52 

44  47 

Long 

1  27 

31  22 

6  LARGE  WOODY  DEBRIS  MANAGEMENT  PROGRAM  (DMP) 

The  relationships  between  LWD  and  channel  processes  have  been  incorporated  into  a  LWl) 
Management  Program  (DMP).  Version  2.0  of  this  program  is  included  with  this  report  on  a 
disk. 

Version  2.0  is  an  updated  version  of  the  program  to  that  included  with  the  Project  R&D  7258- 
EN-09,  submitted  to  the  US  Army  Corps  of  Engineers,  June  1995  (Wallerstein,  1995). 

The  program  predicts  the  likely  jam  type  in  a  given  reaeh,  determines  its  impact  upon  the 
channel  and  outlines  an  appropriate  management  strategy.  Inputs  variables  are  those  found  to 
be  most  critical  in  the  analysis  of  LWD  in  the  fluvial  system  and  include  channel  width 
(determined  from  a  catchment  area  function),  average  riparian  tree  height,  reach  sediment  type 
and  the  ripatian  land-use  type.  The  ratio  of  tree  height  to  channel  width  is  used  to  define  the 
debris  jam  type  likely  to  be  present,  with  the  precise  limits  of  each  classification  determined 
Irom  the  empirical  relationships.  Sediment  diameter  is  used  to  give  an  indication  of  the  jam’s 
potential  to  induce  backwater  sedimentation  or  downstream  bars  Debris  jam  types  are 
classified  using  a  scheme  modified  from  Robinson  &  Bcschta  (1990),  described  in  Wallerstein 
&  Thorne  (1994)  Jani  types  are  divided  into  llnderflow.  Dam,  Deflector  and  Plow  Parallel. 
Program  output  takes  the  form  of  a  text  file  which  describes  indicated  likely  jam-type  and 
olVers  advice  on  appropriate  in-channel  LWD  management  strategics  for  this  situation.  While 
the  management  strategies  are  based  solely  on  theoretical  considerations,  the  program  never- 
the-less  provides  a  framework  for  future  model  development  as  empirical  relationships 
between  the  variables  are  better  characterised.  A  flow  diagram  of  the  computer  program  is 
shown  in  figure  b  I . 


6.1  GIS  Interface 

The  program  has  also  been  linked  to  a  GIS  (Geographical  Information  System)  data  input 
system  which  was  constructed  by  Peter  Cheeseman,  a  masters  student  at  Nottingham 
University  as  part  of  this  project  (Cheesman,  1995).  The  work  was  carried  out  to  demonstrate 
the  potential  for  using  GIS  as  a  platform  for  data  input  to  expert  systems  designed  to  aid 
engineers  with  river  basin  management. 

The  GIS  was  constructed  in  ARC  INFO  using  data  layers,  supplied  by  the  WES  Intergraph 
data-base,  for  the  Abiaca  Creek  watershed  and  provides  both  automatic  data  input  for  the 
necessary  variables  and  a  platform  for  running  the  program.  This  walershed  wa.s  selected 
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because  it  contains  four  of  the  debris  survey  reaches  being  monitored  in  the  current  research. 
The  theoretical  model  could,  therefore,  be  tested  against  the  empirical  data  results  from  the 
field  studies,  validated  and  further  developed. 

A  full  description  of  the  (ilS  model  construction  is  presented  in  Appendix  A. 


The  GIS  user  interface  for  the  DMP  is  also  now  available  for  use  and  can  be  obtained 
upon  request  from  nick@geography.nuttinghatn.ac.uk 

The  program  runs  on  the  UNIX  version  of  ARC  INFO  and  requires  approximately  98  mb  of 
memory  to  install. 


The  GIS  is  composed  of  four  layers:  the  drainage  network;  road  network;  landcover;  and, 
channel  sediment  type.  There  is  also  a  terrain  model  which  is  used  to  calculate  drainage  basin 
area.  The  system  incorporates  a  menu  driven  interface  used  to  display  the  data  layers  and 
perform  analyses  On-line  help  files  are  also  included.  The  analysis  is  performed  by  simply 
placing  the  cursor  on  the  area  of  interest  and  clicking  the  mouse.  The  system  then  e.xtracts  the 
relev  ant  input  data  from  the  database  for  that  location  and  passes  them  to  an  input  file.  The 
debris  management  program  is  then  automatically  activated  and  reads  the  input  file,  calculates 
the  results  and  produces  an  output  text  file.  This  file  is  then  read  back  into  the  GIS  and 
displayed  on  the  computer  screen. 

This  management  model  is  simple  to  operate  and  provides  a  framework  for  future  development 
as  empirical  relationships  between  variables  are  better  characterised. 
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Figure  6.1  LWD  management  program  flow  diagram 
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6.2  Testing  the  Model 

The  model  has  been  tested  against  the  geomorphic  data  obtained  for  each  reach  in  northern 
Mississippi  In  the  test  run  channel  width  was  determined  using  the  linear  regression  function 
obtained  using  the  DEC  data  from  Watson  et.  al.  (1993).  The  formula  for  this  function  is: 

y  =  8.7X“-^-' 

where,  y  =  mean  channel  width  (m),  and  X  =  drainage  basin  area  (km^) 

The  r^  value  is  0,34.  Tliis  result  is  statistically  significant  at  the  95%  level 
The  results  of  this  test  are  displayed  in  table  6.1  and  are  compared  with  the  distribution  found 
in  the  field  data,  No  underflow  type  jams  are  predicted  by  the  model.  Nolehoc  and  Worsham 
Creeks  are  predicted  to  have  predominantly  dam  type  jams.  Ten  creeks  are  predicted  as  having 
deflector-type  jams,  while  the  largest  five  creeks  arc  predicted  to  have  flow-parallel  jams.  This 
distribution  is  in  reasonable  agreement  with  the  distribution  found  in  the  filed  data,  the 
predominant  jam  type  being  predicted  correctly  for  1 1  of  the  17  reaches. 

The  differences  between  predicted  and  observed  jam  types  can,  for  the  most  part,  be  explained 
by  variations  in  local  channel  stability,  which  affect  channel  dimensions,  but  are  not  accounted 
for  in  the  simple  model,  Nolehoc  Creek,  for  example,  has  a  large  number  of  defiector  type  jams 
rather  than  the  expected  undertlo\/  jams  despite  having  the  smallest  drainage  area  and, 
therefore,  the  nanowest  predicted  channel  width  The  reason  for  this  anomaly  probably  lies  in 
the  fact  that  Nolehoe  Creek  is  highly  unstable  (stage  2  type  reach  in  Schumm  ct  al.’s  Channel 
Evolution  Model  (CEM))  and  is,  therefore,  deeply  incised.  Large  trees  close  to  the  channel 
edge  are  therefore  prone  to  collapsing  into  the  channel  to  surcharge,  whicfi  results  in  mass 
bank  failure.  Consequently,  debris  vertical  angle  (beta  angle)  is  greater  than  zero,  making  the 
span  of  fallen  trees  1  ss  than  they  would  otherwise  be  if  the  debris  lay  horizontally. 

In  Long  Creek  there  is  a  disproportionate  number  of  parallel  type  jams  because  much  of  the 
debris  found  in  the  reach  is  ‘paleodebris’.  that  is  material  leintroduced  into  the  channel  from 
erosion  of  old  alluvial  deposits  containing  debris,  which  has  broken  down  over  time  and  is, 
therefore,  found  in  great  abundance  in  individual  elements  of  short  length. 

Abiaca  Ch  eek,  site  3  is  predicted  to  have  predominantly  deflector-type  jams,  but  actually  has  a 
majority  of  dam  type  jams.  The  reason  for  this  discrepancy  is  because  channel  width  varies 
considerably  about  the  mean  within  the  reach,  with  wide  bend  sections  separated  by  much 
shorter,  nanower  straight  reaches  which  contain  all  the  dam-type  jams 
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Coila  Creek,  which  has  an  even  greater  average  width  than  Abiaca  (3),  contains  two  underflow 
type  jams.  These  have  fonned  because  the  channel  reach  has  re-stabilised  (CEM  stage  5)  and 
the  jams  are  located  in  a  short,  narrow  reach  which  has  been  narrowed  by  berm  formation. 

It  is  evident,  therefore,  tliat  for  the  model  to  yield  improved  predictions  of  jam  type,  it  must 
incorporate  further  predictive  variables.  A  channel  instability  factor  should  be  added  to  account 
for  the  deviation  of  channel  geometry  from  that  predicted  by  the  simple  hydraulic  geometry 
relationship  for  stable  conditions  Also,  an  additional  factor  is  necessary  to  indicate  the  rate  and 
distribution  of  lateral  channel  migration. 

Functions  could  also  be  included  to  predict  the  volume  of  debris  and  number  of  debris  jams 
likely  to  be  found  in  a  reach,  using  regression  relationships  derived  from  observations  of  debris 
volume  and  number  of  jams  per  unit  reach  length  as  a  function  of  drainage  basin  area. 

Channel  instability  and  lateral  migration  factors  could  also  be  used  to  predict  the  volume  of 
debris  likely  to  be  input  to  a  reach  by  estimating  the  rate  of  bank  retreat  The  input  rate  should 
be  reflected  in  the  storage  volume  and  frequency  of  jams.  An  attempt  is  made  in  section  6.3.2 
to  incorporate  reach  debris  volume  and  ..umber  of  jams  per  unit  reach  length  as  prediction 
factors 

The  model  also  attempts  to  predict,  in  a  simple  manner,  the  impact  that  debris  jams  have  in 
retaining  sediment  load,  causing  bed  scour  atid  generating  bank  erosion  This  estimale  is  based 
solely  on  the  median  sediment  grain-size  The  net  sediment  balance  volumes  (average  sediment 
stored  at  each  jam  minus  average  volume  of  sediment  scour  by  each  jam)  for  each  reach  (see 
table  5  2),  compared  with  predicted  relative  sedimentation  potential  in  Table  6.1  show  no 
significant  correlation.  Hence  predictive  potential  in  the  model  is  purely  descriptive,  making 
comparison  with  field  data  rather  tenuous.  It  may  be  concluded  that  reach  gruin-size  alone  is 
too  simple  a  measure  of  jam  sedimentation  potential.  A  measure  of  the  energy  available  to 
tran.sport  the  sediment  load,  such  a  stream  power,  m  combination  with  sediment  charactei  isiics 
is  likely  to  be  a  more  successful  predictor  limpirical  models  of  flow  dynamics  and  scour 
around  in-channel  obstructions  are  more  likely  to  give  at  least  order  of  magnitude  predictions 
for  scour  and  sedimentation  potential  at  debris  jams.  I’he  use  of  a  such  model  is  discussed  in 
section  6  3  3. 

Sedimentation  processes  are  also  affected  by  the  vertical  uigle  ol'  the  debris  elements.  If,  Ibi 
example,  debris  is  steeply  inclined  more  bedload  is  likely  to  pass  through  the  jam  than  if  the 
debris  lies  horizontally  on  the  channel  bed.  This  factor  could  be  accounted  foi  in  further 
developments  of  the  model  by  incorporating  a  width-depth  ratio  factor  which  predicts  ‘true’ 
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debris  length  to  channel  width  ratios  and  the  percentage  of  channel  cross-sectional  area 
blocked.  An  attempt  is  made  to  incorporate  a  beta  angle  factor  in  section  6,3,1 , 

Limitations  to  the  model  include  the  fact  that  jam-type  prediction  is  very  much  dependent  on 
the  width  function  selected  and,  therefore,  the  goodness  of  fit  of  the  hydraulic  geometry 
relationship.  To  make  the  model  more  robust  it  may  be  necessary  to  include  channel  width  as  a 
variable  input  by  the  user,  rather  than  as  a  factor  predicted  by  spatial  functions  based  on 
upstream  drainage  area  Alternatively,  it  may  be  more  realistic  to  make  jam-type  a  function  of 
CEM  stage,  because  this  has  strong  control  over  channel  geometry. 


Table  6.1  Model  Results 
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6.3  Revisions  to  the  Model 

The  results  of  the  tests  described  above  indicate  that  the  initial  theoretical  model  is  currently 
limited  as  a  piactical  predictive  tool,  but  has  the  potential  for  improvement  in  predicting  debris 
jam  type  and  associated  channel  geoinorphology  if  additional  factors  are  incorporated. 
Improvements  to  the  model,  outlined  below,  utilise  further  theoretical  concepts  and  empirical 
results  from  the  field  data  results, 

6.3.1  Debris  Beta  Angle 

It  was  noted  in  section  6,2  that  if  the  beta  angle  of  key  debris  was  much  greater  than  zero 
degrees  (debris  lying  horizontal  on  the  channel  bed)  the  jam  type  might  not  be  predicted 
correctly  (such  as  in  the  case  of  Nolehoe  ('reek).  Beta  angle  will  also  affect  the  extent  to  which 
floating  debris  is  trapped  by  the  jam,  the  flow  field  and,  therefore,  scour  and  sedimentation 
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processes.  Beta  angle  can  be  predicted  by  assuming  that  key  debris  elements  stand  close  to  the 
channel  bank  before  they  fall  in  and  then  topple  forward  into  the  channel,  forming  the 
hypotenuse  of  a  right-angle  tri,  ngle  with,  bank  height,  and  the  actual  width  spanned  by  the 
element  being  the  other  two  sides.  A  better  prediction  of  jam  classification  (channel  width  to 
tree  height  ratio)  can  then  be  determined  using: 

w  -  ylo'-  -  (6.1) 


where,  w  =  actual  bottom  width  spanned  by  debris  (m),  T  =  tree  height  (m)  and  B  =  bank 
height  (m) 

This  formula  can  be  directly  incorporated  into  the  model  to  improve  predictions  debris  jam- 
type  Correlation  between  beta  angle,  scour  and  sedimentation  volumes  for  each  reach  could  in 
theory  be  used  to  improve  prediction  of  jam  impact  on  sediment  processes  No  statistical 
correlation  between  beta  angle  and  either  the  retained  sediment  volume  or  the  scour  volume 
was  found,  however,  and  it  is  therefore  not  meaningful  to  add  any  ivlationship  between  these 
factors  to  the  model, 

6.3.2  Prediction  of  Debris  Volume  and  Number  of  Jams 

Prediction  of  the  average  number  of  debris  jams  and  average  volume  of  debris  in  a  reach  would 
be  beneficial  in  directing  LWD  management  activities  to  those  reaches  where  debris  build  up  is 
significant  The  conceptual  relationships  between  these  parameters  and  drainage  basin  area  and 
unit  stream  power  were  outlined  in  Chapter  5  Analysis  of  the  field  data  revealed  statistically 
insignificant,  correlations  for  all  these  relationships,  except  stream  power  and  debris  volume  It 
must  be  concluded  that  debris  volumes  and  number  of  jams  per  unit  reach  length  cannot  be 
accurately  predicted  using  simple  regression  formule. 
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and  the  occurrence  of  jams  because  these  factors  to  a  great  extent  determine  the  rate  of  debris 
input,  Examination  of  the  distribution  of  average  debris  volumes  for  each  stage  of  the  Channel 
Evolution  Model  is  discussed  in  chapter  5  and  summarised  in  Table  5,1, These  values  can  be 
included  in  the  model  to  improve  prediction  of  the  magnitude  of  debris  build-up  relative  to  that 
found  in  stable  reaches.  Reach  stability  has  been  added  to  the  model  as  a  further  input  variable. 
The  user  must,  therefore,  be  familiar  with  the  Cihannel  Evolution  i  'odel  in  order  to  make  an 
a.ssessment  of  the  reach  in  question  based  up(m  width/depth  ratios,  bank  stability  and  depth  of 
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sediment  in  the  channel  bed.  Stability  can  also  be  inferred  using  the  mode’  itself  on  the  basis  of 
the  average  width/depth  ratios  which  characterise  reach  type  (see  width/depth  ratios  in  table 
3.1).  This  classification  is  not  straight  forward,  however,  as  ranges  of  width/depth  ratio 
associated  with  particular  CEM  stages  overlap  and  the  model  must  be  given  a  llirther 
parameter,  such  as  depth  of  sediment  in  the  channel  or  average  bank  angle  to  allow  it  to 
discriminate  between  evolution  stages 
6.3.3  Prediction  of  Scour  ami  Sedimentation  Processes 

When  the  initial,  conceptual  model  was  developed  there  was  no  literature  base  that  dealt 
specifically  with  scour  and  sedimentation  processes  around  LWD.  The  simple  approach 
adopted  in  the  model  for  estimating  the  extent  of  scour  and  sedimentation  processes  proved  to 
be  rather  inadequate  because  of  the  lack  of  an  empirical  basis. 

Recent  research  by  Abbe  and  Montgomery  (1996)  partially  addresses  this  issue  by  attempting 
to  predict  channel  scour  associated  with  jams  based  upon  experimental  and  obscn/ational 
studies  of  scour  around  bridge  piers  and  abutments.  According  to  Raudkivi  (1990),  scour 
processes  can  be  divided  into: 

1 )  General  scour  which  occurs  due  to  increasing  discharge  or  slope,  irrespective  of  an  ; 

2)  Constriction  scour,  which  occurs  due  to  a  reduction  in  channel  cross-sectional  area, 

and, 

3)  Local  scour,  which  occurs  due  to  an  obstructions  direct  etiect  on  the  flow  field. 

Abbe  and  Montgomery  approximate  local  scour  around  debris  in  a  large  river  using  a  model 
developed  by  Liu  et  al.  ( 1961 )  which  predicts  clear  water  scour  around  an  abutment: 
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where,  di,  =  depth  of  local  scour  (m),  h  -  flow  depth  (m).  La  =  abutment  length  (m),  Fr  =  IJ  / 
(gh)'”,  and  U  =  mean  flow  velocity  (m/s'‘). 

Constriction  scour  is  approximated  using  an  empirical  model  presented  by  Laursen  (1963): 
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wh  di:s  =  depth  of  constriction  scour  (m),  .  =  bed  shear  stress  (N/m^),  =  critical  bed 

shear  stress  (N/m^),  wbl  =  unobstructed  flow  width  (m),  and  wb2  =  constricted  flow 
width  (m). 

Total  scour  depth  due  to  an  obstruction  is  calculated  by  Abbe  and  Montgomery'  as  the  sum  of 
scour  values  from  equations  6,2  and  6.3. 

It  was  found  that  the  values  predicted  by  these  equations  significantly  overestimated  scour 
depths  in  the  study  reach  (Queets  River,  Washington  State)  of  an  annoured,  gravel-bed  river. 
Experimental  work  by  Raudkivi  and  Ettema  (1973)  showed  that  scour  depths  were  greatly 
influenced  by  sediment  sorting  and  Raudkivi  (1990)  presents  a  graphical  method  for'  adjusting 
scour  prediction  that  predicts  adjusted  scour  depths  as  a  function  of  the  standard  deviation  of 
the  reach-averaged  sediment  grain  size.  When  this  adjustment  factor  was  applied  to  Queets 
River  data,  the  total  estimated  scour  was  found  to  be  within  3%  to  17%  of  pool  depths 
obser\'ed  in  the  study  reach.  These  two  equations  were  incorporated  into  an  updated  version  of 
the  model  because  they  appear  to  provide  reasonable  approximations  for  prerliction  of  scour 
depths  at  debris  jams.  Given  the  type  of  debris  accumulations  studied  in  the  Queets  River, 
equations  6.2  and  6.3  are  probably  best  suited  to  describing  scour  conditions  around  flow 
parallel  and  deflector-type  jams,  where  the  channel  width  is  considerably  greater  than  the 
obstruction  width.  Follov/ing  the  successful  use  of  these  models  by  Abbe  and  Montgomeiy, 
other  scour  models  were  explored  to  attempt  to  approximate  scour  conditions  at  dam  and 
underflow  type  jams. 

If  the  jam  type  encountered  in  a  reacli  is  of  the  ‘dam’  type  and  fully  blocks  the  channel  width, 
constriction  and  local  scour  are  unlikely  to  be  significant.  Instead,  scout  may  occur 


downstream  of  the  jam  due  to  a  jet  of  water  plunging  over  the  dam  ‘weir’.  Raudkivi  (1990) 
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and  similar  hydraulic  structures.  The  equation  developed  by  Mason  and  Aruumugam  (1985) 
(see  Raudkivi,  1 990)  should  be  appropriate  for  approximating  plunge  scour  at  a  dam-type  jam, 
however,  because  it  satisfies  Froude-  law  scaling.  The  model  is  expressed  as 
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where,  D  =  flow  plus  scour  depth  (m),  q  =  unit  d'scharye  in  the  jet  (m/s‘Vm),  H  =  head 
difference  between  reservoir  crest  and  tailwater  surface  (tn),  h  =  tailwater  depth  (m),  g 
gravitational  constant  (9.81  and  d^  =  mean  grain  size  (m). 

Scour  at  underflow  jams  is  rather  more  difficult  to  predict  as  the  jam  will  only  impinge  on  the 
flow  above  a  certain  flow  stage,  depending  upon  the  height  of  the  debris  element  above  the  bed 
(this  need  not  necessarily  be  bank  top  height).  A  formula  is  proposed  by  Shalash  (1959),  (in 
Raudkivi  (1990)),  to  determine  the  scour  caused  by  an  underflow  jet,  but  this  model  assumes 
that  there  is  a  solid  horizontal  apron  on  the  channel  bed  which  extends  dowrtstream  of  the 
obstruction  by  a  minimum  of  1,5  times  the  approach  flow  depth,  The  model  is,  therefore,  not 
applicable  to  channels  with  erodibie  beds.  Scour  at  underflow  jams  cannot  be  adequately 
represented  by  empirical  formula  that  have  been  developed  to  date,  and  it  remains  for  future 
studies  to  model  this  phcjnornenon. 

Estimation  of  tlie  volume  of  sediment  stored  by  a  debris  jam  is  a  more  complex  problem,  and 
no  research  has  been  found  which  attempts  to  model  local  distribution  and  volumes  of 
sediment  storage  in  alluvial  channels.  Sediment  is  stored  at  two  distinct  locations  around  debris 
jams.  The  fiist  storage  zone  is  in  the  zone  of  low  shear  stress  downstream  of  the  obstruction 
where  streamlines  are  separated.  The  second  storage  zone  is  upstream  of  jams,  but  is  only 
significant  if  the  jam  blocks  a  sufficient  percentage  of  the  channel  cross-section  to  cause  a 
backwater  effect.  Sedimentation  in  this  backwater  zone  is  similar  to  that  in  a  reservoir. 
Reservoir  sedimentation  has  been  found  to  be  dependent  upon  sediment  size  and  grading,  size 
and  shape  of  the  reservoir,  inflow  and  outflow  rates,  and  type  and  location  of  the  outflow.  The 
formation  of  deltas  in  resciweirs,  let  alone  in  local  channel  reaches  is  disTicaU  to  model  and  has 
mostly  been  dealt  within  a  qualitative  manner,  although  a  few  analytical  models  have  been 
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therefore,  a  complex  problem  and  it  is  beyond  the  scope  of  this  thesis  to  develop  analytical 
models  to  adequately  describe  it, 

6.3.4  Rerining  the  iVlodel 

The  foilowing  revisions  were  made  tr?  the  tiurdel  in  the  light  of  the  preceding  analysis  and 
discussion, 

•  V/idth  Input 

Average  channel  width  may  now  be  entered  into  the  mode!  using  one  of  three  methods.  The 
first  option  is  to  use  the  hydraulic  geometry  function  derived  by  Scbumni  et  al.  (1984)  for 


channels  in  northern  Mississippi  The  second  option  is  to  enter  a  user-defined  hydraulic 
geometry'  function  of  a  simple  power  fonn  (Y  =  ax'").  The  third  option  is  to  enter  an  observed 
value  for  the  width  directly. 

«  Debris  Length 

Equation  6.  i  is  used  in  the  modified  model  to  determine  the  effective  debris  length  in  the 
chamiel,  if  the  jam  type  is  determined  to  be  other  than  the  undeiflow  type.  This  calculation 
requires  a  cliantiel  depth  value  which,  like  the  channel  width,  may  be  entered  into  the  model 
using  any  one  of  the  three  options  described  above.  The  hydraulic  geonictiy  Ibnction  delined 
by  Schumm  et  al  (1984)  for  channel  depth  in  quasi-equilibrium  reaches  in  northern  Mississippi 
is: 


Average  depth  (feet)  -  12.27  drainage  area**  ^' 


(6.5) 


<•  Debris  Volume  and  Jam  Frequency 

Reach  stability,  defined  by  the  Channel  Evolution  Model  (Schumm  et  al.,  1984),  is  used  to 
predict  the  approximate  volume  of  debris  and  the  approximate  number  of  debris  jams  to  be 
found  per  100m  ach  (see  table  5.1).  The  channel  evolution  stage  may  be  entered  directly  by 
the  user,  or  the  model  can  predict  reach  stage  on  the  basis  of  the  calculated  width/deptb  ratio 
and  the  mean  depth  of  sediment  on  the  channel  bed.  The  width/depth  ratio  and  sediment  depth 
ranges  for  each  stage,  which  are  incorporated  into  the  model,  are  taken  fiom  the  field  data 
analysis  presented  in  table  3,1.  The  sediment  thickness  variable  is  entered  by  the  user  and  must 
be  based  on  field  measurement  or  estimation. 

•  Scour  Prediction;  deflector  and  flQW"pHr,?HeS  jams 

For  deflector  and  flow  parallel  jams  maximum  local  and  constriction  scour  is  piedicted  using 
equations  6  2  and  6.3.  Total  scour  is  calculated  as  the  sum  of  these  figures.  Mean  flow 
velocity,  required  to  calculate  the  Froude  number  in  equation  6.2,  is  determined  using  the  Maw 
of  the  wall’  equation  for  turbulent  flow  (see  Chow,  1959,  pg.  201); 
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where,  v  ~  time  averaged  velocity  at  elevation  d  above  the  bed  (m/s  ’),  v*  =  shear  velocity  = 
(gds)“’  (m/s’’),  g  =  gravitational  constant  (9.81  m/s'^),  s  =  bed  slope,  k  =  von  Karman 
constant  =  0.4,  and  do  =  boundary  roughness  height  =  30/06?  (appi  ox.)  (m). 

The  depth  value,  d,  is  multiplied  by  0  4  in  equation  6.6  to  obtain  v,  the  mean  depth-  averaged 
velocity.  Do?  is  substituted  by  D50  in  the  model  to  reduce  the  data  input  required  by  the  user. 
La  is  taken  as  the  effective  debris  length  in  the  model  if  the  jam  type  is  dam  or  deflector  and  as 
average  root  wad  diameter  if  the  jam  type  is  flow  parallel.  Scour  is  not  calculated  if  the  jam 
type  is  underflow,  because  the  model  is  only  applicable  to  obstructions  resting  on  the  channel 
bed, 

In  equation  6.3,  the  critical  bed  shear  stress  (Tc)  is  determined  using  the  following  formula, 
suggested  by  Pitlick  (1992): 


Tc  =  0(p«-pw)gD»  (6.7) 

where,  6  =  critical  Shields  Parameter  value  (approx.  0,05),  p,  =  density  of  sediment,  taken  to 
be  approx.  2650  kg/nv\  and  Pw  =  density  of  water  -•  1000  kg/In^ 

Bed  shear  stress  (!„)  for  the  flow  depth  being  modelled  is  calculated  as: 


3;,  pw  g  R  S 


(6.8) 


where,  R  hydraulic  radius  (A/P)  (m),  A  -  mean  channel  cross-sectional  area  (w  x  d)  (m), 
and  P  -  wetted  perimeter  =  (2  x  d  +  w)  (m) 

The  obstructed  flow  witiili  is  calculated  as  channel  width  minus  effective  debris  length  rf  the 
jam  type  is  deflector,  and  r's  chairne!  width  minus  root  wad  diameter  if  the  jam  type  is  flow 
parallel.  If  the  jam  type  is  underflow  constr  iction  scour  is  not  calculated  because  the  model  was 
developed  only  tor  obstructions  resting  on  the  channel  bed. 

it  is  left  to  the  user  to  adjust  pool  scour  atrd  constriction  scour  values  using  the  method 
suggested  by  Raudkivi  (1990,  pg.  251)  if  bed  armouring  processes  are  thought  to  be 
significant. 

•  Scour  Prediction;  Darn-type  jani,s 

ir  dam  type  jams,  scour  depth  is  calculated  using  equation  6  4.  Unit  discharge  is  calculated  as 
V  X  flow  depth.  H,  the  head  difference  across  the  dam,  is  takf n  be  approximately  equal  to 


the  jam  height  minus  tailwater  depth,  Jam  height  is  simply  taken  as  diameter  of  the  key  debris 
element.  The  modified  model,  therefore,  requires  the  input  of  mean  tree  trunk  diameter.  The 
tailwater  depth  (h)  is  very  difficult  to  estimate  and  is  set  somewhat  arbitrarily  in  the  model  to 
be  20%  of  the  approach  flow  depth,  giving  ‘worst  case’  scour  conditions.  The  mean  sediment 
diameter  (d„,)  is  replaced  by  d5()  in  the  model.  Depth  of  scour  is  D  (total  flow  and  scour  depth) 
minus  tailwater  depth. 

•  Scour  Prediction:  Underflow-type  Jam 

For  underflow  type  jams  scour  is  not  calculated  owing  to  the  difficulty  in  applying  any  of  the 
equations  presented  in  the  literature. 
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7  HYORAUMC  SIGNIFICANCE  OF  I,WD 


A  comprehensive  investigation  of  the  hydraulic  cftect  of  LWD  in  rivers  has  not  to  date  been 
documented  However,  some  studies  have  investigated  the  eftect  of  LWD  on  channel  roughness, 
runoff  hydrographs,  velocity  distributions  and  the  water  surface  profile, 

7.1  Effect  of  LWD  on  channel  roughness 

The  Manning's  equation  generates  a  resistance  coefficient  that  represents  all  sources  of  roughness  in 
the  channel  This  equation  is  widely  used  by  river  engineers  who  select  values  of  "n",  ftom  personal 
experience,  tables  in  Chow  (1959)  or  photographs  in  Barnes  (1967).  The  range  of  “n”  factor  in 
normal  channels  is  between  0.025  and  0. 15,  For  heavily  congested  streams  less  than  30m  wide  n 
ranges  Irom  0  075  to  0, 1 5.  Irregular  and  rough  reaches  of  large  streams  have  values  of  ‘  n"  from 
0,035  to  0,10, 
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(7,1) 


where,  n  ^  Mannings  rouglincss  coefficietit,  R  hydraulic  radius  (m),  S  ^  energy  slope,  V  -  mean 
velocity  (ms  '),  1 .49  ^  conversion  to  fps  units. 

Rather  than  Using  Manning’s  “ii”,  the  more  theoretically  based  Darcy-Weisbach  flow  resistance 
equation  was  used  (Richards,  1982)  which  is  expressed  as: 


J/2 


(7.2) 


where,  f  Darcy-Weisbach  friction  factor,  R  hydraulic  radius  (m),  <md  Sw  energy  skipe 
The  effect  of  LWD  on  flow  resistance  varies  as  a  tunction  of  relative  flow  depth  of  Bevari  et  al 
(1979)  found  that  when  LWD  is  large  in  relation  to  flow  depth  the  roughness  ct^efficient  is 
extremely  high  (  Manning's  n  >1 ).  As  I,,W13  becomes  increasingly  submerged  it  exerts  less  intluence 
on  flow  resistance.  Shields  and  Smith  (1992)  measured  a  large  decrease  in  Darcy-Weisbach  friction 
factor  as  discharge  increased,  and  also  observed  that  friction  tactors,  for  cleared  and  uncleared 
reaches  converged  at  high  in-bank  flows.  Indirect  evidence  to  support  lliese  findings  is  provided  by 
investigations  of  downstream  hydraulic  geometry  which  show  that  rougliness  gener  ally  decr  eases  a.s 
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channel  size  increases  (Wolman,  1 955).  Pctryk  and  Bosniajian  (1975)  derived  an  equation  to  predict 
Manning's  “n”  as  a  function  of  density  of  vegetation  in  the  channel,  hydraulic  radius,  Manning’s  “n” 
due  to  boundaiy  roughness  and  a  vegetation  drag  coefficient. 
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where,  nj,  =  Manning's  boundary  roughness  coefficient  excluding  the  etfect  of  vegetation,  Cd  =  drag 
coefficient  for  vegetation  (assumed  to  be  1),  Ai  =  projected  area  of  the  ith  plant  in  the 
streamwise  direction  (m^),  A  =  cross-sectioucil  aiea  of  flow  (m^),  L  length  of  the  channel 
reach  being  considered  (m),  and  P  =  (m), 

In  this  fonnula  the  expression  CdLAi/AL  represents  the  density  of  vegetation  in  the  channel.  Gippel 
et  al.  (1992)  note  that  a  problem  with  this  fonnula  lies  in  selecting  a  value  for  the  drag  coefficient, 
Cd.  Petryk  and  Bosmajian  assumed  a  value  of  1,  but  this  strictly  only  applies  to  cylinders  in  infinite 
flow  In  streams,  interfei'^nce  Ifom  nearby  obstnictions  and  the  efl'ect  of  blockage  on  the  drag 
coefficient  must  also  be  considered. 

Manning’s  equation  is,  liowever,  inapplicable  in  situations  where  there  is  a  high  degree  of 
obstfuetion  in  the  chaiinel,  particularly  where  n  >1  It  was  developed  empirically  to  describe  open 
channel  situations  with  fi'lly  turbulent  flow,  whcie  firiction  is  controlled  primarily  by  skin  friction  at 
the  channel  boundary.  The  equation  attaches  significance:  to  t’  e  hydraulic  radius  which  may  be 
iiTelevant  if  the  channel  is  heavily  choked  with  1  .WD 

Smith  and  Shields  (1992)  studied  the  effects  of  varying  levels  of  1-WD  density  on  the  physical 
aquatic  habitat  of  South  Fork  Obion  River,  Tennessee,  USA  Two  secondary  objectives  in  this 
study  were  to  develop  a  id  demonstiate  a  method  for  quantifying  LWD  in  a  given  reach,  and  to 
relate  the  quantity  of  TVD  to  reach  hydraulics  An  approach  similar  to  that  u.sed  by  Petiyk  arid 
Bosmajian  (1975)  was  uicd  to  calculate  the  effect  of  LWD  on  channel  roughness.  The  LWD 
density  in  a  reach  was  calt  ilated  using  the  following  formula: 

nA  -  I  '^'7,  .  (i  /  (7.4) 

,  I  ^  I  k  \ 
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where,  DA  =  LWd  density  (m^),  n  -■  total  number  of  LWD  formations  in  the  reach,  Ai  =  area  of  the 
ith  debris  formation  in  the  plane  perpendicular  to  How  (m^),  A  =  reach  mean  flow  cross- 
sectional  area  (m^),  U  =  reach  length  (m),  =  formation  type  weighting  factor  for  jth 

foimation  type,  Nj  k  =  numbei'  of  type  j  LWD  formations  in  Kth  width  icategory,  and  = 
weighting  factor  based  on  t^WD  formation  width  category. 

Appendix  B  presents  the  Smith  &  Shields  (1992)  survey  form  and  shows  how  to  calculate  the 
weighting  tactois, 

In  a  channel  reach  where  LWD  plays  a  major  role  in  flow  resistance,  total  resistance  can  be 
expressed  as: 


fl  -  fb  +  fi 


(7,5) 


where,  tj  ==  total  Darcy- Weisbach  friction  factor,  Ij,  =  boundary  friction  factor  excluding  I.WD 
effects,  and  fa  =  friction  factoi  due  to  LWD. 

Total  head  loss  is  the  sum  of  a  boundary  uiction  loss  and  a  LWD  blockage  loss,  as  follows: 


(7,6) 


where,  hj  -  total  head  loss  (m),  Si;  =  slope  of  the  energy  gradient,  and  Kj  =  dimensionless  loss 
coeflicient  (dependc.nt  upon  LWD  density). 

The  energy  gradient  (Si.)  can  be  calculated  using  a  total  friction  factor  from  the  Darcy-Weisbach 
equation: 


(8k«) 


Substituting  this  expression  foi  Si;  into  equation  6  gives: 
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The  ratio  Kj/L  may  be  expressed  in  terms  of  the  LWD  density  as; 


Ka/L  =  DA  (7.10) 

Smith  and  Shields  calculated  values  for  lb  using  curves  developed  by  Alam  and  Kennedy  (1969), 
hydraulic  parametas  determined  from  dye  tracer  tests  in  the  LWD  reaches  (which  provide  direct 
discharge  and  velocity  estimates  (Richards  1 982)),  and  the  median  bed  grain  size  cletermined  from 
sieve  analysis.  Values  for  fj  were  then  calculated  using  equaiions  7  3,  7  9  and  7  10.  They  then 
compared  computed  values  of  with  values  measured  using  dye  tests 

The  results  of  their  study  showed  a  reasonable  positive  correlation  between  the  measured  and 
computed  friction  factors  However,  they  recognised  that  considerable  refinement  and  site-specific 
adaptation  may  be  required,  and  that  the  method  does  not  account  for  local  energy  loss  because  of 
bends  or  flow  expansion  and  contraction  at  bridges,  debris  dams  or  riffles  The  method  docs  have  a 
sound  theoretical  basis,  however,  and  could  be  uselully  employed  in  future  research  into  the 
hydraulics  effect  of  LWD 


7.2  Effect  of  LWD  on  velocity  distribution 

LWD  clearly  influences  the  direction  and  magnitude  of  flows  currents  within  stream  flow,  but  few 
data  have  been  diKumented  in  the  literature.  Swanson  (1979)  produced  detailed  maps  of  debris 
jams  indicating  flow  orientation  with  directional  arrows.  Smith  and  Shields  (1990)  reported  that  the 
removal  of  LWD  from  a  river  18-23m  wide  3  5  to  4  5  m  deep  produced  more  unitbrni  flow,  with 
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flow  field  and  velocity  distribution  are  not  well  established 


7.3  Effect  of  LWD  on  stage/discharge  relationships,  the  hydrograph  and  flood  frequency 

LWD  i,,  often  removed  because  it  is  believed  that  this  will  achieve  a  significant  reduction  in  channel 
roughness  which  will  produce  a  higher  mean  velocity  and  thereby  increase  in-bank  channel  flow 
capacity.  There  is  some  evidence  to  suppoit  this  assumption.  For  example  Smith  and  Shields 
(1990)  measured  the  mean  velocity  in  two  cleared  reaches  of  a  river  to  be  0.04  m/s  and  0.34  m/s 
respectively.  In  an  uncleared  reach  of  the  same  river  the  mean  velocity  was  0  27  m/s.  MacDonald 
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and  Keller  (1987)  also  found  that  there  was  a  local  increase  in  velocity  by  up  to  250%  as  a  result  of 
LWD  removal  and  a  decreased  sinuosity  of  the  low-flow  thalweg.  According  to  Gippel  et  al  (1992) 
the  Murray-Darling  Basin  Commission  calculated  a  theoretical  reduction  in  water  level  of  0. 3-0.4  m 
alter  the  removal  of  approximately  200  snags  per  kilometre.  However,  later  analysis  of  flow  records 
indicated  a  reduction  of  only  0,2  m.  In  theoiy,  there  should  be  a  statistical  reduction  in  the 
magnitude  and  frequency  of  overbank  flooding  where  debris  is  removed  from  a  channel  because  of 
the  increased  channel  capacity.  Bodron  (1994),  used  a  dynamic  routing  model  to  demonstrate 
changes  in  both  stage  and  duratiott  of  flood  events  before  and  after  LWD  removal,  using  Manning 
“n”  values  calculated  in  the  study  by  Smith  and  Shields  at  South  Foi  k  Obion  River,  west  Tennessee 
Despite  the  fact  that  the  increase  in  channel  cross-sectional  area  due  to  LWD  removal  was  ignored, 
small  reductions  in  flood  height  and  duration  were  calculated  based  solely  on  the  change  in 
Manning’s  “n”.  Bodron  also  noted  that  flood  stage  would  be  reduced  further  if  sediment 
accumulations  at  each  jam  site  had  been  removed  However,  according  to  Gippel  et  al  (1992)  many 
claims  that  this  cftect  has  been  achieved  lack  any  hard  supporting  evidence,  Counterclaims  also  lack 
supportive  evidence,  because  of  the  difficulty  of  isolating  the  hydraulic  effect  of  LWD  removal  It  is 
even  possible  that  LWD  removal  might  increase  downstream  flood  peaks,  because  in  the  smoother 
channel  the  flood  wave  is  less  attenuated, 

Gregoiy  et,  al,  (1985)  found  that  LWD  ponds  water  which  results  in  an  increase  in  water  depth  aiid 
a  decrease  in  velocity  wliich,  at  low  flows,  influences  travel  time  significantly.  At  high  flows, 
however,  the  ponding  effect  of  l.WD  is  drowned  out.  Shields  rmd  Nunnally  (1984)  noted  that 
because  large  accumulations  of  I.WD  have  a  damming  e!fa4  on  the  flow  which  locally  elevates  the 
base  level  they  can  be  treated  as  geometric  elements  within  the  channel,  rather  than  simply  as 
roughness  elements,  in  backwater  profile  computations. 

7.4  IVlodelling  the  hydraulic  effect  of  LWD 

Most  studies  of  resistance  to  flow  in  livers  liave  concentrated  on  small-  scale  roughness,  especially 
skin  friction  offered  by  bed  sediments,  where  the  size  of  the  roi  ghness  element  is  small  compared  to 
the  flow  depth  L,WD,  on  the  other  hand,  represents  a  form  of  large-scale  r  oughness  for  which  skin 
fnetion  is  small  compared  to  fonn  drag  (Petryk  and  Bosmajian,  1975).  Flow  conditions  associated 
with  the  presence  of  l.WD  in  streams  varies  from  sub-critical  to  super-critical  depending  on  the 
dimensions  of  the  LWD  and  the  depth  of  water. 

Gippel  et.  al.  (1992)  used  the  momentum  principle  to  detemiine  the  hydraulic  effect  of  LWD,  the 
effect  being  quantified  in  terms  of  afflux  or  backwater  effect.  If  flow  is  subcritical  (Froude  number  ' 
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1),  then  apart  from  local  distui'bancc  of  the  velocity  profile,  LWD  only  has  an  influence  in  the 
upstream  direction  Quantifying  the  backwater  effects  is  problematical  because  of  the  practical 
difficulties  of  directly  measuring  the  afflux  at  debris  jams.  An  alternative  to  direct  measurement  is 
prediction  on  the  basis  of  a  knowm  relationship  between  afflux  and  more  easily  measured 
parameters,  Gippcl  et  al  used  the  results  of  a  laboratory  hydraulic  study  to  develop  a  method  of 
detemiining  the  afflux  caused  by  LWD  (see  Fig'are  7, 1 ). 

They  propose  the  use  of  the  following  equation  to  calculate  afflux  : 


A/i  - 
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(7.11) 
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where.  Ah  =  afflux  ==  hi  -  h,i  (m),  and  the  drag  coefficient  (Cp)  is  given  by: 


plI^Ld 


(7,12) 


where,  I  jj  drag  force  (N),  p  -  density  of  water  (lO(X)  kg/nv'),  Ui  -  mean  velocity  at 
section  upstream  of  object  (tn/s’’),  L-  -  projected  length  of  l,WD  in  flow  (m),  d  -  diameter 
of  LWD  (m) 

and  the  Froude  number  (F)  is  defined  by: 
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where,  mean  velocily  at  section  downstream  of  object  (m/s'),  hi  -  water 
depth  downstream  of  l,WD  (m) 
and  the  blockage  i  atio  (B)  is  given  by: 


B  -  l..d  /  A 

where,  A  W.h|  cross  sectional  area  ol’llow  (m'’) 
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(7.14) 


Figure  7.1  Definition  sketch  of  LWD  model  by  Cippel  et  al.  (1992) 
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Thus,  the  afflux  depends  on  !<,  Cd  and  B  The  hroude  number  can  be  calculated  from  direct 
measurement  or  from  flow  records.  B  can  be  found  from  a  field  survey  'I'he  remaining  problem 
centres  on  selecting  an  appropriate  drag  coefficient  The  drag  characteristics  of  a  cylindei’  in  infinite 
flow  are  well  known  (Potiyk  and  Bosmajian,  1975)  l  .ess  is  known  about  drag  on  cylinders  within 
boundaries  where  the  “blockage  cflect”  is  significant  and  the  drag  coefficient  is  consequently 
increased. 

Gippel  et  al,  conducted  experiments  on  l.WD  models  to  detcrmiiK  drag  force,  using  a  towing 
carriage  and  water  tunnel,  Froi  number,  1,WI )  length  to  diameter  ratio  and  LWD  depth  from  the 
bed  all  affected  drag  coeffl.cient,  but  were  much  less  important  than  the  blockage  effect,  angle  of 
orientation  to  the  flow  and  the  shielding  effect  (of  one  piece  of  I.WD  behind  another),  A  suitable 
drag  coefficient  (C’d)  for  the  LWD  in  question  can  be  selected  from  their  experimental  results 
((  ippel  et,  al  1992,  figures  1.8  or  1,12)  on  the  basis  of  its  overall  shape  and  angle  of  orientation, 
1  he  drag  coefficient  should  then  be  adjusted  for  the  blockage  elfect,  which  can  be  calculated  using 
the  following  equation  developed  by  Gippel  et.  al,  using  their  empirical  data  from  flume  studies: 

C,v-<'i)(I-B)-'  (7.15) 

where,  ('u  adjusted  drag  coefficient,  and  G'n  drag  coefficient  in  infinite  flow. 


These  data  are  then  substituted  into  equation  7. 1 1  to  calculate  the  afflux. 

Ptedicted  and  measured  afflux  values  resulting  from  the  flume  study  were  very  closely  correlated, 
and  they  conclude  that  the  flume  conditions  did  not  seriously  violate  any  of  the  assumptions  in 
equation  7,11.  The  proposed  method  of  afflux  estimation  was  then  applied  to  data  collected  from 
the  Thomson  River,  Victoria  and  revealed  that  de-snagging  there  would  produce  a  reduction  in 
stage  of  only  0  01m  at  bankfull  flow.  Appendix  C  presents  the  method  for  predicting  flow  afflux  due 
to  L,WD,  developed  by  Gippel  et  al,  (1092), 

In  conclusion  this  method  of  backwater  or  afflux  calculation  due  to  individual  items  of  LWD  could 
be  used  as  a  tool  to  help  determine  whether  the  afflux  reduction  due  to  LWD  removal  would  have  a 
significant,  positive  impact  according  to  the  perceived  management  requirements  or  whether  LWD 
could  be  left  in  place  perhaps,  re-orientated,  lopped  oi  even  re-introduced  where  sympathetic 
rehabilitation  management  is  desirable,  without  significant  effect  on  high  in-bank  stages. 
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8  IMPACT  OF  LWD  AT  BRIDGES 


ij  S.l  Introduction 

In  an  investigation  of  bridge  scour  research  needs,  Jones  et.  al.  (1991)  cited  the  affect  of 
debris  on  pier  scour  depths  as  a  subject  of  pressing  concern  that  required  model  studies  and 

field  observations  to  characterise  debris  build-up  Doheny  (1993)  observed  scour  conditions  at 

Ij 

876  l.ighway  bridges  in  Maryland  and,  amongst  other  relationships,  found  that  bridges  in 
forested,  urban  and  pasture  basins  were  more  prone  to  bloekage  than  those  in  basins  with  row 
crops  or  swamp,  Diehl  &  Biyan  (1993)  assessed  potential  debris  volumes  that  could  be 
')  transported  to  bridge  sites  in  the  West  Harpeth  River  basin,  Tennessee  and  found  bank 

instability  to  be  the  channel  characteristic  most  useful  in  identifying  channel  reaches  with  high 
potential  for  production  of  LWD.  Instability  through  channel  migration  and  mass  failure  or 
il  fluvial  erosion  can  be  detected  on  maps  and  aerial  photographs  (Diehl  &  Bryan,  1993).  A  study 

by  Parola,  Fenske  &  Hagerty  was  initiated  to  investigated  the  basin-wide  impact  of  the  1 993 
Mississippi  River  Basin  flooding  on  damage  to  the  highway  infrastructure.  Structural  geometry 
information,  as  well  as  hydraulic  information,  was  collected  at  two  sites  where  bridges  had 
'  collapsed  at  least  partly  as  a  result  of  debris  loading  and  was  noted  to  be  a  contributing  factor 

in  the  lateral  load  and  scour  of  many  bridges.  Plate.  7  shows  the  Missouri  1 1 3  bridge  over 
Florida  Creek,  where  floating  debris  was  a  key  factor  in  its  collapse. 

I 

Plate  7  :  Bridge  113  over  Florida  Creek,  .Skidmore  Missouri.  Failure  due  to  debris 
loading.  Source  ;  Parola,  Fenske  &  Hagerty  (1994). 
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7.2  Methsds  for  Managing  Floating  Debris  at  Bridges 

Only  one  paper  has  been  found  ti.at  directly  addresses  debris  management  at  bridges.  Saunders 
&  Oppenheirner  (1993)  believe  that  conventional  methods  of  protecting  piers  from  floating 
debris  are  inadequate.  The''  comment  that  the  use  of  pilings,  or  some  other  barrier  upstream  of 
a  bridge,  can  actually  e:...-eiDate  the  problem  because  the  debris  accumulated  may  all  be 
releasea  at  one  time  in  the  form  of  a  raft  which  cannot  be  pass  under  the  bridge.  They  describe 
a  novel  deflector,  a  iurate  shaped  hydrofoil  which  generates  counter-rotating  streamwise 
vortices  iii  its  wake,  positioned  below  the  surface  so  that  it  is  not  impacted  by  debris  upstream 
of  the  piers  and  so  ihat  the  '’ortices  migrate  to  the  surface  ahead  of  the  pier. 

The  ptin''ip!c  is  that  the  near  surface  flow  induced  by  the  vortices  deflects  debris  safely  around 
the  pier.  Figure  8.1  shows  the  hydrofoil  in  elevation  and  planform.  The  foil  is  mounted  on  a 
tether  or  pylon  at  a  depth,  d,  below  ihe  surtace  and  a  distance,  Zo,  upstream  cf  the  pier  and  is 
inclined  at  an  angle  such  that  the  three  on  the  foil  is  downwards  and  the  reaction  on  the  water 
causes  a  local  motion  upwards  lowards  the  surface.  After  intt^racting  with  the  vortex,  debris  is 
deflected  at  the  angle,  a,  and  is  displaceo  sideways  by  a  distance,  D,  by  the  time  it  reaches  the 
pier. 

A  flimu  inoeiel  constructed  by  Saunders  &  Oppenheirner  indicated  that  the  vorticity  remains 
highly  concentrated  tl'r  a  distance  of  about  20  tinns  the  span  of  the  hydrofoil,  b,  when  b  = 
O.bxh  (depth  of  flow).  The  problem  is  characterised  by  a  bridge  pier  width  w  and  by  the  size  of 
debris  An  .rvesge  debris  size  is  utilised  with  diameter  Dd  and  length  f„  The  vortex 
produced  hy  'he  devite  has  a  characteristic  diameter,  Dv,  of  order  b  (hydrofoil  span).  If  Dd  > 
Dv  then  the  v  utex  will  noi  impart  a  net  motion  to  the  debris,  so  they  recommend  a  value  o+'b 
>  2Dd  or  b  -  'v  (pier  vvidtli)  as,  they  assume,  the  majority  of  debris  will  have  a  diameter  less 
than  ilir  pw  ’vidtli  and  this  scaling  will  ensure  that  the  vortex  is  positioned  correctly  with 
respect  to  tlie  pier.  U  is  also  suggested  that  the  device  be  tetliciod  so  that  it  can  oscillate 
trunsveisei.  ♦ :  the  flow,  with  tiic  result  that  the  vortiec;  will  tend  to  destabilise  any  debris  tha^. 
.night  t'-tve  acc-urnulated  irr  djc  free  of  the  pier. 

In  fii.r.ie  tests  the  hydrofor  ts  Tep>'irted  to  work  very  otTecrively  and  the  device  would  appear  to 
ofl’ei  a  possible  approach  to  rnit’iaging  floating  debris  at  bridges  However,  if  the  average 
dcbiis  Icttgoi  is  greater  than  the  pier  spacing  debris  floating  with  their  lo.ng  axis  transserse  to 
the  How  arc  still  ..xely  to  be  trapped  and  the  vortices  might  even  turn  flow  parallel  debris 
(luoi!;.,li  '>0  degrees  so  it  is  more  lif.ely  to  become  iammeJ  betuzecn  adjacent  piers. 


8,3  DEBRIS  AT  BRIDGE  PIER  PREDICTION  PROGRAM  (DBP3) 


A  computer  model  has  been  developed  to  calculate  the  probability  of  debris  build-up  at  bridge 
piers,  and  the  associated  debris  induced  scour,  based  upon  modified  theoretical  equations 
published  by  Melville  and  Dongol  (1992)  and  Simons  and  Li  (1979)  Version  2.0  of  the  model 
(D6P3)  is  also  included  on  the  enclosed  disk. 

There  are  only  a  limited  number  of  studies  that  have  addressed  the  problem  of  debiis 
accumulations  at  bridges.  Melville  &  Dongol  ( 1 992)  look  at  the  problem  of  pier  scour  due  to 
debris,  while  Simons  &  Li  (see  Callander,  1980)  have  used  a  probabilistic  approach  to  quantify 
the  rate  of  bridge  span  blockage  by  debris  and  the  subsequent  backwater  effect,  pressure  forces 
generated  on  the  piers 

Local  scour  at  bridge  piers  has  been  extensively  investigated.  However  the  impact  of  debris 
rafts  at  piers  ryhich  create  additional  ffow  obstruction  and  therefore  increase  scour  depths  has 
beer,  largely  neglected.  A  design  method  for  estimation  of  scour  depths  at  piers  is  presented  by 
Melville  and  Sutherland  (1988),  based  on  envelope  curves  from  laboratory  data.  The  design 
curve  for  pier  scour  without  debris  accumulations,  developed  by  Melville  and  Sutherland  is 
described  by  the  following  two  equations: 


<Jx 

l.872| 

frr"’ 

fr  ^ 

1  ~ 

—  <  2.6 

1  (8.1) 

1)  ~ 

i.J)J 

U)  y 

ds 

( Y  ^ 

1) " 

2.4 

—  >  2.6 

KJ)  y 

(8.7.) 

where,  ds  =  depth  of  scour  (m),  D  =  bridge  pier  diaraeter(m),  and  Y  -  approach  flow  depth 
(ni) 

These  relationships  are  affected,  however,  by  additional  factors  where  dear-water  •’cour 
conditions  exist,  the  (low  is  relatively  shallow,  and  the  sediment  size  relatively  eoarse.  In  the 
case  of  non-cylindrical  (.ners,  additional  rriuifiplying  factors  arc  applied  to  account  for  piers 
shape  and  alignment. 

Scour  depth  due  to  the  pier  is  calculated  m  the  program  using  the  full  design  method  presented 
by  Melville  and  Sutherland  (1988)  rather  than  the  simplified  version  stated  in  equations  8. 1  and 
8.2.  Total  scour  at  a  bridge  pier  is  defined  as. 


ds  =  kikdkyk„k„D 


(8,3) 


where,  ds  =  total  scour  depth  (m),  ki  =  flow  intensity  factor,  kj  =  sediment  size  factor,  ky  = 
flow  depth  factor,  =  pier  alignment  factor,  k*  =  pier  shape  factor,  and  D  -  pier 
diameter  (m). 

•  Flow  intensity  factor  (ki) 
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where,  LI  mean  approach  flow  velocity  (ro/s'*),  Ua  ='  mean  approach  velocity  at  the  armour 
peak  (m/s  '),  and  Uc  ^  mean  approach  flow  velocity  at  threshold  condition  (m/s’'). 


l/c  -  ti.^  5  751oB^%J+5.53«... 


where,  u*c  ~  critical  shear  velocity  (m/s''),  y  -  total  flow  depth  (m),  djo  --  median  particle 
size  (m). 

=  V^(/w  -  ]  (^  '^) 

where,  9  =  Shield’s  parameter/the  critical  shear  stress  (taken  to  be  0  05  in  the 
program),  p,  =  density  of  the  sediment  (taken  as  2650  kg/m'  in  the  program),  p.,  = 
density  of  water  (1000  kg/m^),  and  g  gravitational  acceleration  (9.81  m/s'^) 


Ua  =  0.8  //,,„5.751og(-J^'^^_^)  +  5.53/ j  if 


UasUc 


Ua  ^  Uc 


iflJa>Uc 


"7^ 


t 


where,  '  critical  shear  velocity  of  annoured  bea,  calculated  as  per  equation  8,7  using 
dso.  rather  than  dso  (m/s’’),  dso.  =  dso  size  of  the  coarsest  a'mour  layer  =  dm«  /  1.8 
m),  and  dmax  =  majdmum  particle  size  for  a  nonuniform  sediment  (m). 

Sediment  size  factor  (kd) 


U-aia-Uc)  , 

it - - - >  1 

Uc 


ku=  l.O 


-0,57  log  2, 


if  D/dso  >25 


if  D/dso  <  25 


(8.10) 


U-{Ua-Uc) 


<  1  and  Og  <  1,3  (where  Og  =  standard  deviation  of  the  grain  size 


distribution  =  d84  /  djo) 


EquationsS.  10  and  8. 1 1  are  calculated  using  d5o» 


(8.12) 


How  depth  factor  (ky) 


ky=  1.0, 


k..  =  0.7; 


if  y/D  >  2.6 
if  y/D  <2,6 


(8.13) 


(8.14) 


•  Pier  aSignment  factor  (L.) 

This  is  determined  using  figure  8,2. 


•  Pier  shape  factor  (k,) 

Pier  shape  factors  are  presented  in  figure  8.3.  The  program  contains  a  look-up  table  of  pier 
shape  factor. 
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Figure  8.2  ;  Alignment  Factor  Ka  (after  Melville  S'  Sutherland,  1988) 
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L=  pier  length;  b  =  pier  diameter 
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Figure  8.3  Pier  Shape  Factors  (after  Melville  &  Sutherland,  1988) 


Sliapc  in  plan 

Length/ 

width 

Tison 

(1940) 

Laursen  and 
Toch(1956) 
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Melville  and  I^ongol  ( 1 992)  present  a  method  for  delernitning  total  scour  due  to  a  bridge  pier 
and  floating  debris  raft,  based  upon  flume  model  data.  The  experimcjital  arrangement  used  by 
Melville  and  Dongol  is  shown  in  Figure  8.4. 

The  trend  in  equation  8.1  was  found  to  al.se  de.se .ibe  srrour  due  to  piers  with  debris 
accumulations  for  values  of  Y/D  <  At  higher  values  of  Y/D  scour  depths  decreased  again 


77 


because  the  proport  ion  of  pier  length  covered  by  debris  decreased,  For  deep  flo  ws  the  effect  of 
debris  became  insignilicant  and  tended  towards  the  value  ds/D  -2.4. 

Analysis  fi'orn  tiic  model  produced  the  following  equation  for  bridge  scour  affected  by  a  debris 
raff  based  upon  an  effective  pier  diameter  (De) : 

1)^.^ - - (8  15) 

where,  Td*  ‘  0,52Td,  Td  -  debris  raft  depth  (m),  Dd  =  debris  rafl  length  parallel  to  the  flow 
(m).,  Y  ■■=  flow  depth  (rn),  and  D  -  pier  diameter  (m), 

T  he  factor  0.5,2  was  determined  by  evaluating  the  limits  of  Td  and  Dd/D  for  the  hypothetical 
case  where  D  is  assumed  to  be  zero  and  the  debris  is  assumed  to  extend  to  the  base  of  the 
scour  hole. 

D  can  therefore  be  substituted  for  De  to  calculate  scour  depth  at  piers  with  debris 
accumulations  using  the  Melville  and  Sutherland  design  method.  Conversely  a  maximum 
allowable  Td  and  Dd  can  be  calculated  by  specifying  an  upper  scour  depth  within  an  acceptable 
factor  of  safety  for  a  given  pier  size. 

Consideration  of  the  likelihood  of  debris  build-up  is  not  addressed  by  Melville  and  Dongol 
(1992)  but  they  do  note,  however,  that  single  cylindrical  piers  are  the  least  likely  shape  to 
accumulate  debris,  and  that  the  free  space  bet  ween  columns  is  seldom  great  enough  to  pass 
debris.  Prediction  of  the  size  of  possible  debris  rafts  accumulations  remains  the  biggest  problem 
for  accurate  factor  of  safety  calculations 
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Figure  8.4  Experimental  Set-up  (after  Melville  and  Dongol,  1992) 


8.3.1  Probability  based  debris  build-up  model 

The  rate  of  acoris  accumulation  <!'  a  bridge  is  difficult  to  quantify.  The  only  method  found  in 
the  literature  is  that  presented  by  Simons  &  I,i  (1979)  in  an  MSc  thesis  by  Callander  (1980) 
entitled  "Fluvial  processes  occurring  at  bridge  sites”. 

According  to  Simons  Ik  Li,  the  trapping  etHcicncv  of  a  bridge  is  determined  by 

I )  (dcarance  beneath  the  bridge 
:i)  Span  lei:gths 

.5)  Size  and  concentration  of  debris  elements 

rhe  following  possible  consequences  are  identified  which  can  result  from  debris  blockage: 

1 )  Mackwate:  etVeets 

2)  Potential  local  ilovv  diversion 

3)  Channel  avulsion 
/  bridge  failure 

Sf  iiions  k  1,1  e.spres'i  tile  v4ume  '  debris  as  a  Fraction  of  the  sediment  yield  (mass),  and 
suggcM  vegetalioii  'le  :is  yi  M  '  !'>■  Ft  an  attempt  to  estimate  the  number  and  volume  of 

tr<:;cs  arrivi.ig  ’*  a  bridge  the  -  •  v(  n  ■(■  of  floodplain  er(i,sion  necessary  to  yield  a  tree, 

and  use  a  repieseiiiali  v'C  tree  SI/  lot  .  iei.'.i  ' 


♦ 


Trees  are  assumed  to  be  cylindrical,  with  a  diameter  Dt  and  a  height  Ht,  The  span  between 
piers  is,  Ls  and  the  clearance  between  the  water  surface  and  the  underside  of  the  bridge  is  C. 
The  chance  that  a  tree  will  be  trapped  depends  on  a  larger  diameter  however,  Db,  which 
represents  either  the  canopy  dimension  or  the  root  zone,  whichever  is  greater  (see  figure  8.5). 
IfHt  >  Ls  the  probability  of  at  least  one  average  tree  being  trapped  is  100%.  The  blocked  area 
is  then  estimated  to  be,  NHtDt,  where  N  is  the  equivalent  number  of  average  trees  assumed  to 
be  trapped  against  the  upstream  face  of  the  bridge  IfHt  <  Ls  a  probabilistic  approach  is  used 
Pt  is  the  probability  of  a  tree  being  trapped  and  as  the  blockage  beneath  a  span  increases  the 
chance  of  other  trees  being  trapped  increases  The  probability  of  the  first  tree  being  trapped  is 
assumed  to  be  a  ratio  of  half  the  tree  diameter,  Db,  to  the  total  wateiway  ai  ea  beneath  a  span, 
Ls(-  Hence: 


ji  Dh^ 
Ls<'  ''  8  7,.v(’ 


(8,16) 


where,  PTl  =  probability  of  the  first  tree  becoming  trapped,  Db  canopy/root  wad  diameter 
(m),  Ls  =  distance  between  piers  (m),  C  =  distance  between  water  surface  and  the 
underside  of  the  bridge  (m),  and  t:  -3.14 

Li  (see  Callander,  1 980)  observed  that  a  tree  caught  on  a  nier  will  in  general  lie  with  its  tmnk 
in  the  direction  of  flow.  A  tree  thus  trapped  offer::  an  area  ol 
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(8.17) 


lO  trap  fuilhcr  uebiiis. 

In  general  when  (m-1)  trees  are  trapped  bvineath  a  span  the  probability  of  an  mill  tree  becoming 
trapped  (PTm)  is: 


PTm 


_ nPh^  /  8 
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(8.18) 


'I'he  probability  of  passing  all  N  I'  trees  from  the  watershed  is: 


(8,19) 


(1-PTl)^ 

The  probability  of  at  least  one  tree  being  trapped  at  a  span(Pl)  is: 

PI  =  l-(l-PTl)''  (8,20) 

where,  N  is  the  equivalent  number  of  average  trees  arriving  at  the  span.  According  to  Li  most 
trees  will  stay  close  to  the  bank,  thus: 


N  =  NT  /  2 


(8.21) 


The  probability  that  m  trees  will  be  trapped  (Pm)  is: 


Pm  -=  [i-(l-PTm)'‘'^'‘'"*‘’']P(m-l)  (8.22) 

On  this  basis  the  probability  of  a  least  m  trees  being  trapped  (for  any  m  <  N)  can  be  estimated. 
To  calculate  Td  and  Dd  requires  an  estimate  of  the  blockage  area.  It  is  assumed  that  debris 
elements  stack  up  and  that  trees  overlap  by  Db/2.  Thus  for  m  trees  trapped  the  percentage  of 
the  waterway  area  which  is  blocked  is: 


%H  lockage  — 
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(8.23) 
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Figure  8.5  Schematic  diagram  of  debris  accumulation  at  bridge  piers 


8.3.2  Program  ("onstrurtion  and  Input  Variables 

The  program  applies  the  Simons  and  I,i  probability  model  and  then  calculates  the  potential  pier 
scour  due  to  the  predicted  debris  mat  size  based  on  the  blockage  area  and  assuming  all  the 
trees  available  in  the  reach  upstream  become  trapped  beneath  a  span. 

Initially  average  (mid-height)  tree  trunk  diameter  (Dt),  a  maximum  tree  diameter  (Db),  (eithei- 
root  wad  or  canopy,  whichever  is  the  larger),  and  average  tree  height  (Ht)  values  arc  entered. 
Next,  the  number  of  trees  approaching  the  bridge  .span  (NT)  is  entered.  Although  Simons  and 
Li  suggest  using  N  NT/2  ii  the  probability  calculatiirns,  the  model  assumes  that  all  the  trees 
available  in  the  upstream  reach  will  pass  through  the  span  in  question.  However,  the  number  of 
spans  (S)  between  piers  (P)  that  are  set  in  the  channel  will  normally  be  S  P  '  1  (counting  the 
two  spans  between  pier  and  river  bank).  It  is  therefore  necessary,  Ibr  an  accurate  assessment  of 
blockage  potential  and  debris  related  scour,  to  calculate  probabilities  for  each  span  individually, 
perhaps  using  a  .simple  division  rule  (N  N'l'  /  S)  Ibr  N  trees  arriving  at  each  span.  It  is  let)  up 
to  the  user  to  make  the  appropriate  adjustments  for  each  span.  NT  can  either  be  estimated  in 
the  held  and  entered  as  a  total  potential  tree  supply  or  it  can  be  c.stimatod  throug,h  calculation 
ol'  potential  bank  failure  in  the  upstream  icach  To  calculate  the  latter  estimate  requires 
knowledge  ol’ riparian  tree  den.sily,  the  length  ol’lhe  reach  in  quc.slion  and  the  potential  bank 
iailuie  width  I'he  lailiire  width  value  can  lic  determined  using  an  appiopiiatc  b.uik  stability 


model  such  as  BURBANK  (developed  under  the  DEC  Monitoring  Project).  The  potential 
number  of  trees  that  will  reach  the  span  is  then  calculated  as: 

NT  =  tree  density  x  failure  width  x  reach  length  x  2  (two  banks)  (8.24) 

Finally,  the  bridge  pier  diameter  (D),  span  between  piers  (Is)  and  average  flow  depth  (Y) 
values  are  entered. 

Calculations  in  the  model  then  proceed  as  follows: 

1)  If  tree  height  is  less  than  the  pier  spacing  the  probability  of  the  first  tree  becoming  caught  is 
calculated,  followed  by  the  probability  of  the  next  tree  becoming  caught,  consecutively.  This  is 
repeated  for  n  trees  up  to  NT. 

In  the  calculation  of  trapping  potential  it  is  considered  that,  the  use  of  the  ratio  of  tree  area  to 
the  entire  area  under  the  span  as  suggested  by  Simons  and  Li,  is  somewhat  inappropriate  as 
tree  capture  is  dependent  only  upon  the  length  of  span  and  diameter  of  tree  given  that  the 
water  level  is  constant.  Deck  elevation  above  the  water  (C)  has,  therefore,  been  substituted  by 
maximum  tree  diameter  (Db)  in  this  model. 

2)  If  tree  height  is  greater  than  span  width  it  is  assumed,  as  outlined  in  the  Simons  and  Li 
model,  that  at  least  one  tree  will  become  trapped  and  thus  all  subsequent  trees  arriving  at  the 
span  will  also  be  caught. 

3)  The  percentage  of  the  channel  cross-sectional  area  that  is  blocked  if  all  the  trees  supplied  to 
the  reach  become  trapped  is  calculated,  as  outlined  in  the  theoretical  model  (see  figure  8.5) 
for  Ht '  Is.  However  if  Ht  >  Is,  Dt  is  substituted  for  Db  and  the  blockage  area  is  calculated 
as: 

square  root  x  blockage  area  =  blockage  depth  (assuming  debris  builds 

up  as  a  square)  x  tree  height  /  span  width  x  flow  depth  x  100  %  (8.25) 

This  calculation  assumes  that  for  Ht  >  Is  all  trees  will  build  up  in  a  square  raft,  when  view 
orthogonal  to  the  flow,  as  oppose  to  parallel  with  the  flow  when  Ht  -  Is. 

4)  The  hydrostatic  pressure  force  on  each  pier  per  unit  width  is  calculated  as  : 


pressure  force  =  bulk  weight  of  water  x 

blockage  depth  x  1  (unit  width)  x  (blockage  depth  /  2). 


(8,26) 


5)  Dynamic  pressure  force  on  the  pier  is  calculated  using  the  following  equation: 


Fd-0.5  pw  ACd/1000 


(8.27) 


where,  Fd  =  Force  on  pier  (kN),  pw  =  density  of  water  (1000  kg/m^),  U  =  approach  flow 
velocity  (m/s  '),  A  -  frontal  area  of  debris  raft  normal  to  the  flow  (ni^),  and  =  drag 
coefficient  of  debris  raft,  A  estimated  value  of  Cu  =  1 .5  is  used  in  the  program. 

6)  Bridge  pier  scour  with  the  debris  accumulation  is  then  calculated  using  the  Melville  and 
Dongol  model  If  Ht  >  Is  the  debris  raft  depth  is  taken  as  the  square  root  of  the  blockage  area 
(assuming  debris  build-up  is  in  a  square)  If  Ht  <  Is,  debris  raft  diameter  is  assumed  to  be  Ht 
because  the  debris  is  aligned  parallel  with  the  direction  of  flow.  Factors  (Ki),  (Ky),  (Kd)  and 
(K,)  are  incorporated  into  the  model.  The  program  requires  the  following  data  to  calculate 
these  factors:  mean  approach  velocity  for  the  design  flood  (U);  median  particle  size  (djo);  the 
largest  particle  size  (d„,4x);  standard  deviation  of  the  particle  distribution  (Og)  =  dsVdso,  pier 
diameter  (D);  angle  of  flow  attack,  pier  dimensions;  and  pier  shape.  The  alignment  factor,  Ku, 
must  be  determined  by  the  user  using  figure  8,2, 

It  should  be  noted  that  the  formula  developed  by  Melville  and  Dongol  for  calculating  debris 
related  pier  scour  was  only  developed  for  floating  debris  accumulations.  However,  it  is 
considered  that  this  formula  can  be  extended  to  debris  accumulations  which  have  their  base 
resting  on  the  channel  bed,  as  the  critical  factor  in  the  calculation  method  is  an  effective  pier 
diameter  which  is,  in  any  case,  extended  to  the  channel  bed  in  the  situation  where  the  debris  is 
floating 

8.3.3  Testing  the  Model 

The  pier  scour  component  of  this  model  has  been  tested  using  field  data  collected  from  a 
number  of  bridges  in  the  survey  reaches  in  northern  Mississippi.  Scour  depths  were  measured 
at  each  site  (during  low  flow  conditions)  on  piers  which  had  significam  debris  accumulations. 
Plate  8  shows  debris  build-up  against  the  piers  of  a  county  road  bridge  over  Fannegusha  Creek. 
The  parameters  required  for  the  model  calculations  were  also  collected.  Bankfull  discharge  and 
channel  dimension  values  obtained  from  Watson  et  al  (196,3)  were  used  in  the  model  to 
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simulate  critical  conditions  because  accurate  discharge  measurements  could  not  be  made  at  the 
time  of  the  survey.  Logistical  constraints  also  prevented  the  calculation  of  realistic  upstream 
debris  loadings  that  could  arrive  at  each  pier  and  so  these  were  represented  by  the  dimensions 
of  the  debris  accumulations  actually  present  at  the  time  of  survey.  Table  8.1  shows  a  summary 
of  the  surveyed  scour  values  and  the  model  results  Figure  8.6  shows  a  plot  of  these  results. 
The  diagonal  line  in  this  plot  represents  a  perfect  match  between  the  actual  measured  scour 
depths  and  those  predicted  by  the  model.  The  graph  shows  how  debris  rafts  increase  scour  by  a 
factor  of  2  to  4.  It  is  evident  that  the  model  significantly  overestimates  scour  due  to  both  the 
pier  and  debris  raft  and  slightly  overestimates  scour  due  to  the  pier  alone.  The  predicted  results 
would,  therefore,  appear  to  be  rather  conservative,  although  this  does  create  a  good  factor  of 
safety. 

The  discrepancy  between  measured  and  predicted  values  may  be  explained  by  the  fact  that 
scour  hole  depths  are  much  greater  under  bankfull  flow  conditions,  than  those  at  low'  flow 
(when  the  measurements  were  made)  due  to  general  scour  hut  are  subsequently  reduced  as 
flows  recede  owing  to  the  deposition  of  the  highly  mobile  sand  and  silt  sediment  load.  To  fully 
validate  the  model  it  will  therefore  be  necessary  to  undertake  fiirther  fieldwork  to  measure  the 
necessaiy  parameters,  including  scour  depth  values  during  bankfull  flow  conditions.  The  model 
must  also  be  validated  using  different  channel  environments  such  as  armoured,  gravel-bed 
rivers. 


TaHc  8.1  Pier  scour  summary  table 
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Figure  8.6  Relationship  between  Measured  and  Predicted  Pier  Scour  Depths 
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9  CONCLUSIONS 


The  following  conclusions  can  be  drawn  from  the  data  and  results  presented  in  this  paper  about 
the  impact  of  LWD  in  degrading,  sand-bed,  rivers: 

•  The  main  debris  input  mechanisms  (outer  bank  erosion  in  bends  and  bank  mass-wasting  in 
degrading  reaches)  are  sp  tially  predictable  processes  which  can  be  used  to  predict  and  map 
the  distribution  of  major  debris  sources  in  the  watershed 

•  Stream  power  per  unit  bed  area  was  found  to  be  a  statistically  significant  predictor  of  debris 
volume.  However,  watershed  area  was  found  not  to  be  a  good  predictor  of  debris  volumes  and 
jam  frequency,  probably  because  the  channels  surveyed  are  undergoing  degradational  cycles 
which  disrupt  simple  hydraulic  geometry  relationships. 

•  Observed  difference.,  in  the  volume  of  debris  and  frequency  of  jams  in  each  Stage  of  the 
Schumm  (1984)  Channel  Evolution  Model  (CEM)  can  be  explained  by  the  suite  of  processes 
which  occur  in  the  CEM  evolution  sequence 

•  The  conceptual  jam  classification  model  (Underflow,  Dam,  Deflector,  Flow  Parallel)  shown 
schematically  ..i  Figure  5,2  has  been  found  to  predict  with  reasonable  accuracy  the  observed 
distribution  of  jam  types  encountered  when  moving  downstream  through  the  drainage  network, 
despite  distortion  of  the  simple  underlying  relationship  in  the  model  (ratio  of  average  channel 
width  to  ‘key’  debris  length)  by  channel  evolution  processes. 

•  Jam-induced  sedimentation  and  scour  values  have  been  found  to  have  statistically  significant 
trends  when  correlated  with  watershed  area  These  trends  can  be  explained  because  they 
correspond  to  the  changing  flow  processes  and  consequent  scour  and  sedimentation  patterns, 
induced  by  each  jam  type  encountered  through  the  watershed  IJndeiilow  jams,  in  small 
watersheds,  interfere  very  little  with  the  flow  and  therefore  do  not  have  a  high  scour  and 
sediment  retention  potential.  Dam  type  jams,  found  further  downstream,  cause  large  volumes 
of  seuiiiiciU  to  be  stored  in  backwaters,  but  also  cause  plunge  pool  scour  Downstream  from 
Dam  type  jams  at  Deflect  ;  jams,  sediment  is  stored  in  the  lee  of  jams  as  bar  deposits,  but  the 
jams  also  cause  flow  to  impinge  on  one  or  both  banks,  resulting  in  localised  scour  and  bank 
collapse.  Further  downstream  still,  Flow-Parallel  debris  block  the  flow  much  less  so  that 
energy  dissipation  and,  therefore,  sediment  retention  and  scour  are  less  significant 

•  The  net  sediment  balance  results  shown  in  Table  5.2  support  the  arguments  that  debris  helps 
to  accelerate  sedimentation  processes  in  these  highly  unstable  channels  and  that  jams  may  even 
act  as  temporary  grade  control  structures  l-WD  therefore  accelerates  channel  equilibrating 


processes  in  stages  4  and  5  of  the  Schumm  et  al  (1984)  Channel  Evolution  Model.  Jams 
should  not,  therefore,  be  removed  unless  they  threaten  structures  and  addition  of  LWD  is 
recommended  to  aid  sediment  retention  and  improve  aquatic  habitat  in  debris-starved 
degrading  reaches 

•  Meaningful  conclusions  concerning  short-term  jam  stability  can  only  be  made  when  the 
survey  data  is  compared  with  corresponding  discharge  data  for  the  one  year  interval,  which 
will  show  the  magnitude  and  fiequency  of  events  that  the  jams  were  able  to  survive.  The 
relevant  data  are  not  available,  as  yet,  to  permit  this  analysis.  Given  the  warm,  humid  climate, 
and  highly  erodible  nature  of  these  channels  it  is  likely,  however,  that  jams  are  more  transient 
features  in  North  Mississippi  than,  for  example,  in  the  Pacific  North  V/est  because  debris  will 
decompose  more  quickly  and  become  bypassed  by  bank  erosion  and  bed  scour.  A  long-term 
monitoring  programme  is  required,  however,  to  verify  this  assertion 

•  Four  LWD  analysis  models  are  presented  in  this  work  which  have  been  tested  and  verified 
using  field  data  and  are  recommended  for  channel  design  and  management  They  are: 

1)  The  Shields  &  Smith  (1992)  LWD  flow  resistance  model. 

This  model  predicts  Darcy-Weisbach  friction  factor  due  to  LWD  jam  accumulations.  It  has  a 
sound  theoretical  basis  and  has  shown  a  reasonable  correlation  between  measured  and 
computed  friction  factors. 

2)  The  Gippel  et  al.  (1992)  LWD  induced  flow  alflux  model. 

This  analysis  method  has  been  developed  from  flume  experiments  and  theoretical  hydraulic 
considerations,  The  analysis  method  is  presented  in  a  manual  format  in  Appendix  C.  This 
method  of  backwater  or  afflux  calculation  due  to  individual  items  of  LWD  can  be  used  as  a  tool  to 
help  determine  whether  the  afflux  reduction  due  to  LWD  removal  would  have  a  significant,  positive 
impact  according  to  the  perceived  management  requirements  or  whether  LWD  could  be  left  in  place 
perhaps,  re-orientated,  lopped  or  even  re-introduced  where  sympathetic  rehabilitation  management 
is  desirable,  without  significant  effect  on  high  in-bank  stages 

3)  The  Debris  Management  Program  (DMP). 

This  computer  model,  presented  on  the  disk  enclosed,  is  designed  as  a  management  too!  for 
predicting  the  geomorphic  impact  of  LWD  on  charmel  processes  and  morphology  throughout  the 


<» 
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v/atershed.  The  model  has  been  developed,  tested  and  modified  using  empirical  analysis  of  field  data 
gathered  from  23  river  reaches  in  Northern  Mississippi  which  are  undergoing  channel  adjustment 
caused  by  base  level  change.  This  tool  is  recommended  for  geomorphic  analysis  and  management  of 
LWD  in  channel  environments  similar  to  those  from  which  the  model  has  been  developed. 

4  4)  The  Debris  at  Bridge  Pier  Prediction  Program  (DBP3). 

This  computer  model,  also  included  on  the  enclosed  disk,  determines  the  probability  of  debris  build 
up  at  bridge  piers,  the  consequent  maximum  debris  enhanced  potential  pier  scour  and  the  change  in 
^  hydrostatic  and  dynamic  pressure  forces  on  the  pier.  The  model  is  developed  from  research 

presented  by  Melville  and  Sutherland  (1988),  Melville  and  Dongol  (1992)  and  Simons  and  Li  (see 
Callander,  1980).  The  model  has  been  tested  using  field  data  collected  from  a  number  of  bridges  in 
Northern  Mississippi  and  demonstrates  that  the  empirical  and  theoretical  assumptions  of  the  model 
provide  a  good  factor  of  safety  for  bridge  pier  design  analysis  where  LWD  build-up  is  a  potential 
hazard 
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APPENDIX  A 

1 

LWD  MAMAGEMENT  PROGRAM  GIS  DATA  INPUT  SYSTEM 

i 

1.1  Introduction 

i:: 

This  GIS  (Geographical  Information  System)  user  interface  demonstrates  how  a  GIS 

(ARCANFO)  can  be  used  in  conjunction  with  external  programming  to  assist  in  management 

i 

1 

practices  associated  with  erosion  control. 

The  GIS  is  used  to  calculate  variables  of  sediment  type,  land  use,  and  sub- watershed  are  which 

are  the  input  variables  for  the  DMP.  The  GIS  has  been  constructed  for  the  Abiaca  Creek 

t 

w'atershed,  located  within  the  DEC  study  area  in  Northern  Mississippi.  The  operation  of  the 

program  is  set  up  to  facilitate  interactive  selection  of  stream  section,,  within  the  watershed. 

The  soil  type  and  landuse  variables  are  calculated  by  combining  polygon  coverage’s  and 

1.. 

overlaying  them  onto  the  stream  network.  The  watershed  area  is  calculated  through  the  use 

t 

i;'  • 

of  a  digital  elevation  model,  which  is  also  used  to  derive  the  stream  network.  The  system 

macro  language  (AML)  of  the  GIS  is  used  to  automate  data  input  in  the  DMP.  A  pull  down 

menu  is  employed  to  create  a  user  friendly  windows  style  environment  which  guides  the  user 

t 

1 

through  the  steps  necessary  to  obtain  management  information. 

1.2  Installation 

t 

The  GIS  is  in  UNIX  ARCANFO  format,  and  is  contained  in  a  directory  named  “Project” 

1 

The  program  can  be  obtained  upon  request  from  nick(l^geography.noUingham.ac.uk  and 

will  be  sent  via  FTP  or  on  DAT  tape  “Project”  requires  about  98  MB  of  hard  drive  space  to 

install 

( 

r 

The  “Project”  directory  should  be  copied  into  the  same  directory  as  ARC  if  ARC  is  not  in  the 

command  path 

To  run  the  GIS  load  ARC  At  *he  ARC:  prompt  type  &RUN  Project  to  activate  the  AMI. 

(' 

(ARC  Macro  Language)  driven  user  interface. 

1.3  Methodology 

r 

A  fundamental  requirement  for  this  utility  was  too  extract  these  variables  automatically  with 

limited  user  interaction  and  run  the  DMP,  all  within  the  framework  of  the  GIS 

Figure  1 . 1  outlines  the  method  of  approach  from  variable  extraction  to  the  final  output  from 

the  DMP, 

1 

1(H> 

The  section  of  stream  which  the  user  is  querying  is  interactively  selected  by  using  the  screen 
cursor  operated  by  the  mouse.  This  provides  the  co-ordinates  for  the  analysis  from  other  data 
layers  which  are  all  geometrically  aligned.  AJIC/INFO  provides  interactive  selection  facilities 
which  enable  the  user  to  select  points  on  the  display  screen.  The  stream  network  coverage 
contains  attribute  information  detailing  the  soil  types  and  landuse  for  each  arc  section  of 
stream.  The  original  stream  coverage  data  contained  no  attribute  information  relevant  to  the 
study.  A  function  command  called  IDENTITY  enables  arcs  (lines)  to  be  overlaid  onto  polygon 
coverages  and  relate  the  attributes  from  the  polygon  coverage  to  the  arcs  coverage,  where  the 
arcs  overlay  and  intersect  the  polygons.  By  converting  the  grid  layers,  soil  types  and  landuse  to 
polygon  coverages  this  IDENTITY  command  is  then  used  to  attach  the  attribute  information 
from  these  two  data  layers  to  the  stream  netv/ork.  ARC/INFO  facilitates  conversions  from 
raster  to  vector,  and  the  GRIDPOLY  command  is  used  to  convert  the  two  raster  files,  Using 
the  IDENTITY  command  the  stream  network  can  be  overlain  onto  the  soil  polygon  coverage, 
then  the  resulting  new  stream  network  coverage  'vith  soil  attributes  is  overlain  onto  the 
landuse  polygon  coverage  resulting  in  a  stream  coverage  which  has  attached  soil  type  and 
landuse  attributes,  The  methodology  of  this  process  is  summarised  in  figure  1 .2. 

An  important  requirement  for  the  success  of  this  project  is  that  the  stream  network  is 
geometrically  true  to  the  elevation  grid  for  the  calculation  of  the  watershed  area  upstream  of 
the  query  point.  The  interaction  with  the  stream  network  is  the  sole  method  of  producing  a  co¬ 
ordinate  reference  point,  which  is  used  as  the  starting  point  or  watershed  pour  point  for  the 
area  calculation  from  the  DEM.  However,  the  likelihood  that  eveiy  section  of  the  original 
vector  scanned  stream  coverage  being  a  precise  match  to  the  DEM  is  small.  If  it  does  not 
match  precisely  the  analysis  method  would  failure  because  the  watershed  area  calculated  would 
be  incorrect  For  this  reason  the  stream  network  is  produced  from  the  DEM,  and  converted  to 
a  vector  coverage  through  the  command  GRIDLINE.  The  same  procedure  is  used  to  produce 
a  stream  coverage  with  soil  and  landuse  attributes. 
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Figure  1.1  Flow  diagram  of  GIS  operation 


1.3.1  The  Digital  Terrain  Model 

A  common  problem  in  using  DEM  matrices  for  detecting  linear  features  such  as  stream 
networks  is  that  of  depressions  (sinks)  in  the  digital  surface  caused  by  noise,  The  noise  may 
result  fiom  short  range  variations  on  the  digital  land  surface  or  as  a  result  of  the  quantitization 
of  the  original  data.  Ridges  and  drainage  courses  may  be  missed  because  the  grid  is  too  coarse, 
ARC/IMFO  can  detect  these  sinks  and  fill  them  A  sink  is  defined  as  a  cell  whose  elevation  is 
less  than  or  equal  to  that  of  its  eight  neighbours  when  passing  a  kernel  of  3x3  cells  over  the 
elevation  grid  (Marks  1984).  ARC/INFO  uses  Marks  method  to  detect  these  sinks  and  by 
setting  a  minimum  elevation,  which  in  this  case  was  83ft,  these  sinks  can  be  filled  up  to  that 
minimum  height.  The  command  FILL  was  used  with  the  options  of  SINKS  and  a  z  value  set  at 
83,  From  here  the  depressionless  DEM  was  used  to  create  matrices  necessary  for  further 
analy  ;es.  The  command  FLOWDIRECTION  creates  a  flow  direction  matrix  from  the  DEM, 
which  is  achieved  by  computing  the  local  gradient  and  aspect  for  every  cell.  This  flow  matrix  is 
used  tc  determine  the  boundary  of  catchment  or  selected  sub-catchment.  To  find  the  stream 
network  from  the  DEM  a  flow  accumulation  grid  needs  to  be  produced.  The  method  behind 
this  is  an  algorithm  for  which  each  cell,  compares  its  altitude  with  its  eight  neighbours  within  a 
3x3  kernel.  The  lowest  neighbour  is  flagged,  and  the  amount  of  water  (which  is  expressed  as  a 
function  of  the  number  of  cells  traversed  and  the  area  of  the  cell)  is  carried  over  to  that  cell. 
The  kernel  is  moved  to  the  lowest  neighbour,  and  the  process  is  repeated.  The  ARC  command 
that  is  used  to  carry  out  this  process  is  FLOW  ACCUMULATION,  A  stream  network  can  then 
be  delineated  using  the  output  from  this  flowaccumulation  grid.  Flow  accumulation  in  its 
simplest  form  is  the  number  of  up  slope  cells  that  flow  into  each  cell.  By  applying  a  threshold 
value  to  the  results  of  this  FLOWACCUMULATION  using  a  grid  algebraic  expression,  a 
stream  network  can  be  defined.  An  expression  was  used  to  create  a  grid  where  the  value  1 
represented  a  stream  network  on  a  background  of  no  data,  the  expression  was: 

stream  network  =  con  (flowaccumulation  >  100,  1) 

This  assigns  the  value  1  to  all  cells  with  more  than  100  cells  flowing  into  them,  and  all  other 
cells  are  assigned  no  data.  A  final  process  that  was  carried  out,  was  to  apply  Strahler’s  method 
of  stream  ordering,  achieved  through  the  use  of  the  command  ORDER.  The  reason  for 
creating  this  grid  was  to  create  a  visual  guide  to  the  user,  the  network  is  colour  coded,  white 
representing  a  stream  order  of  1 ,  while  bright  pink  was  a  stream  order  of  6. 
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Once  the  necessary  data  had  been  constructed  irs  the  database  a  set  of  Arc  Macro  Language 
(AML)  statements  was  used  to  extract  the  variables  necessary  for  input  into  the  DMSP.  A 
command  named  RESELECT  enables  the  user  to  interactively  select  a  point  on  the  stream 
network  displayed  on  the  screen.  This  function  then  selects  all  the  attributes  associated  with 
that  one  arc  which  has  been  selected.  These  attributes  are  stored  in  ;  named  FNFO  file  (TRY) 
created  by  the  RESELECT  function.  The  next  step  is  to  calculate  the  upstream  watershed  area 
from  that  selected  point.  ARCANFO  command  WATERSHED  allows  inter;  ;tive  selection  of  a 
pour  point  from  the  screen.  By  selecting  at  the  same  point  as  before  (the  RESELECT  function 
leaves  a  small  box  on  the  screen  marking  the  previous  selection  point,  the  user  then  selects  the 
point  again),  the  WATERSHED  function  creates  a  new  file  which  represents  the  sub¬ 
watershed  only.  The  count  attribute  in  the  new  sub-watershed  file  represents  the  number  of 
grid  squares  which  make  up  the  area  of  the  sub- watershed.  In  the  ARC  module  TABLES,  new 
items  can  be  created  for  this  file,  then  using  the  CALCULATE  function  within  TABLES  a  data 
value  is  assigned  to  represent  the  number  of  squares  (one  grid  square  -  0.000036  miles^) 
Two  files  are  therefore  created  in  the  analysis,  one  containing  the  values  representing  the  soil 
type  and  landuse  at  the  selected  point  on  a  stream,  and  one  file  containing  the  upstream 
watershed  area  from  that  point. 

In  TABLES  using  the  function  UNLOAD  the  relevant  items  from  these  files  are  selected  and 
put  into  a  text  file  which  list  the  landuse,  watershed  area  and  soil-type  respectively.  This  text 
file  is  the  input  file  into  the  Debris  Management  Program. 

At  the  arc  prompt  the  function  TASK  activates  a  program  outside  the  ARC/INFO 
environment.  This  TASK  command  is  used  to  execute  the  Debris  Management  Program  which 
imports  the  text  file  containing  the  input  variables.  The  DMP  creates  an  output  text  file  which 
detculs  the  soil  type,  landuse  and  stream  width  (calculated  from  the  watershed  area)  and  the 
form  of  management  which  should  take  place  at  that  section  of  the  stream  network  in  the 
Abiaca  Watershed. 

1.3.2  The  Menu  Driven  Interface. 

The  final  part  to  the  method  for  this  project  was  to  tie  all  the  data  layers  and  interactive 
analyses  functions  together  in  a  menu  driven  interface  The  type  of  menu  chosen  was  the 
pulldown  version  which  is  perhaps  the  most  familiar  to  Personal  Computer  users.  An 
ARC/INFO  menu  file  is  a  text  file,  created  in  the  operating  systems  text  editor.  The  format  of 
the  menu  determines  how  the  menu  is  to  be  displayed  and  which  choices  are  to  be  included, 
and  what  action  is  taken  when  a  selection  is  made.  Table  1.1  shows  the  format  for  this  menu 
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file  with  explanatory  text  written  in  red  The  AML  files  are  described  is  section  1 ,4.  Figure  1,3 
shows  a  schematic  diagram  of  the  outlining  operation  of  the  menu  interface  and  how  it  links 
wth  the  GIS  database. 


Table  1.1  Menu  Interface  Commands 


Draw 

""'‘main  menu  choice,  with  sub  menu  choices  beiow** 

Landuse 

&R LAND 

**  runs  the  AML  land  (displays  landuse)** 

Soils 

&R  SOILS 

**  mns  the  AML  soils  (displays  soils)** 

Streams 

&R  STREAMS 

**  runs  the  AML  streams  (displays  streams)** 

Roads 

&R  ROADS 

**  runs  the  AML  roads  (displays  roads)** 

List  Attribute 

**main  menu  choice,  with  sub  menu  choices  below** 

Landuse 

LIST  AB-LAND-GRID.VAT 

**lisls  attribute  file** 

Soils 

LIST  AB-SOiL-GRID  VAT 

**lists  allribute  111c** 

Streams 

LIST  FINAL-NET  AAT 

**!ists  attribute  file** 

Check  Count 

LIST  TEST.VAT 

**lists  attribute  file** 

Soil  &  Landuse 

LIST  TRY 

'**lisls  attribute  tile** 

Analysis 

**mcnu  choice,  with  ,sub  menu 

choices  below** 

Help 

&POPUP  ANL-HELP.TXT 

*  ""displays  Itelp  text  file** 

Zoom 

&R  ZOOM.  AML 

**runs  the  AML  /.oom** 

Calculate  Variables 

&R  del,  AML 

**runs  the  AML  del** 

Run  Program 

**main  menu  choice,  with  sub  menu  choices  below** 

Run 

*R  FILE,  AML  **runs  the  AML  file** 

Clear 

**nuiin  menu  choice  v/hich  clears  current  display** 

Quit 

&RETURN  **qiiits  the  menu  display** 

Draw: 

t 

provides  a  list  of  ctwerages  to  be  drawn  on  the  screen 

List  Attributes: 

provides  a  list  of  coverages  form  which  the  attributes  associated 
with  them  can  be  chosen. 

Analysis: 

• 

activates  the  stream  analysis  siection  of  the  project,  provides  a 
series  of  help  menus  to  guide  the  user  through  the  process. 

Run  Program: 

after  tlie  stream  variables  have  been  collected  the  DMSP  can  be 
activated 

'  Clear: 

clear  the  current  display  on  the  screen 

Quit: 

Quits  the  menu  and  returns  to  the  arc  prompt. 

i 
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Figure  1.3  Schematic  of  menu  operation 


The  menus  is  set  up  in  such  a  way  that  choosing  a  subject  activates  an  AML  program.  This  is 
set-up  so  that  the  user  need  not  be  aware  vof  the  processing  that  is  taking  place  when  they  make 
a  choice,  The  AML  programs  are  a  list  of  on  line  commands  which  activate  various 
ARC/INFO  functions  These  AMLs  could  all  be  accomplished  by  typing  in  the  numerous 
commands  line  by  line  Instead  an  ARC/INFO  fianction  &WATCH  was  used  to  record  all  iine 
entries  made,  This  creates  a  watch  file  which  contains  all  the  commands  entered  at  that  time, 
these  watch  files  are  then  converted  to  AML  through  the  command  CWTA  (convert  watch  file 
to  AML),  By  setting  up  watch  files  whilst  doing  the  individual  process  by  line  entry,  separate 
small  command  files  were  created  which  would  carry  out  a  whole  process  just  by  activating  the 
particular  AML  program  file  in  question  These  AML  files  are  then  brought  together  through 
the  menu  interface. 

1.4  AML  Commands 

Listed  below  are  some  of  the  key  ARC/INFO  commands  used  in  the  project  with  brief 
descriptions  of  their  function. 

IDENTITY :  computes  the  geometric  intersection  of  two  coverages. 

All  features  of  the  input  coverage,  as  well  as  those  features  of  the 
identity  coverage  that  overlap  the  input  coverage,  are  preserved 
in  the  output  coverage. 

RESELECT;  selects  a  set  of  records  from  the  specified  object  i  .e.  arc 

attributes. 

creates  a  grid  identifying  all  sinks,  or  areas  of  internal 
drainage 

fills  sinks  to  the  value  of  the  lowest  boundary  cell  in  the 
watershed. 

FLOW  ACCUMULATION;  creates  a  grid  of  accumulated  How  to  each  cell,  by 

accumulating  the  wsight  for  all  cells  that  flow  into  each 
downslone  cell 

FLOWDIRECTION;  creates  a  grid  of  flow  direction  from  each  cell  to  its 

steepest  downslope  neighbour. 

WATERSHED;  creates  a  grid  of  the  upslope  area  contributing  flow  to  a 

given  location 

STREAMORDER;  creates  a  grid  of  streams  characterising  the  stream 

network  based  upon  their  number  of  tributaries. 
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HILLSHADE:  creates  a  shade  relief  grid  from  a  grid  by  considering  the 

illumination  angle  and  shadows. 

GRIDPOLY :  converts  a  grid  to  a  polygon  coverage.  Polygons  are 

built  from  groups  of  contiguous  cells  having  the  same 
cell  value. 

1.4.1  Arc  Macro  Language  Programs 

Programs  1  to  4  are  display  functions 


1. **land** 
clear 

mape  ab-land-grid 
map  land,  map 
map  end 
&return 

2.  **soil** 
clear 

mape  ab-soil-grid 
map  soil,  map 
map  end 
&relum 

3.  * ’"streams'"* 
mape  ab-soil-grid 
arcs  final-net 
&return 


4*  "roads** 
mape  ab-soil-grid 
arcs  ab-roads 
&return 


""sets  the  map  extent  to  the  landuse  grid  file** 

*  "displays  a  map  file  which  has  title  and  key** 

*  "closes  map  file** 


*  "displays  arcs  coverage** 


Programs  5  to  9  are  data  handling  functions 
5  *  "del. ami** 

&if  [exists  TRY  -INFO]  &then  ""this  is  an  if-then  type  statement  program 

&do  IF  an  INFO  file  called  TRY  exists  delete  it, 

&sv  delstat  =  [DELETE  TRY  -INFO]  if  not  (ELSE)  continue  and  run 
analysis,  ami** 

&TYPE  TRY  FILE  DELETED 

&end 

&else 

&r  analysis. ami 

&return 

&ranalysis 


OS) 
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6. **file.aml*8 

q 

tables 

additem  test  vat  area  14  16  n3 
test,  vat** 
select  try 

urJoad  dump.txt  land-code 

land-code** 

select  test,  vat 

calculate  area=  count  *0,00036 
unload  dump.txt  area 
q  stop 

&sv  X  =  [task  3res  ‘file  =0’] 
&pc)pup  output.txt 
&stop 

7. **zoom.aml** 

&POPUP  ZOOM-HELP,  TXT 
mape  ab-order 

mape* 

clear 

gridshades  ab-order 
&RETURN 


**this  opens  up  the  module  TABLES** 

*  *  creates  new  item  AREA  in  the  file 

**  selects  INFO  file  try** 

**creates  text  file  fump.txt  inputs 

**selects  INFO  file  test.vat** 

**calculates  area  value** 

**inputs  area  into  dump  txt** 

** quits  Tables** 

**  activates  DMSP** 

*  *displays  the  output  file  from  the  DMSP*  * 


**zooms  in  on  analysis  section  of  streams 
the  *  indicates  this  is  a  user  interaction** 


8.**ptojc-ct.aml** 

&termind  9999 

rm  dump  txt 

rm  output.txt 

&Fullscreen  &Popup 

display  9999  3  position  1 1  screen  1 1 

DISPLAY  POSITION  LR 


**this  is  the  starting  AML  which  sets  up  the 
display  environment  in  ARC/INFO  deletes 
previous  text  files  and  then  activates  the  menu 
file  first  menu** 


grid 

mape  ab-land-grid 

&Menu  first  menu  &PULLDOWN  &  position  &ul  &screen  &ul  &size  &frame  700  75 
&  STRIPE  ‘Stream  Erosion  Management  System’ 


9  ♦♦analysis  ami** 

&POPUP  ANALYSIS  TXT 
CLEARSELbCr 
SEARCHTOLERANCE  1 50 
RESELECT  FINAL-NET  ARCS  ONE* 

INFOFILE  FINAL-NET  ARCS  TRY  SOIL-CODE  TYPE  LAND-CODE  '.AND-TYPE 
LIST  TRY 

&POPUP  SHED  TXT 
KILL  TEST  ALL 

TEST  =  WATERSHED  (ab-net-flow,  SELECTPOINT  (ab-elevation,*) 

LIST  TEST  VAT 
&POPUP  CHECK  TXT 
clearselect 
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clear 

MAPE  ab-order 

SE ARCHTOLERANCE  AUTOMATIC 
&RETURN 


1.5  Data  Source 

The  original  DEC  Project  data  resides  on  an  Intergraph  6080  v/orkstation,  and  access  to  the 
database  is  made  with  Intergraph  GIS  software  The  purpose  of  the  engineering  database/GIS 
is  to  serve  as  a  repository  for  all  design,  analysis,  and  monitoring  data  collected  on  the  DEC 
Project.  It  is  still  in  a  development  stage  but  when  completed  it  is  anticipated  that  the  database 
will  contain  design  data  for  all  project  features  such  as  low  and  high-drop  structures,  bank 
stabilisation  structures,  floodwater  retarding  structures,  channel  improvements,  levees,  riser 
pipes,  and  box  culverts.  This  design  data  will  be  complemented  with  date  coverages  suitable 
for  various  hydrological  analyses  such  as  DBMS’s  and  soil  type  data. 

The  database  contains  1:24000  digitally  scanned  USGS  quadrangle  maps  and  OEM’s  for  all 
the  DEC  Project  watersheds  All  major  tributaries  and  highways,  which  have  been  obtained 
from  1:100000  USGS  Digital  Line  Graph  (DLG)  files  are  incorporated  into  the  database  in 
vector  format,  these  maps  were  again  digitally  scanned  rather  than  digitised.  Spot-view  satellite 
photography  has  been  incorporated  into  the  database  and  is  used  as  a  visual  reference  for  all 
DEC  Project  features.  Landuse  and  soil  type  data  for  the  DEC  watersheds  are  incorporated  in 
the  database  on  a  1-acre  grid.  Elevation  data  is  present  on  a  100ft  grid 
The  data  that  were  made  available  for  this  project  were  the  elevation,  landuse  and  soil  type 
grids,  which  were  in  Geographical  Resource  Analysis  Support  System  (GRASS)  format.  The 
vector  files  were  visual  files  only  and  had  no  specific  attributes  attached  to  them.  The 
watershed  chosen  for  this  study  was  the  Abiaca  Creek  Watershed  which  is  one  of  the  fifteen 
watersheds  within  the  DEC  Project  It  was  chosen  because  a  large  number  of  field 
investigation  si^es  were  located  in  this  watershed,  enabling  the  model  to  be  tested  against 
ground  data. 

To  import  the  GRASS  files  into  ARC/INFO  the  image  (GRASS)  file  is  converted  to  GRID 
format  by  initialising  the  command  IMAGEGRID,  which  is  one  of  ARC/INFO’s  raster 
conversion  programs  The  IGDS  files  were  translated  into  ARC/INFO  coverages  with  the 
IGDSARC  command.  Before  the  IGDS  ARC  command  was  implemented  the  command 
IGDSINFO  was  used  which  provided  a  summary  of  what  the  IGDS  file  contains.  IGDS  files 
can  be  multi-layered  coverages  and  it  is  therefore  necessary  to  identity  the  specific  layers 
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contained  in  the  file  for  conversion.  IGDSINFO  provided  a  breakdown  of  the  total  number  of 
elements,  the  names  of  the  elements  and  how  many  of  each,  plus  a  breakdown  of  level,  style, 
colour,  and  weight  and  the  number  of  occurrences  for  each.  The  road  and  the  stream  networks 
were  all  present  in  one  IGDS  file,  which  also  contained  data  on  hydrological  structures.  All 
the  layers  and  options  that  were  entered  during  one  IGDSARC  conversion  were  converted  to 
the  same  ARC/INFO  coverage.  The  roads  were  in  separate  layers  with  respect  to  thiii  class 
i.e.  major  or  minor.  All  the  layers  which  represented  these  roads  v/ere  entered  during  one 
IGDSARC  conversion  and  thus  made  up  one  ARC/INFO  coverage  representing  the  roads  in 
the  Abiaca  watershed  The  same  process  was  carried  out  to  create  the  stream  network,  shown 
in  plate  1 .  Plates  2,  3  and  4  display  the  soil  type,  landuse,  and  DEM  respectively. 
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Plate  1  :  Stream  Network 
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I'ht'  roads  which  arc  in  vector  format  arc  literally  rrnly  a  visual  Ic  for  indicating  where  on  a 
stream  a  bridgi  occurs  (an  overlay  of  the  two  networks  indicates  a  bridge  at  any  intersection) 
These  are  potentially  vulneiabic  areas  when  the  debris  build  up  is  indicated  causing  possible 
.lilure  of  the  bridges  due  to  excessive  stream  bed  scour  By  carrying  mit  the  stream  analysis 
upstream  of  the  intersection  (loinl  i  e  a  bridge,  these  bridge  sites  can  be  analysed 
The  stream  network  lias  been  created  from  the  DbM  It  matched  with  the  original  digitally 
scanned  vector  network  in  general  shape,  the  mam  channels  being  very  much  alike  bar  mote 


minor  tributaries  were  apparent  in  the  derived  network,  although  tli  ■  vector 


scanned  from  a  1  lOOOtM)  llSCiS  map  which  would  not  have  |iicked  up  these  very  small 


tiibutaries 


The  soil  type  map  shows  only  three  types  ol  soil  in  the  watershed  These  soil  types  are 
classitu'd  by  their  particle  si/e  which  is  the  requirement  for  the  Debiis  Management  Support 
l^rogram  In  the  top  two  thirds  of  the  watershed  the  soil  type  is  of  gravel  '  oniirosltion  with  an 
average  particle  si/e  of  2mni  b'ollowing  this  /one  is  a  band  of  line  sandy  soil  with  jrarticle  si/e 
oft)  25  mm,  followed  by  a  narrow  band  of  coarser  sarulv  soil,  paiticle  si/e  I  5mm,  then  again  a 
reg.ion  ot  Tmci  sandy  soil  before  finishing  with  a  small  area  oi’ coarse  saiuly  soil  at  the  outlet  oi' 
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the  watershed.  To  compile  the  attributable  information  in  the  stream  network,  this  grid  was 
converted  into  a  polygon  coverage  of  five  polygons. 

The  landuse  map  shows  three  types  of  landuse  in  the  watershed.  The  small  nose  iike  extension 
from  the  main  body  of  the  watershed  in  the  western  region  is  almost  entirely  agricultural  land; 
this  is  where  the  watershed  reaches  the  end  of  the  bluff  and  stretches  out  into  the  floodplain. 
The  main  b^  dy  of  the  watershed  is  generally  forest  interspersed  by  agricultural  land  and  small 
lakes.  The  original  raster  image  contained  twelve  landuse  types,  fores*,  lakes  and  ten  differing 
types  of  agricultural  land.  When  this  original  image  was  converted  to  a  polygon  coverage, 
ARC/INFO  was  unable  to  complete  the  conversion  because  the  polygons  created  exceeded 
10,000  which  is  ARC/TNFO’s  limit  The  raster  image  was  therefore  reclassified,  all 
agriculti-ral  types  were  classified  as  one  type.  The  significant  variable  value  in  the  DMP  is 
whether  the  landuse  is  forested  or  not,  different  types  of  agricultural  land  are  not  distinguished 
and  therefore  one  value  was  sufficient.  The  grid  to  polygon  conversion  then  created  a  more 
manageable  number  of  1 25 1  polygons. 

The  DEM  is  colour  shaded  based  on  elevation  which  enhances  the  readability  of  the  map,  It  is 
very  clearly  shown  where  the  flood  plain  gives  way  to  the  bluff,  where  the  pale  greens  change 
to  pale  pinks.  A  marked  feature  of  this  map  is  the  apparent  visibility  of  the  main  drainage 
channels,  beginning  from  the  bluff  line  where  the  Abiaca  Creek  protrudes  out  into  the  flood 
plain 

1.5  Program  Operation 

Plate  5  shows  the  menu  interface  displaying  the  landuse  map  with  the  streams  (in  blue)  and 
roads  (in  yellow)  overlaid.  This  is  achieved  under  the  DRAW  option  on  the  menu  title  bar.  This 
provides  the  user  with  the  ability  to  carry  out  some  visual  interpretation  such  as  identifying 
bridging  points  before  they  choose  to  analyse  a  specific  area.  When  the  ANALYSIS  choice  is 
selected  on  the  menu  title  bar  a  ext  box  appears  encouraging  the  user  to  zoom  in  on  t  MV 
for  analysis,  because  it  is  difficult  to  select  an  arc  at  such  a  small  visual  scale  When  zooming  in 
on  an  area,  the  user  interactively  drags  out  a  small  box  around  the  chosen  section,  once  the  box 
is  defined  the  screen  clears  and  then  redisplays  the  selected  area  as  shown  in  Plate  6,  After 
zooming  in  on  the  chosen  section  the  user  then  selects  the  ‘analyse  streams’  choice  in  the  sub 
menu  of  the  ANALYSIS  option.  The  text  box  shown  in  Plate  6  explains  to  the  user  how  to 
activate  the  analysis  by  placing  the  cross  hairs  onto  the  stream  section  of  their  choice.  Once  a 
choice  is  made  a  small  box  indicates  the  position  and  a  text  box  is  displayed  which  explains  to 
the  user  to  re-select  the  box  and  then  press  “9”  to  activate  the  watershed  calculation  Further 


text  boxes  are  displayed  after  this  if  the  selection  process  has  been  unsuccessful,  i.e.  the 
positioning  of  the  selection  point  failed  to  select  an  arc,  or  the  watershed  area  was  far  too  low 
which  would  indicate  that  the  selection  point  was  not  in  the  grid  square  of  the  stream  channel 
on  the  DEM.  Once  this  process  has  been  completed  the  user  activates  the  DMP  by  selecting 
RUN  PROGRAM  on  the  menu  title  bar.  This  activates  processing  in  TABLES  before 
activating  the  DMP  through  the  TASK  function.  Once  the  DSP  has  processed  the  variables  a 
results  text  box  is  automatically  displayed  on  the  screen  (Plate  7) 
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APPENDIX  B 


Large  Woody  Debris  Formation  Survey  used  by  Smith  and  Shields  (1992) 


Wictth-Pcipcndiciilar  to  Flow  Direction 
W<B/4  lB/4<W<B/2  I  B/2<W''B 


Width-Pcrocndicular  to  Flow  Direction 
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TYPE  B  :  RAMP 
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Direction 
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TYPE  D  :  STREAMBANK  TREES 


APPENDIX  C 


% 

(i 

Method  for  predicting  afflux  due  to  LWD,  developed  by  Gippei  et.  al.  (1992)  • 


The  recommended  procedure  for  predicting  the  hydraulic  effect  of  managing  large  woody  debris 
from  a  lowland  river  is  as  follows: 

1 )  Measure  the  LWD; 

•  projected  length  of  LWD  (L*) 

•  mean  diameter  of  LWD  in  flow  (d) 

«  angle  of  orientation  of  the  LWD  in  the  flow  (a) 

2)  Measure  the  chamiel  morphology: 

•  cross-sectional  area  of  flow  at  selected  discharge  (A) 

3)  Measure  or  estimate  flow  characteristics  at  selected  discharge: 

•  depth  of  flow  downstream  of  LWD  (hi) 

•  'ocity  downstream  of  LWD  (Uj) 

4)  Select  a  drag  coefficient: 

•  based  on  angle  of  orientation  and  snag  form  (C  d)  using  Figure  Al  or  A2 

5)  Calculate  the  following; 

•  Froude  number  downstreain  of  LWD 

F  = 

•  blockage  ratio  of  LWD 

B  =  L.d/A 


JZl 


6) 


7) 

8) 


o  drag  coefficient  corrected  for  blockage 
Cl)  =  C  i)(l-B)'^ 

Calculate  afflux  due  to  1  .WD: 


Ah  =  - 


-\)  +  ^(f^  -if 


Calculate  the  upstream  extent  of  the  afflux  using  a  backwater  procedure 

Repeat  the  calculations  for  various  management  strategies  such  as  lopping  and  rotation 


•  i 


,20 


« 


Figure  1  :  Yarution  in  drag  coefficient  with  angle  of  rotation  to  the  flow,  measured  for  a 
model  LWD  complete  with  trunk,  branches  aeid  butt,  and  for  other  combinations  of  these 

« 

components 


•  Cain,>l»t»  LWb 
o  T(unkandba»ichc» 


Figure  2  :  Variation  in  drag  coeflicient  with  angle  of  rotation  for  cylinders  of  various  lengths 
and  diameters,  lioemer’s  (1958)  relationship  is  for  infinitely  lon^,  cylinders. 


i 
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APPENDIX  D  ;  FIELD  DATA  AND  ANALYSIS  CALCULATIONS 
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